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We compare the rigorous equations describing the motion of spinning test particles in gravitational 
and electromagnetic fields, and show that if the Mathisson-Pirani spin condition holds then exact 
gravito-electromagnetic analogies emerge. These analogies provide a familiar formalism to treat 
gravitational problems, as well as a means for a comparison of the two interactions. Fundamental 
differences are manifest in the symmetries and time projections of the electromagnetic and gravita- 
tional tidal tensors. The physical consequences of the symmetries of the tidal tensors are explored 
comparing the following analogous setups: magnetic dipoles in the field of non-spinning/spinning 
charges, and gyroscopes in the Schwarzschild, Kerr, and Kerr-de Sitter spacetimes. The implications 
of the time-projections of the tidal tensors are illustrated by the work done on the particle in vari- 
ous frames; in particular, a reciprocity is found to exist: in a frame comoving with the particle, the 
electromagnetic (but not the gravitational) field does work on it, causing a variation of its proper 
mass; conversely, for "static observers", a stationary gravitomagnetic (but not a magnetic) field does 
work on the particle, leading to a quantitative explanation of the Hawking- Wald spin interaction 
energy. The issue of hidden momentum, and its counterintuitive dynamical implications, is also 
analyzed. Finally, a number of issues regarding the electromagnetic interaction are clarified, namely 
the differences in the dynamics of electric and magnetic dipoles, and the physical meaning of Dixon's 
equations. 
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I. INTRODUCTION 

The analogies between the equations of motion for gy- 
roscopes in a gravitational field and magnetic dipoles in 
an electromagnetic field have been known for a long time, 
and were presented in many different forms throughout 
the years. This is the case for both the force equation 
(center of mass motion) and the spin evolution equation 
of these test particles in external fields. The former was 
first found by Wald |T] in the framework of linearized the- 
ory, who showed that the gravitational force exerted on a 
spinning pole-dipole test particle (hereafter a gyroscope), 
whose center of mass is at rest in a stationary field, takes 
the form Fq = KW(H-S), where H is the so-called "grav- 
itomagnetic field", K is some constant (depending on the 
precise definition of H, e.g. [2r0]), and S is the particle's 
angular momentum. This formula is similar to the for- 
mula for the electromagnetic force on a magnetic dipole, 
Fem = V(B ■ ft), where B is the magnetic field and Jl 
is the dipole's magnetic moment. The analogy was later 
cast in an exact form by one of the authors in [5] , using 



the exact gravitoelectromagnetic (GEM) inertial fields 
from the so-called 1+3 "quasi-Maxwell" formalism. The 
force was seen therein to consist of an electromagnetic- 
like term in the form above plus a term interpreted as the 
weight of the energy of the gravitomagnetic dipole, and 
the limit of validity of the analogy was extended to arbi- 
trarily strong stationary fields and when the gyroscope's 
worldline is tangent to any time-like Killing vector field 
(which comprehends e.g. circular trajectories with arbi- 
trary speed in axisymmetric spacetimes) . Recently it was 
shown that there is actually an exact, covariant and fully 
general analogy relating the two forces; this analogy is 
made explicit not in the framework of the "gravitoelec- 
tromagnetic" vector fields, but by using instead the tidal 
tensors of both theories, introduced in [B] . 

The analogy between the so-called "precession" of a 
gyroscope in a gravitational field and the precession of a 
magnetic dipole under the action of a magnetic field was 
noticed long ago, in the framework of linearized theory, 
by a number of authors [3] 0] [7J [5] , who pointed out that 
the spin vector of a gyroscope at rest in a stationary field 
evolves as dS/dt — KS x H. This formula is similar 
to the formula for the precession of a magnetic dipole 
in a magnetic field, dS/dt — /2 x B. The analogy was 
later cast in an exact form in the framework of the GEM 
inertial fields, e.g. [HOSEI]; it is not covariant, holding 
only in a specific frame comoving with the particle, but, 
in the formulation in [5] [TO] (as we shall see in detail) , 
the test particle can be moving with arbitrary velocity in 
an arbitrary field. 

These analogies provide a familiar formalism to treat 
otherwise complicated gravitational effects, as well as a 
means to compare the two interactions. In this work we 
explore these analogies, exemplifying their usefulness in 
some applications, and especially the insight they provide 
into fundamental, yet not well know, aspects of both in- 
teractions. 

We will also make use of a third exact gravito- 
electromagnetic analogy (see e.g. [lTMl3| ). this one a 
purely formal one (see [33]), relating the quadratic scalar 
invariants of the Maxwell and Weyl tensors [TT] [E?] Q5] . 
The insight provided by this analogy will be crucial in 
the physical applications of Sec. 



IV B 



A. The equations of motion 

In this work we start, in Sec. [TTJ, by writing the rigor- 
ous, general relativistic equations describing the motion 
of spinning test particles with gravitational and electro- 
magnetic pole-dipole moments, subject to gravitational 
and electromagnetic external fields, in terms of physi- 
cally meaningful quantities. This turns out not to be a 
straightforward task, as the covariant equations for this 
problem are still not generally well understood, with dif- 
ferent methods and derivations leading to different ver- 
sions of the equations, the relationship between them 
not being at all clear. Perhaps more surprising is the 
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fact that it is the electromagnetic (not the gravitational) 
field that has been posing more difficulties, with a num- 
ber of misconceptions arising in the physical interpreta- 
tion of the quantities involved in the equations (some au- 
thors [THl [T7] have even concluded that such a covariant 
description is not possible) . The main difficulties concern 
the interpretation of the time components of the forces 
and the mass of the test particles, and the fact that in 
most works, e.g. [T8Tf2"3"] . equations of motion are pre- 
sented (Eqs. (|A2|-( A3 1) which are symmetric in terms of 



electric and magnetic dipoles; that would, given the dif- 
ferent nature of the two dipole models, seemingly lead (in 
many ways) to physically inconsistent predictions. This 
is discussed in detail in Appendix |Aj where it is shown 
that those equations are correct and consistent with the 
equations used in this work, Eqs. |8|-([9]); the difficulties 
arise in the physical interpretation of the quantities ap- 
pearing therein. 

In order to form a determined system, the equations 
of motion need to be supplemented by a spin condition; 
the latter is even today still regarded as an open ques- 
tion, with a long history of debates concerning which 
one is the "best" condition (see Ref. 24J for a compre- 
hensive review and list of references). In Sec. II A and 
Appendix [C] we briefly discuss its meaning, and the re- 
lated problem of the relativistic definition of center of 
mass; a detailed discussion of the subtleties involved is 
given in the companion paper [25) . This subject is of 
relevance in this work because in order for the two phys- 
ical gravito-electromagnetic analogies mentioned above 
(both for the force and for the precession of the spin vec- 
tor) to be exact, the Mathisson-Pirani spin condition is 
required. This condition is poorly understood in most 
literature; namely, its degeneracy and the exotic heli- 
cal solutions it allows. A lot of skepticism has been 
drawn into it because these helical solutions have been 
deemed unphysical. These assessments arise from mis- 
conceptions dissected in [25], where it is clarified that 
there is nothing wrong with Mathisson-Pirani spin con- 
dition; it is as valid as the Tulczyjew-Dixon [2"2l |21)I 12"7] . 
the Papapetrou-Corinaldesi [28 , or any other physically 
reasonable condition (there are infinite possibilities), the 
choice between them being just a matter of convenience. 
Many applications in this work are examples where the 
Mathisson-Pirani condition is the best choice; that is dis- 
cussed in detail in Appendix [C] 

Also related with the spin supplementary condition is 
an issue central to the understanding of the dynamics 
of a spinning particle: the decoupling of the 4-velocity 
from the 4-momentum, discussed in Sec. |IID| and Ap- 
pendix [D] The 4-velocity of a spinning particle is in gen- 
eral not parallel to its 4-momentum; it is said to possess 
"hidden momentum", which, for a pole-dipole particle, 
consists of two parts: 1) the hidden momentum of elec- 
tromagnetic origin and 2) a pure gauge term originated 
by the choice of spin condition, which we dub "inertial 
hidden momentum"; in the framework of the 1+3 for- 
malism (with the appropriate condition) , it is cast as the 



exact inertial analogue of 1). The hidden momentum is 
known to lead to counter-intuitive behavior of the spin- 
ning particles; examples are the bobbings studied in [2"5] . 
and also (as shown in [S3]) Mathisson's helical motions 
themselves, where a particle accelerates without the ac- 
tion of any force. In this work we present another, per- 
haps even more surprising consequence in Sec. |IV A| due 
to the hidden momentum, a magnetic dipole with radial 
initial velocity in the field of a point charge accelerates 
in approximately opposite direction to the force. 



B. The main realizations 



Most of our applications, Sees. |III||V| of this paper, will 
deal with the tidal tensor formalism introduced in [5] , and 
the exact analogy it unveils: both the electromagnetic 
force on a magnetic dipole and the gravitational force on 
a gyroscope are given by a contraction of a rank 2 mag- 
netic type tidal tensor (B a p, H a| g), with the dipole/spin 
4-vector. Here B a p gives the tidal effects of the mag- 
netic field and H Q| g is the magnetic part of the Riemann 
tensor, both measured in the particle's rest frame. This 
makes this formalism especially suited to compare the 
two forces — it amounts to simply comparing the two 
tidal tensors. Such comparison is done through Einstein 
and Maxwell equations, as they also can be written in 
terms of tidal tensors. Apart from the non-linearity of 
W a 0, the tensorial structure differs essentially when the 
fields vary along the test particle's worldline, since this 
endows B a p with an antisymmetric part (non-vanishing 
even in vacuum), and a non- vanishing time projection 
along that worldline, whereas its gravitational counter- 
part is spatial and, in vacuum, symmetric. This signals 
fundamental differences in the interactions, as the anti- 
symmetric parts of the electromagnetic tidal tensors (as 
well as their time-projections) encode the laws of elec- 
tromagnetic induction. We discuss them in two separate 
sections. 

In Sec. |IV| we explore the physical consequences of the 
different symmetries of the gravitational and electromag- 
netic tidal tensors. They are seen to imply e.g. that par- 
ticles moving in a non-homogeneous electromagnetic field 
always measure a non- vanishing B a p (thus feel a force), 
which is not necessarily the case in gravity. The fol- 
lowing analogous setups are compared: magnetic dipoles 
in the field of non-spinning/spinning charges, and gyro- 
scopes in the Schwarzschild, Kerr, and Kerr-dS space- 
times. It is seen that in the cases where B a p reduces 
to B[ a j3] , we have H Q ^ = in the gravitational analogue 
(i.e., gyroscopes feel no force). Geodesic motions for spin- 
ning particles are found to exist in Schwarzschild (radial 
geodesies) and in Kerr-dS (circular equatorial geodesies) 
spacetimes. 

In Sec. [V]we explore the physical content of the time 
projections of the forces in different frames. These are 
related with the rate of work done on the test particle 
by the external fields; in order to obtain that relation- 
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ship, we start by deriving the general equation yielding 
the variation of energy of a particle with multipole struc- 
ture with respect to an arbitrary congruence of observers. 
This generalizes the power equation in [5] for test par- 
ticles possessing hidden momentum and varying mass, 
and non-spatial forces. We then show that the electro- 
magnetic force has a non-vanishing projection along the 
particle's 4- velocity U a , which is the power transferred 
to the dipole by Faraday's induction, reflected in a vari- 
ation of its proper mass to. The projection of the gravi- 
tational force along U a , by contrast, vanishes (since H Q) g 
is spatial with respect to U a , manifesting the absence of 
an analogous induction effect), leading to the conserva- 
tion of the gyroscope's mass. Also of particular interest 
in this context are the time projections as measured by 
"static observers", analyzed in Sec. |VB| For these ob- 
servers, the time projection of the electromagnetic force 
vanishes, meaning that the total work done on the mag- 
netic dipole is zero. This reflects the well known fact 
that the work done by the magnetic field is zero, which 
in this framework is seen to arise from an exchange of 
energy between three forms: translational kinetic energy, 
proper mass to, and "hidden energy", which occurs in a 
way such that their variations cancel out, keeping the 
total energy constant. In the gravitational case, since 
the induction effects are absent, such cancellation does 
not occur and, unlike its electromagnetic counterpart, a 
stationary field does work on mass currents, there exist- 
ing an associated potential spin-curvature potential en- 
ergy. Considering, as a special case, the gyroscope to 
be a small Kerr black hole, falling axially into a large 
Kerr black hole, it yields Wald's [T] gravitational spin- 
spin interaction energy, which quantitatively explains the 
spin-dependence of Hawking's upper bound [30] for the 
energy released by gravitational radiation when two black 
holes collide. 

It is worth noting that the above mentioned differences 
between the two interactions are not negligible in the 
weak field and slow motion regime, as is shown in detail 
in Sec. Ill (this is commonly overlooked in the literature 



concerning this regime). 



C. Beyond the pole-dipole 

In Sec. |VI| we go beyond the pole-dipole approximation, 
including the moments of quadrupole order, to clarify the 
mechanism by which the proper mass of a spinning par- 
ticle in a electromagnetic field varies, and solve an ap- 
parent contradiction of the former approximation: On 
the one hand (as stated above), we show in Sec. V]that 
the projection of the force on a magnetic dipole along its 
worldline is the rate of work done on it by the induced 
electric field Ei nc i (as measured in the particle's frame), 
which is manifest in a variation of the particle's proper 
mass to. For a quasi-rigid spinning particle, the latter 
is a variation of rotational kinetic energy. Likewise, the 
angular momentum of the particle (and its magnitude 



S = v S 2 ) must also vary, as expected from the fact that 
Eind, by having a curl, should torque the spinning body. 
On the other hand, however, no trace of that torque (fi n d) 
is found to dipole order, where the torque on a magnetic 
dipole is the textbook expression f — fl x B, leading to 
the conservation of S. In Sec. |VI A| this apparent para- 
dox is clarified; indeed a torque T\ n d generated by the 
induced electric field is exerted on the particle (and it 
is governed by the time projection of the magnetic tidal 
tensor B a ^U^), but, it involves second moments of the 
charge (i.e., of quadrupole order), which is the reason 
why it does not show up in pole-dipole approximations; 
however the rate of work it does, T; n( j • Q (and the corre- 
sponding variation of to), can be written in terms of the 
the magnetic dipole moment, yielding exactly the time 
projection of the dipole force along its worldline, as ob- 
tained in Sec. |V A| Then in Sec. |VI B| we study the anal- 
ogous gravitational problem, showing that, as expected 
(as H Q( 3 is spatial), no analogous torque exists. In partic- 
ular it is shown that the gravitational field cannot torque 
a "spherical" body (case in which the electromagnetic 
torque reduces essentially to find), which is seen to lie 
at the heart of the frame dragging effect. 



D. Notation and conventions 

Signature and signs. We use the signature h 

++; £af}o-{ = y/— g[ct/3~fS) denotes the Levi-Civita 
tensor, and we follow the orientation [1230] = 1 
(i.e., in flat spacetime, 61230 = 1). tijk = eyfco is 
the 3-D alternating tensor. We use the convention 



for the Riemann tensor: R a , 



P/J.U 



pa 

1 p V ,y.- 



pa 

1 p»y 



2. Sometimes we use the abbreviation e Q( 9 7 = 
CafijsU 5 , where U a is the 4- velocity of the test par- 
ticle's CM. 

3. We use bold fonts to denote tensors T (including 4- 
vectors v) , and arrows for 3- vectors v. Greek letters 
a, ft, 7, ... denote 4-D spacetime indices, Roman 
letters i,j,k,... denote 3-D spatial indices. Follow- 
ing the usual practice, sometimes we use compo- 
nent notation T a/3 to refer to a tensor T. 

4. By u a we denote a generic unit time-like vector, 
which can be interpreted as the instantaneous 4- 
velocity of some observer. U a = dz a /dr is the tan- 
gent vector to the body's representative worldline 
z a (r), which in this work is taken to be a suitably 
defined center of mass; U a is thus the 4-velocity of 
the test particle's center of mass. 

Time and space projectors. (T")^ = —u a up and 



(h u y 



u a u/3 denote, respectively, the pro- 



jectors parallel and orthogonal to a unit time-like 
vector u a ; may be interpreted as the time and space 
projectors in the local rest frame of an observer of 
4-velocity u a . 
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6. Tensors resulting from a measurement process. 
(A") 01 " " denotes the tensor A as measured by 
an observer 0{u) of 4-velocity u a . For example, 
(E u ) a = F«uP, (E u ) a/ 3 ee Faj.pu-1 and (E u ) a/3 = 
Ravpi/U^u 11 denote, respectively, the electric field, 
electric tidal tensor, and gravito-electric tidal ten- 
sor as measured by 0{u). Analogous forms apply 
to their magnetic/gravitomagnetic counterparts. 
For 3- vectors we use notation A(u); for example, 
E(u) denotes the electric 3- vector field as measured 
by 0(u) (i.e., the space part of (E u ) a , written in a 
frame where u 1 = 0). When u a = U a (i.e., the par- 
ticle's CM 4-velocity) we drop the superscript (e.g. 
(E u ) a ee E a ), or the argument of the 3-vector: 
E(U) ee E. 

7. Electromagnetic field. The Maxwell tensor F a/3 and 
its dual *F a ° decompose in terms of the electric 
(E u ) a ee F^u? and magnetic (B u ) a = *F a p u^ 
fields measured by an observer of 4-velocity u a as 

F aP = 2u [a {E u ) 0] +e aPlS u 5 {B u Y ■ (1) 
*F a/3 = 2u [a (B u ) p] ~e a ^ a u' 7 {E u y . (2) 

8. Static observers. In stationary, asymptotically flat 
spacetimes, we dub "static observers" the rigid con- 
gruence of observers whose worldlines are tangent 
to the temporal Killing vector field £ = d/dt; 
may be interpreted as the set of points rigidly 
fixed to the "distant stars" (the asymptotic iner- 
tial rest frame of the source). For the case of 
Kerr spacetime, these correspond to the observers 
of zero 3-velocity in Boyer-Lindquist coordinates. 
This agrees with the convention in e.g. [311 132] . 
Note however that the denomination "static ob- 
servers" is employed with a different meaning in 
some literature, e.g. [33], where it designates rigid, 
vorticity-free congruences (existing only in static 
spacetimes). In the case of the electromagnetic sys- 
tems in flat spacetimes, by static observers we mean 
the globally inertial rest frame of the sources. 

9. GEM. This is the acronym for "gravitoelectromag- 
netism". By "inertial GEM fields", we mean the 
fields of inertial forces that arise from the 1+3 split- 
ting of spacetime: the gravitoelectric field G, which 
plays in this framework a role analogous to the 
electric field of electromagnetism, and the gravit- 
omagnetic field H, analogous to the magnetic field. 
These fields are discussed in detail in the compan- 
ion paper [3~4"j . 

II. EQUATIONS OF MOTION FOR SPINNING 
PARTICLES IN A CURVED AND 
ELECTROMAGNETIC BACKGROUND 

In most of this work we will be dealing with the dy- 
namics of the so-called pole-dipole spinning test parti- 



cles. We consider systems composed of a test body plus 
background gravitational and electromagnetic fields. Let 
(T tot ) Q ^ = + (T matter ) a P denote the total energy 
momentum tensor, which splits into the electromagnetic 
stress-energy tensor ® a P and the energy-momentum of 
the matter {T matter ) afi . Moreover, let T a ? and j a de- 
note, respectively, the energy momentum tensor and the 
current density 4-vector of the test body. We also con- 
sider that the only matter and currents present are the 
ones arising from the test body: (T matter ) a P = T a/3 , 
(jtot) Q = j a - In this case (see Appendix [D]for details) the 
conservation of total energy-momentum tensor yields: 

(T tot ) a % = T a % = -G a % & T a % = F^jp , 

(3) 

where F a P is the Maxwell tensor of the external (back- 
ground) electromagnetic field. 

In a multipolc expansion the body is represented by 
the moments of j a (its "electromagnetic skeleton") and 
a set of moments of T Q ^, called "inertial" or "gravita- 
tional" moments (forming the so called [35] "gravitational 
skeleton"). Truncating the expansion at dipole order, the 
equations of motion for such a particle involve only two 
moments of T a/3 : 

P a ee / T aP dSp , (4) 

S af> ee 2 f r^T^d^ , (5) 
and two electromagnetic moments |21) : 

q ee f j a dX a , (6) 

Q a P ee / (r^/'t/ 7 -r la U 0] p)dE 7 . (7) 

Here the motion of the test particle is described by 
a reference worldline z a (r) relative to which the mo- 
ments are taken, where r is the proper time. The pre- 
cise choice of this worldline will be discussed below. 
U a ee dz a /dr denotes its (unit) tangent vector; P a {r) 
is the 4-momentum of the test particle; S oi3 (t) is the an- 
gular momentum about the point z a (r) of the reference 
worldline; r Q (r) ee x a — z a (r), where {x a } is a chart 
on spacetime; E(r, U) = S(z(t),J7) is the spacelike hy- 
persurface generated by all geodesies orthogonal to U a 
at the point z Q (r); dS 7 ee n 7 dS, where n 7 is the (past- 
pointing) unit normal to S(r, u) and c?S is the 3- volume 
element on E(r, u); q denotes the total charge, which is 
an invariant (making E in this case arbitrary); and Q a ^ 
is the electromagnetic dipole moment about the point 
z q (t). Expressions Q, ^ and ^ are integrals of rank 
2 and 3 tensors over E (i.e., amount to adding tensor 
components at different points), which requires a justi- 
fication: the pole-dipole approximation amounts to con- 
sidering T a P and j a non-vanishing only in a very small 
region around z a (r), so that only terms linear in r are 
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kept; since to first order spacetime can always be taken 
as flat, these integrals are meaningful mathematical op- 
erations and do indeed define tensors, just like in flat 
spacetime. 

The equations of motion are [2TJ H7J OH] ( see Ap- 
pendix [A] for a discussion) : 



DP a 
dr 



DS a(3 

dr 



1 

2 " r " p " 2' 
Dd? 



(8) 



= 2P^ a U p] + 2^F% + 2d [a F p \W , (9) 



where F a/3 is the background Maxwell tensor, d a and \i a p 
are the time and space projections of Q a @ with respect 
to U a : 

d a = _Q^ U0 . = {h u ) a ^h u f 5 Q~< 5 , (10) 

where {h u ) l3 a = S@ + U^U a . Introducing also \i a by 

= e a ^U 5 ; ii a = ^Vr^V 5 , (H) 

the vectors d a and [i a have the form (0,d) and (0,/2) in 
the frame where U l — 0, and correspond to the usual 
notion of electric and magnetic dipole moments with re- 
spect to an observer of 4- velocity U a . With respect to 
an arbitrary observer of 4- velocity u a we can decompose: 



Q"' 9 = 2{d u ) [a u^ +e a ^ s {^ u ) 1 u s , 



(12) 



(d u ) a and {^ u ) a = \^ a p 1 ^ i being the electric and 
magnetic dipole moments as measured by u a . In some 
applications we will assume fi a @ to be proportional to the 
spin tensor: fi a/3 — aS al3 , where a is the gyromagnetic 
ratio. 

The first term in ^ is the Lorentz force; the sec- 
ond term, \F llu '' a [i^ v = F£ M , is the force due to the 
tidal coupling of the electromagnetic field to the mag- 
netic dipole moment; and the third, — \FL OL [ illv S tiV = Fg 
is the Mathisson-Papapetrou spin-curvature force. The 
last two terms are the force exerted on the electric dipole, 



consisting of a tidal term F°L.gU 1 d^ and of a 



tidal 



term F^Dd^/dr. 

Up until now, the reference worldline z a (r), relative to 
which the moments in Eqs. (|8|-([9| above are taken, is 
still undefined. Had we made an exact expansion keep- 
ing all the infinite multipole moments as in [371 13H] , such 
worldline would be arbitrary. Herein, however, it must 
be assumed that it passes through the body (or close 
enough), so that the pole-dipole approximation is valid; 
it will be chosen as being prescribed by a suitably defined 
center of mass of the test particle. As discussed in the 
next section, that is done through a supplementary con- 
dition S a ^up = 0, for some time-like unit vector field u a . 
If F a P = there are 13 unknowns in Eqs. (l8|)-(|9j) {P a , 3 



independent components of U a , and 6 independent com- 
ponents of S a P) for only 10 equations. The condition 
S a "up = 0, for a definite u a , closes the system as it 
kills off 3 components of the angular momentum. But 
in the general case F a/3 ^ one needs also to give the 
laws of evolution of Q a " in order for the system to be 
determined, cf. [55] . 



Center of mass and spin supplementary 
condition 



In relativistic physics, the center of mass of a spinning 
particle is observer 1 dependent. This is illustrated in Fig. 
1 of [25]. Thus one needs to specify the frame in which 
the center of mass is to be evaluated. That amounts 
to supplement Eqs. (fS >— dsl) (which, as discussed above, 
would otherwise be undetermined) by the spin subsidiary 
condition S a ^up — 0, as we will show next. The vector 



(d u G . 



U 6 



(13) 



yields the "mass dipole moment" (i.e. the mass times 
the displacement of the reference worldline relative to 
the center of mass) as measured by the observer O of 
4- velocity u a . In order to see this let us, for simplicity, 
take F & = so that the integral ^ does not depend 
(to pole-dipole order) on E, by virtue of the conservation 
equations T™^ = 0. And let us set up, at the reference 
worldline (cf. |40j ) . a locally nearly Lorentz frame e a mo- 
mentarily comoving with 0; i.e., eo = u a and the triad 
e~i spans the instantaneous rest space of O. In this frame, 
u" = and S % ^u^ = S zQ Uq; also r a = x a , and thus from 
Eq. (§ we have: 



S M = 2 



S(t,u) 



x'T M d i x = m{u)x l CM {u)\ 
(14) 

note that r° = 0, since the integration is performed 
in the hypersurface £(t, it) orthogonal to u a . Here 
m(u) = J T 00 d 3 x = —P a u a denotes the mass as mea- 
sured by observer O. We see that S M is by definition the 
mass dipole in the frame e a : S M — m(u)x l CM (u) = (rf^) 1 , 



1 We will make reference to quantities such as the center of mass, 
mass dipole, etc as measured by a given observer 0(u) with some 
4-velocity u a . This demands an explanation. These quantities 
involve integrals over a spacelike 3-D hypersurface E; E is defined 
as the hypersurface generated by all geodesies orthogonal to the 
observer's worldline; in flat spacetime this corresponds to the 
notion of rest space of 0(u), defined as the 3-space orthogonal 
to the vector field obtained by parallel transporting u a to every 
point. In any case it is the 4-velocity u a , at the given point, that 
entirely defines the hypersurface of integration, hence we refer to 
this as the quantities measured by a given observer of 4-velocity 
u a , in place of quantities measured in a given frame (e.g. the 
rest frame of u a ). 
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and x l CM (u) — S M /m{u) is the center of mass position. 
Thus the condition S a/3 up = 0, implying in this frame 

S M = =>■ x l CM (u) — 0, is precisely the condition that 
the reference worldline is the center of mass as measured 
in this frame (or, equivalently, that the "mass dipole mo- 
ment" vanishes for such an observer). The condition 
S a Pup = allows us to write S af} = e aP jg S^u s , where 
S a is the spin 4- vector, defined as the 4- vector with com- 
ponents (0, S) in the rest frame of O (where S is the 
angular momentum vector measured by this observer). 
For details on how the center of mass position changes in 
a change of observer, we refer the reader to [25] . 

In general one wants the equations of motion not to 
depend on quantities (the center of mass) measured by 
a particular observer, but instead on a center of mass 
defined only in terms of "intrinsic" properties of the par- 
ticle. Two conditions for accomplishing this are com- 
mon in the literature 2 : the Frenkel-Mathisson-Pirani 
condition S a ^Ufj — (hereafter Pirani's condition, as 
it is best known) and the Moller-Tulczyjew-Dixon con- 
dition S a PPp — (hereafter Dixon's condition). The 
former, which is analogous to the magnetic side con- 



dition n a,3 Ui3 = (if we assume 



a/3 _ 



aS a ?, they 



are the same), is the one on which the exact gravito- 
electromagnetic analogies presented in this work rely. It 
seems also the most natural choice, as it amounts to com- 
puting the center of mass in its proper frame, i.e., in the 
frame where it has zero 3-velocity. Such center of mass 
is dubbed a "proper center of mass" [IT]. This condition 
arises also in a natural fashion in a number of treatments 
[121 143] (see also [44"]): it turns out, however, that it does 
not determine the worldline uniquely. For instance, in 
the special case of Minkowski spacetime in the absence 
of electromagnetic fields, Eqs. ([8|-([9]), supplemented by 
S a "Up = 0, are known to lead, in addition to the ex- 
pected straightline motion, to an infinite set of helical 
motions, which were notably found by Mathisson |4"5] . 
These solutions have been deemed unphysical by some 
authors, with the argument that the radius of the helices 
can be arbitrarily large (see e.g. [371 13S1 HS1 EH), which 
would be contradicted by experiment. As shown in detail 
in |25j . those statements originate from a misconception 
in the derivation of the radius. The reason for the in- 
determinacy of this condition is seen from the fact that 
for a spinning particle not only arc there infinite possi- 
ble centers of mass, depending on the observer, but there 
are also infinite proper centers of mass |41j ; the motion of 
the latter yielding the helices. Their radius is finite and 
contained within a tube of radius i? m ax = S/M, which is 
actually the minimum size that a classical spinning par- 
ticle can have, if it is to have finite S and positive mass 
without violating the speed of light limit [J5] (see also 
[UJ 0H] for a different argument). The different world- 



2 A review (with a comprehensive list of references) on the litera- 
ture regarding this subject may be found in |24| . 



lines compatible with this condition are just equivalent 
descriptions of the same motion. 

The Dixon condition S a ^Pp = amounts to comput- 
ing the CM in the frame where it has zero 3-momentum. 
This condition determines uniquely the center of mass 
worldline [221 HOI [5TJ ; there is no ambiguity in this case, 
since P a is given in advance by ([4]) , and for this reason is 
preferred by a number of authors. Note that in flat space- 
time the worldline specified by S a/3 Pp = corresponds to 
Mathisson's non-oscillatory solution; but in the presence 
of gravitational/electromagnetic field, P a cannot in gen- 
eral parallel to U a if we assume S a/3 Pp = (see Eq. (35) 
of [25] for a general expression for U a in terms of P a 
and S a °). In the presence of an electromagnetic field 
it presents some counter-intuitive properties as well (see 
Appendix C 2 1 . 



The fundamental point that we want to emphasize here 
is that (contrary to the point of view in some literature) 
these two conditions, as well as others explored in [10] 
(including non-covariant ones such as the Corinaldesi- 
Papapetrou condition |28j). are equivalent descriptions 
of the motion of the test particle, the choice between 
them being a matter of convenience. This is one of the 
main points made in |25j . It is clear that for a free par- 
ticle in flat spacetime the Dixon condition gives the sim- 
plest, non-degenerate, description (straightline motion); 
but in some more complicated systems the Pirani condi- 
tion turns to be more suitable (examples are the appli- 
cations in Sees. |IV A[ |IVB"1 and [V B[ this is discussed in 
detail in Appendix \G\. In the case of massless particles, 
it has been argued in [521 [53] that the latter condition is 
actually the only one that can be applied. 

In the present work we will be mainly using Pirani's 
condition, since, as we shall see, it is the one on which 
the exact gravito-electromagnetic analogies rely. The dif- 
ferences between the two conditions, as well as their im- 
pact on the gravito-electromagnetic analogies, are dis- 
cussed in Appendix [C] With this choice, we have S^" — 
e^ TX S T U\, where S a is the spin 4-vector, defined as be- 
ing the 4-vector with components (0, S) in the CM frame. 
Substituting in Eqs. Q-Q (and also performing the con- 
tractions with U a ) we obtain 



DP a 
oIt 



dr 



e^uU v {d a Ef } +/AE^) , (16) 



where Dp/dr denotes the Fermi-Walker derivative, E a = 
F af3 Ufs and B a = *F a ^U are the electric and mag- 
netic fields as measured by the test particle, and E a p = 



F a ^U^ B af3 = *F afl;0 U» and 



are, respectively, the electric, magnetic and gravitomag- 
netic tidal tensors as defined in 6, 34], measured by the 
test particle. 



B. Force on gyroscope vs. force on magnetic dipole 
- exact analogy based on tidal tensors 

Herein we are interested in purely magnetic dipoles, 
i.e., dipoles whose electric moment vanishes in the CM 
frame; this is ensured by the condition d a = =>■ Q a ° = 
li ap ~ , see expressions (10). In this case Eq. (15) simplifies 
to: 



DP° 
~dV 



= qF a PU p + B pa n, 3 



0a, 



(17) 



These equations manifest the physical analogy B a p 
M a ff, summarized in Table llj (i) both the electromag- 
netic force on a magnetic dipole and the gravitational 
force on a gyroscope are determined by a contraction 
of the spin/magnetic dipole 4- vector with a magnetic 
type tidal tensor. B a p is a covariant derivative, keep- 
ing U a fixed (covariantly constant), of the magnetic field 
B a = -kF a pUP measured by the test particle: B a p — 
B a ;0\u=const.', i.e., it is a derivative of B a as measured in 
the inertial frame momentarily comoving with the parti- 
cle. For this reason it is dubbed B a p the magnetic tidal 
tensor, and its gravitational counterpart H Q| a the grav- 
itomagnetic tidal tensor [6]. (ii) It turns out that B a p 
and H Q; 3 obey the formally similar equations (Tl|2) and 
(|l}3) in Table [!) which in one case are Maxwell equa- 
tions, and in the other are exact Einstein equations. That 
is: the traces (|l]2) are, respectively, the time projection 



of the electromagnetic Bianchi identity *F' 



Ct0 



and 



the time-time projection of the algebraic Bianchi iden- 
tities *i? 7Q 7j g = 0; the antisymmetric parts (jlj3a) are, 
respectively, the space projection of Maxwell's equations 



F"^ = 47T? Q and the time-space projection of Einstein's 
equations i? p „ = 87r(T M1 , — \g^vT a a ). The electromag- 
netic equations take a familiar form in an inertial frame: 
Eq. §2a) becomes V • B = 0; the space part of §3a) 
is the Maxwell- Ampere law V x B = dE/dt + Airj. The 
latter means that the space part of B\ a m encodes the 
curl of B a , which is actually a more general statement, 
holding in arbitrarily accelerated frames; in such frames 
[assumed non-rotating and non-shearing, cf . Eq. ( 93 ) be- 
low] U a -/3 — —a a Up, and we have: 



^IsB^U 6 = e P l s B T ^ => <?*>B jk = (V x B)\ 

(18) 

Expressing also the second member of |lj3a) in terms of 
the electric and magnetic fields E a and B a measured in 
this frame, we obtain, in 3-vector notation, 

DE - 

VxB = a x E + 4irj 

dr 

which is the generalization of Maxwell-Ampere law for 
accelerated frames (cf. Eq. (19) of [13]). This equation, 
as well as Eq. ( 21 ) below, is of use if one wishes to express 
Eqs. (|lj3a) and ( 19 \) in terms of the fields measured in 
the particle's CM frame (as it in general accelerates). 



According to Table [TJ both in the case of the electro- 
magnetic force on a magnetic dipole, and the case of the 
gravitational force on a gyroscope, it is the magnetic tidal 
tensor, as seen by the test particle of 4- velocity U a , that 
determines the force exerted upon it. The explicit anal- 
ogy in Table [I] is thus ideally suited to compare the two 
forces, because in this framework it amounts to compar- 
ing B a p to H Q ^. The most important differences between 
them are: i) B a p is linear in the electromagnetic poten- 
tials and vector fields, whereas H a/ 3 is not linear in the 
metric tensor, nor in the GEM "vector" fields; for a de- 
tailed discussion of this aspect, we refer the reader to the 
companion paper |34j ; ii) in vacuum, W\ a m = (symmet- 
ric tensor), whereas -B[ a /3] = \ *F a i3. n U 1 ^ (generically 
not symmetric, even in vacuum); iii) time components: 
M. a is a spatial tensor, whereas B a p is not. These last 
two differences, which are clear from equations (|TJ 3 ) - ([TJ 4 ) , 
are the ones in which we are mostly interested in the 
present work. In what follows we explore some physical 
consequences. 

Note this important aspect of Eq. ([l]3a), considering 
for simplicity the vacuum case j a = 0: it tells us that 
when the field F a p varies along the worldline of the ob- 
server U a , that endows B a p with an antisymmetric part, 
meaning also that the tensor itself is non- vanishing. Now, 
since B a p in the force (|ljla) is the magnetic tidal tensor 
as measured by the particle (i.e., U a is the test particle's 
4- velocity), this means that whenever the particle moves 
in a non-homogencous field, a force will be exerted on it 
(except possibly for very special orientations of fl). In 
the inertial frame momentarily comoving with the par- 
ticle, this can be interpreted as the varying electric field 
measured in this frame generating, via the induction law 
VxB = dE/dt, a curl in the magnetic field B, implying 
that the particle sees a non- vanishing magnetic tidal ten- 
sor. The fact that its gravitomagnetic counterpart H Q ^ 
is symmetric in vacuum tells us that no analogous induc- 
tion phenomenon occurs in gravity; this observation is 
crucial in some applications in this work (see Sec. IV). 

There is an electric counterpart to this analogy, relat- 
ing the electric tidal tensor E a p with the electric part of 
the Riemann tensor: 

E a p = F aiM -pU^ < — > E Q( 3 = RaupvU^U" , 

which is manifest in the worldline deviations of both 
theories, see [6], and together they form the gravito- 
electromagnetic analogy based on tidal tensors [SI |3"1] . 
These tensors obey the following equations, which will 
be useful in this work: 



E [<*f3] = ^Fup.^U 1 (a); 



E M ] = (b) (19) 



Eq. ( 19 1) results from the space projection of the identity 



*F ap o = 0; and Eq. (19d) from the time-space projec- 
tion of the identity ★it 7T ^ ) ,^ = 0. Contracting (19i) with 
the spatial 3-form e a p 1 sU 5 yields Eq. |l]4a) of Table I 
Again, for a (non-rotating and non-shearing) arbitrarily 
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Table I: Analogy between the electromagnetic force on a magnetic dipole and the gravitational force on a gyroscope 



Electromagnetic Force 


Gravitational Force 


on a Magnetic Dipole 


on a Spinning Particle 



F§ 



D P a 



B a « 



Eqs. Magnetic Tidal Tensor 

B a a = |I 

IP - 2-Ke aPal fU~< (I 



Fq = -WjS a ; H" p = *R a ^ v U»U» (|T]lb) 
Eqs. Gravito-magnetic Tidal Tensor 



B[aP] — 2* F, 



B aP U a = 0; 



2a) 
3a) 







(I 



y.p*\ 



= — 47T6 



(i 



U S <[T|4a) 



Had 



Had 



2b) 
3b) 

(04b) 



accelerated frame we have: 



or, equivalently: 



z l S E lP U s = e^E^U 5 =* e^E Jk = (V x E)\ 

(20) 

Expressing also the second member of ( 19 \) in terms of 
the fields E a and B a measured in this frame, we obtain, 
in 3-vector notation: 



V x E 



DB 
~dV 



a x E, 



(21) 



which is a generalization of Maxwell- Faraday equation 
V x E = —dB/dt for accelerated frames, cf. Eq. (20) 

of pg. 

The fact that the gravitoelectric tidal tensor E Q( g is 
symmetric again means that there is no analogous gravi- 
tational induction effect, and this is a key observation for 



the applications in Sees. V A and VI 



Eqs. (|Tj4) are the time projections of the tidal ten- 
sors with respect to the observer U a measuring them 
(more precisely, since U a is the particle's 4-velocity, it 
is the time projection in the CM frame); they are zero 
in the gravitational Hq,^ is a spatial tensor, and 

non-zero in the electromagnetic case, which again is re- 
lated to electromagnetic induction, as the right Eq. (|lj4a) 
corresponds to the spatially projected Eq. (19i). This 
means that Fq is spatial with respect to U a , and F£ M is 
not, which has important implications on the work done 
by the fields on the particle, which will be discussed in 
Sec. El 

Finally, note that F£ M = B^ a [i^ is the covariant 
generalization of the familiar textbook 3-D expression 
-Fem = V(/J • B), the latter being valid only in the in- 
ertial frame momentarily comoving with the particle. 



C. Spin "precession" - exact analogy based on 
inertial GEM fields from the 1+3 formalism 



For purely magnetic dipoles, Eq. ( 16 ), for the spin evo- 



lution under Mathisson-Pirani condition, simplifies to: 



■PfS/i 
dr 



(22) 



dr 



(23) 



where a a = DU a /dr denotes the acceleration, and B 13 is 
the magnetic field as measured by the test particle. The 
first term in ( 23 ) embodies the Thomas precession; these 



equations tell us that, in the absence of electromagnetic 
field, S a undergoes Fermi- Walker transport. The second 
term is a covariant form for the familiar torque r = /2 x 
B causing the Larmor precession of a magnetic dipole 
under the influence of a magnetic field; such precession 
is generically described by Eq. p2| . 

Consider now an orthonormal frame e& carried by an 
observer of 4-velocity U a , such that U = eg, comoving 

with the test particle. In such frame, 5° = and U a — 
S?, and equation (23 1 reduces to (see [3"4"]): 



DS* _ 5 .| dS l 



S x £1 + fl x B 



(24) 



where is angular velocity of rotation of the spatial axis 
relative to the tetrad Fermi- Walker transported along 
the center of mass worldline. Note the analogy between 
the two terms of the second equation. We cast the anal- 
ogy between the evolution of the spin vector of a gyro- 
scope and a magnetic dipole subject to a magnetic field 
i n t hree levels; the similarity between t he t wo terms of 
Jiil is the first level. If B l = 0, what Eq. Jill in fact tells 



is that S does not precess relative to the local spacetime 
geometry, since DS l /dr = 0. Thus, the so-called "pre- 
cession" of the gyroscope is an artifact of the reference 
frame, with no local physical meaning; it is simply the 
spatial rotation of the chosen frame relative to the local 
spacetime geometry, or the local "compass of inertia" [5] 
(note that gyroscopes are objects that oppose changes in 
direction; hence torque-free gyroscopes are supposed to 
follow the local compass of inertia). That does not mean, 
however, that it is necessarily meaningless, as it can tell 
us about non-local physical effects. Up until now, the 
spatial angular velocity of rotation of the tetrad, Q, is 
arbitrary since the only condition we imposed was eg to 
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be parallel to the 4-velocity of a given observer. Of spe- 
cial interest for the study of the frame dragging effect is 
the case where the tetrad field is locked to the so-called 
"frame of the distant stars". Such frame, for a station- 
ary, asymptotically flat spacetime, is set up as follows. 
Consider the congruence of static observers (cf. point 
H of Sec. ~ 



I D I : the rigid congruence of observers whose 



worldlines are tangent to the time-like Killing vector field 
that, at infinity, corresponds to observers at rest in the 
asymptotic inertial rest frame of the source. These ob- 
servers are interpreted as being "at rest" with respect to 
the distant stars; since the congruence is rigid, it may be 
thought of as a grid of points rigidly fixed with respect 
to them. This fixes the times axis of the local tetrads 
of the frame. Now one needs to ensure that each local 
spatial triad is locked to this grid; that is done by 
demanding the rotation 51 of the local tetrads (relative 
to Fermi-Walker transport) equals the vorticity of the 
congruence (this defines the frame adapted to the con- 
gruence [541 155| ). In this way the tetrads point towards 
fixed neighboring observers, see [33] for details; and the 
local tetrads themselves are locked to one another and to 
the tetrads adapted to inertial frames at infinity, whose 
axes define star-fixed directions. Hence f2, which is minus 
the angular velocity of precession of a gyroscope relative 
to such frame, yields minus its precession rate relative to 
the distant stars, or to the inertial frame at infinity. It is 
this construction that allows one to mathematically de- 
fine the "precession" of a gyroscope relative to the distant 
stars (i.e., to determine the rotation, relative to the local 
guiding gyroscope, of the frame whose axis are locked to 
the distant stars), which was the effect measured by the 
Gravity Probe B mission [SB]. This is an effect which is 
physical, clearly distinguishing, e.g., the Kerr from the 
Schwarzschild spacetimes, even though it has no local 
significance and can only be measured by locking to the 
distant starts, by means of telescopes (notice that the 
choice of tetrad above is precisely equivalent to requiring 
that static observers see the distant stars fixed in their 
skies) . 

When fl equals the vorticity of the congruence, then 
d = H/2, where H is the "gravitomagnetic field"; this 
designation is due 3 to the fact that it plays in the exact 
geodesic equations (see [34]) the same role as the mag- 



3 In the framework of the GEM fields from the 1+3 formalism, the 
force on the gyroscope [5] 134] and the equation for the geodesic 
motion of a (non-spinning) test particle (see e.g. 5, 9 ) can be ex- 
actly described by equations analogous to the ones from electro- 
magnetism. These results are not however as general as Eq. «24B : 
in the case of the force on the gyroscope, the analogy holds only 
for stationary fields, and if the gyroscope it at rest with respect 
to some stationary observer (i.e., if the gyroscope's worldline is 
tangent to a time-like Killing vector); in the case of geodesies, 
the analogy holds only (again) for stationary fields and station- 
ary observers (i.e., the observer's worldline must be tangent to 
a time-like Killing vector). This is discussed in detail in 134) . 



netic field B in the the electromagnetic Lorentz force. 
This is the second level of the analogy in Eq. (24 1. 



The third level of the analogy can be read from the 
field equations for H, which exhibit striking similarities 
with Maxwell equations for B in an accelerated, rotat- 
ing frame; for a detailed account of this analogy we refer 
the reader to the companion paper |34) . A well known 
manifestation is the correspondence between the Lense- 
Thirring effect and electromagnetism: the gravitomag- 
netic field produced by a spinning mass, as measured by 
the congruence of static observers, is similar to the mag- 
netic field produced by a spinning charge. 

If we assume jl = aS, then the quantity S 2 = S a S a — 
S a "S a p/2 is a constant of the motion, which is immedi- 
ately seen contracting (23) with . This result might 



be puzzling, as it apparently contradicts our assessments 
in Sees. HE and [vj where we conclude that dm/dr is 
a variation of kinetic energy of rotation, via electromag- 
netic induction. If the kinetic energy of rotation varies, 



S must also vary. However, as explained in Sec. [VI] such 
contradiction is but an artifact of the pole-dipole approx- 
imation, which takes into account the work done by the 
induced electric field on the test particle (which involves 
its magnetic dipole moment /I), but neglects the torque 
it exerts (which depends on the particle's second moment 
of the charge, that is of quadrupole order). This also ex- 
plains why the analogy manifest in Eq. ( 24 ) is valid for ar- 



bitrary fields, while most gravito-electromagnetic analo- 
gies based on GEM vector fields (not tidal tensors!), such 
as the geodesic equation/Lorentz force, e.g. jSJH], or the 
forces on gyroscopes/magnetic dipoles [SJI33], break down 
when one considers time-dependent fields [34], [57] (an- 
other exception is the hidden momentum analogy, pre- 
sented in the next section). 



D. Momentum of the spinning particle - hidden 
momentum and exact analogy based on inertial 
GEM fields from the 1+3 formalism 

The momentum Q of a spinning particle is not in 
general parallel to its center of mass 4-velocity U a ; one 
may obtain an expression for P a contracting (p5j) with , 
yielding 



P a = mU a 



DSa p 
dr 



U (j + n^E. 



(25) 



where m = — XJ a P a and E a = F%U^ is the electric field 
as measured by the test particle of 4-velocity U a . We 
split P a in its projections parallel and orthogonal to the 
CM 4-velocity U a : 



pa po 



Pi 



pa _ Tja pa 
-'kin — " lu ' Mii 



hi.:- ' kin — ""' - ' hid — (^ C/ ) Q fl-P' 3 - 

(26) 

We dub the parallel projection P£ in = mU a "kinetic mo- 
mentum" associated with the motion of the center of 
mass. This is the most familiar part of P a , formally 
similar to the momentum of a monopole particle. The 
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component P5L orthogonal to U a is the so-called "hid- 
den momentum" [29]. The reason for the latter denomi- 
nation is seen taking the perspective of an observer 0(U) 
comoving with the particle: in the frame of 0(U) (i.e., 
the frame U l — 0) the 3-momentum is in general not 
zero: P = Phid 7^ 0; however, by definition, the particle's 
CM is at rest in that frame; hence this momentum must 
be somehow hidden in the spinning particle. 

With the Mathisson-Pirani condition S a0 Uf} = 0, the 



momentum P a in Eq. ( 25 ) becomes 



P k « n = mU a ; P h Q id = -e^S^a'U 5 + e a ^E^U 5 . 

(27) 

It is useful to introduce the following notation for the two 
parts of PQ lA : 



pa _ a ,,/9 

MiidEM = e 0i ' 



e^E^U 5 ; P£ 



hidl 



/3 7 <5 



S^a'U 8 (28) 



which we shall now explain. Eq. (27) tells us that 



Pfi ld ^ if the particle is accelerating and/or subject 
to an electric field. In the particle's CM frame (where 
U l = 0), and in vector notation, its space part reads 
OP h °id = 0): 

P hid = P=-Sxa + fIxE = SxG + flxE 1 (29) 

where G = —a is the gravito-electric field as defined in 
El [34] . The term PhidEM — ft x E is a still not well 
known feature of relativistic electrodynamics (despite its 
discovery by Schockley & James [55] dating back from 
the 60's, and having since been discussed in number of 
papers, e.g. [58H6T] . including recent ones [29l [62] ) . It is 
a form of hidden momentum that is originated by the ac- 
tion of the electromagnetic field; but it is however purely 
mechanical in nature, as explained in Appendix[D]using a 
simple model. PhidEM equals minus the electromagnetic 
momentum P x generated by the dipole when placed in 
an external electromagnetic field, which, in the particle's 
frame, reads (see e.g. [SUES] and Appendix |D| : 



Pv 



E x B, 



dipolc 



-fix E 



-Pi 



hidEM- 



This means that in flat space-time, for example, the 
spatial momentum of a dipole whose center of mass is at 
rest is in general not zero in the presence of an electro- 
magnetic field; it is actually PhidEM which allows for the 
total spatial momentum P tot = Pmattcr + Pem to van- 
ish in this stationary configuration, as required by the 



\a/3 



= 0. 



conservation equations (Ttot) .« 

G is a field of "inertial forces"; hence the term Phidi = 
S x G is the "inertial" analogue of fix E (for this reason 
hereafter we will dub it "inertial hidden momentum"), 
and the formal similarity between the two can be un- 
derstood through an analogue model, as discussed in de- 
tail in Appendix [DJ However, an important difference 
between these two forms of hidden momentum must be 
noted: whereas P"idEM ^ s a physical, gauge invariant 



quantity, P£ idI is pure gauge, and depends only on the 
frame where the center of mass is computed. That is 
what we are going to show next. 

Consider a generic spin condition S a ^up = 0, where u a 
denotes the 4- velocity of an arbitrary observer 0(u). As 
discussed in Sec. |II A[ this condition means that we take 
as reference worldline the center of mass as measured by 
0(u). For simplicity, let us consider F a ^ = (so that 
Pj°[ d = Ph" idI ); in this case, contracting ^ with up, and 
using S a *up = 0, leads to [25] 



or 



m(u)U° 



(30) 



where j(U,u) = —U a u a and m(u) = —P a u a . Hence, 
if Du a /dr ^ 0, then in general the momentum is not 
parallel to the 4-velocity: P^ a U^' ^ 0, and thus the par- 
ticle will have hidden momentum. This comes as a nat- 



ural consequence of what we discussed in Sec. II A (see 
also [S3] for further details) : the position of the CM of a 
spinning body depends on the vector u a relative to which 
it is computed; if that vector varies along the reference 
worldline it is clear that this is reflected in the veloc- 
ity U a of the CM, which in general will accelerate even 
without the action of any forces. Also U a will in gen- 
eral no longer be parallel to P a , and thus the centroid is 
not at rest in the frame P % = 0. But if we take a field 
u a such that Du a /dr = 0, i.e., if we take as reference 
worldline the center of mass as measured with respect 
to a field u a that is parallel transported along it, then 
P a = mU a = P^ in , and there is no hidden momentum. 

Thus we see that indeed Ph idI is pure gauge (it is worth 
noting that the condition P a \\ U a is actually cast in [30] 
as one of the possible spin supplementary conditions). 
This also means that the motion effects induced by it 
(such as the bobbings studied in [2S]) must be confined 
to the worldtube of centroids (i.e. the worldtube formed 
by all the possible positions of the center of mass, as mea- 
sured by every possible observer; for flat spacetime, it is 
a tube of radius S/M centered around the center of mass 
measured in the zero 3-momentum frame, see (25J, [40] for 
details), so that they can be made to vanish by a suit- 
able choice of reference worldline. For this reason it is 
dubbed in [29 , where it was first discussed, "kinemati- 
cal hidden momentum" (by contrast with gauge invariant 
hidden momentum P^dEM' dubbed "dynamic" therein). 
In flat spacetime, if one measures the center of mass with 
respect to an inertial frame (since u a is in this case co- 
variantly constant), then Ph idI = 0. 

The inertial hidden momentum also plays a key role in 
the helical solutions allowed by the condition S a ^Up — 0, 
as discussed in detail in the companion paper |25j ; namely 
it explains how the motion of a spinning particle can be 
consistently described by an helical solution without vio- 
lating any conservation principle. And using the analogy 
in Eq. ( 29 ), we show in [2S] that these solutions are a phe- 



nomena which can be cast as analogous to the bobbing 
of a magnetic dipole in an external electric field, studied 
in Sec. III.B.l of [Ml. 
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E. Mass of the spinning particle 



We take the scalar m = —P a U a as "the proper 
mass" 4 [49] of the spinning particle. It is simply the time 
projection of P a in the particle's proper frame, i.e., the 
particle's energy as measured in the frame where it has 
zero 3-velocity. Whereas for a non-spinning particle m is 
a constant of the motion, for a spinning particle, that is 
not true in general, as we shall see. It follows from the 
definition of m that: 



dm 
~dr~ 



D pa 

dr 



U a - P a a a 



D F P a 
dr 



(31) 



i.e., dm jr is the time projection, in the CM frame, of 
the Fermi- Walker derivative of the momentum. Not- 
ing that P a a a = -Ph idEM a a , and using the orthogonality 
^hidEM^* = 0, we can rewrite this equation as: 



dm 
~dr~ 



( DP° 



V dr 



DP a 
dr 



Ua 



(32) 



Thus dm/dr equals also the time projection, in the CM 
frame, of the force DP a /dr subtracted by the derivative 
of the "electromagnetic" hidden momentum DP^^^/dr. 
Let us see the meaning of the first term. Contracting ( 17 ) 
with U a we obtain: 



DP a 
dr 



dr 



U 7 Hf3 (33) 



showing that the force has a time projection if the 
Maxwell tensor varies along the CM worldlinc. That 
arises from the force on a magnetic dipole, Eq. (|ljla) 
of Table IT] (by contrast with its gravitational counter- 
part (|l]lb), which is spatial with respect to U a ). As 
we shall see in Sec. |V A| from the point of view of the 
inertial frame momentarily comoving with the particle, 
this time-projection has the interpretation of the rate of 
work done by the electric field induced by the varying 
magnetic field. Hence, in the absence of an electromag- 
netic field, U a DP a /dr = P a a a — 0, and therefore m 
is conserved. If an electromagnetic field is present, in 
general m varies. Now, noting from Eqs. (|27[ ) and ([2 ) 
that P a a a = -kF^ciyHp, an( i putting Eqs. (31 ) and (33) 
together, we see that: 



dm 



(kF^U 7 )^ = -^ 



dr 



d,T 



(34) 



4 This is the most natural definition of the body's mass if one uses 
the Pirani spin condition, since it is the quantity which is con- 
served in the absence of electromagnetic field, as we shall see 
next. However, it should be noted that if one uses Dixon's con- 
dition S a PP/3 = instead, then it is the quantity M = \/-P a P a 
(not m), i.e., the energy of the particle as measured in the zero 
3-momentum frame, that is conserved when F a P = (see Ap- 
pendix [Cj. 



If we take \x a = aS a , being a a constant, and, since from 
Eq. p3j), B^DS^/dr = 0, we have [HI EQl EH [46]: 



dm 
~d^ 
=> m 



D 



m 



m 



(35) 

aS ■ B (36) 



where mo is a constant. Thus, if jl = o~S, m is the sum of 
a constant plus the quantity —p, ■ B which is commonly 
dubbed in elementary textbooks as the "magnetic poten- 
tial energy" of the dipole; for this reason some authors 
(e.g. [22], 06] p. 18, [36] p. 1622) have interpreted this 
term as meaning that the magnetic potential energy con- 
tributes to the mass. We argue (in agreement with the 
analysis in [64H67] ) . that what it actually means, in rigor, 
is that there is no such thing as magnetic potential en- 
ergy. —pi-B is internal energy of the test particle; in fact 
it will be shown in Sec. |VI A 2| (cf , also [55H8T] ). that, for 
a "rigid" body, it is but rotational kinetic energy, associ- 
ated with the rotation of the body around its center of 
mass. What the term —jl ■ B actually does is to ensure 
that the work done by the magnetic field on a magnetic 
dipole is zero (hence no potential energy can be assigned 
to it), as is discussed in detail in Sec, |v} 

Note that this solves an apparent paradox that has 
for long been discussed in the literature [391 EH ESI EI] : 
consider a magnetic dipole in a static non-homogeneous 
magnetic field B; a net force Fem = V(/2 • B) will be 
exerted on it, therefore it gains kinetic energy; but, on 
the other hand, we know that the magnetic field can do 
no work, because if we think about the dipole as a current 
loop (cf. Fig. [3^, below) and consider the force exerted 
in each of its individual moving charges, we see that the 
magnetic force F = q(v x B) is always orthogonal to 
the velocity v of the charges, so that no work can be 
done. Consider herein, for simplicity (the general case 
will be discussed in detail in Sec. IV B|, a setup where 



Pfi id = 0, and a globally inertial frame in flat spacetime; 
in that case, the variation of the particle's kinetic energy 
of translation is exactly canceled out by a variation of the 
particle's proper mass m, so that the particle's energy 
E = — Pn (as measured in a frame where the the field is 
time-independent 5 ) is constant. If we take, as discussed 
above, a "rigid" spinning body as our dipole model, we 
can say that the variation of translational kinetic energy, 
associated to the body's CM motion, is exactly canceled 
out by a variation of kinetic energy of rotation (about 
the body's CM) . This is in agreement with the arguments 
in |64H67j . which are re-formulated and generalized, in a 
relativistic covariant framework, in Sees. |VB| and |VI A| 



5 If in the reference frame chosen the magnetic field is time- 
dependent, the energy E = —Po of the particle, as measured 
in that frame, varies in general, but not through the work of the 
magnetic field; it is due to the work of the induced electric field, 
see Sec. [V] and e.g. [64l 168^ ■ 
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below (see in particular Sec. VI C I. In order to make these 



important points more clear, it is useful to compare with 
the cases of a monopole charged particle, and an electric 
dipole subject to an electromagnetic field, which is done 
in Appendix |B 4| There is also a force on the particle, 
which is set into motion gaining kinetic energy; but, in 
these cases, the electric field is doing work, there is a 
potential energy involved, and the gain in kinetic energy 
of the center of mass is not canceled out by a variation of 
the particle's proper mass (m is constant for a monopole 
particle, and also for an electric dipole if one assumes 
that the dipole vector is parallel transported 6 ). 

Another point we would like to emphasize is that the 
varying mass m, and the contribution — fl ■ B to it, are 
something real and physically measurable, not just a mat- 
ter of semantics or definition of the particle's mass (i.e. 
not an issue that goes away by simply calling instead m , 
in Eq. (361 above, the particle's mass): m is the inertial 
mass of the particle. In order to see that, take, again for 
simplicity, the case that P£ id = 0; we have 

— — = ma" + --U ; 
dr dr 

i.e., the projection of the force, orthogonal to U a , is ma a 
(thus, in the CM frame, DP/dr — ma). This inertial 
mass is measurable, for instance in collisions. Also (since, 
as mentioned above, in the case of a quasi-rigid spinning 
body, — pj ■ B is kinetic energy of rotation) the angular 
velocity of rotation of a spinning body is measurable. 

In the purely gravitational case, by contrast, the proper 
mass is a constant (m = mo); this means, conversely, 
that, unlike its electromagnetic counterpart, the gravito- 
magnetic field does work on mass currents. This point 
will be further explored and exemplified in Sec |VB| 



F. Center of mass motion 

Eqs. (|TJl) yield the force on the spinning particle in 
the electromagnetic and gravitational case; not the ac- 
celeration a a = DU a /dr, as P a ^ m U a . Putting 
m = mo + m' in Eq. (27), and noting that: 



e a ^yE->U° = *FP a fip+^BpU a (37) 

where we have used the decomposition ([2]), we can write 
P a as: 

P a = m Q U a - e^S^a^U 6 + (m +fi a B a )U a + *F p a ^ 



6 In a more general case, the proper mass of an electric dipole is not 
constant, but its variation substantially differs from its magnetic 
counterpart (which can be traced back to the intrinsic difference 
between the two dipole models, being a magnetic dipole a current 
loop, and the electric dipole a system of two close monopoles of 
opposite charges); it does not depend on the time derivative of 
the electric field seen by the test particle, and does not cancel 
the work done by the electric field. The case of electric dipolcs 
is discussed in detail in Appendix IB] 



This is simplified if we assume fi a = aS a \ in that case, 
as described in the previous section, m' = —[i a B a and 
thus: 

P a = m U a - e^S^a^U 5 +*F /3 V 
Therefore: 

n pa f) 

m a a = — + £ a w -C/ s (SV)+*^,[/y 



dr P10 dT 



*F 



dr 



(38) 



In a region where the charge current density j a is zero 
(most of the applications in this paper deal with vac- 
uum) we have, according to Eq. (IT 3), *F a p. T U T — 2B[ a ^. 
Therefore, if j a = 0, and using (fl7|): 



m a a = qF ap Uf) + 5 a % -M Pa Sp - *F p ' 
D 



dr 



fi-,8 



u° 



dr 



(SV) 



(39) 



Note the reversed indices in the second term, comparing 
to the expression for the force (|TJla). This leads to a 
counter-intuitive dynamical behavior, as exemplified in 

Sec. nvn 



III. WEAK FIELD REGIME: GRAVITATIONAL 
SPIN-SPIN FORCE 

If the fields do not vary along the test particle's world- 
line (so that -kFafj.^U 1 = 0), Eqs. (|TJ 3) allow for a simi- 
larity between F£ M and Fq , since B a p and W a p have the 
same symmetries. Consider the simple example of anal- 
ogous physical apparatus: a magnetic dipole in the elec- 
tromagnetic field of a spinning charge (charge Q, mag- 
netic moment fj, s ), and a gyroscope in the gravitational 
field of a spinning mass (mass M , angular momentum J), 
asymptotically described by the Kerr solution. In this 
application we shall be interested in the far field limit 
— where the non-linearities of the gravitational field are 
negligible, and one may expect a similarity between the 
gravitational and electromagnetic interactions — and in 
the slow motion regime. 

We start by briefly describing the approximations we 
will use. We take as gravitational field the usual lin- 
earized Kerr metric, Eq. (13) of 0, with the gravitational 
"potentials" $ = M/r, A = J x r/r 3 , and 6^ = 
(herein for simplicity we will be using nearly Lorentz co- 
ordinates, so that cjij in [S] becomes Sij). The electromag- 
netic potentials by their turn are exactly <fi = Q/r and 
A = p s x r/r 3 . The gravitational tidal tensors, are, con- 
sistently, linearized in the "fields"; in other words, these 
tensors consist of i) derivatives of the metric potentials, 
ii) products of the latter by the metric potentials, and iii) 
products between the derivatives; the terms of the type 
ii) and iii) are dropped, only the type i) are kept. The 
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electromagnetic tidal tensors are linear in the potentials, 
hence no weak field assumption is made in the forces ( 40 1 , 



(42 1- (43). The expression for the acceleration (461, how- 



ever, involves a term of second order in the fields, which 
is to be neglected in order to make a coherent compar- 
ison with linearized gravity. In the computation of the 
electromagnetic and gravitational tidal tensors involved 
in the forces (|42l-(|43l, Jill) -(145]), exerted on slowly mov- 



ing test particles (velocity v) , only terms up to first order 
in v are kept (as usual in slow motion approximations). 
Terms of second order in S 2 . involved in accelerations 



(46), are neglected, as usual practice in Post Newtonian 
and weak field treatments, e.g. [53 EH] (often without a 
clear justification), as they of the same order of magni- 
tude as the quadrupole contributions, as shown below. 

The general expressions for the magnetic B a p and the 
linearized gravitomagnetic H a( g tidal tensors are given 
by Eqs. (12) and (17)-(20) of E]; the expressions corre- 
sponding to the problem at hand are obtained by substi- 
tuting therein the potentials above. Let us first consider 
stationary setups; i.e. the test particle to be at rest rel- 
ative to the static observers (may be thought as the ob- 
servers at rest with respect to the central source; see point 
|8]of Sec. ID ); that amounts to set the 4- velocity u a ~ 5" 



in Eqs. Jl2), (17)- (20) of E]. For these observers, the 
gravitational tidal tensors asymptotically match the elec- 
tromagnetic ones, identifying the appropriate parame- 
ters [S7]: 



(E")i 



^3 



(E u ) 



13 j 



(H u ), 



,7-H-Ms 



(f.J) 



r {t J 3 ) Jr-SyiTj 



(all the time components vanish identically for these ob- 
servers). Therefore, for large r, the force exerted on a 
gyroscope whose center of mass is at rest relative to the 
central mass is similar to its electromagnetic counterpart, 
identifying ji s f-> J and /j, «-» S: 



G 



-B. ji Sj 



{J,S}++{» S ,»} 



-F, 



EM 



(40) 



Thus, for stationary setups, we have a gravitational 
spin-spin force similar to its electromagnetic counterpart, 
which was first found by Wald in [I]. 

Manifestation of the different symmetries — in the 
general case that the dipole/gyroscope are allowed to 
move, however, Eqs. (|lj3) make clear that the two forces 
differ significantly, due to the different symmetries of the 
tidal tensors. These differences (which lead to key differ- 
ences in the dynamics, as exemplified in next sections) 
are manifest even in the weak field and slow motion ap- 
proximation, as we shall now show. Consider now the 
test particles to be moving with 3-velocity v relative to 



the central sources 7 . The magnetic tidal tensor as seen 
by the moving dipole, B a p, can be obtained in terms of 
the tidal tensors {E u ) a p, (B u ) a p measured by the static 
observers, using the decomposition 

* Fafi;y = ^U[ a Bp] 7 — e Q/ 3 MO -Z7°' . (41) 

The force (|l]2a) exerted on it reads, to first order in v 

/| x: - * (B^r^-iE-fe^^fij, (42) 

(43) 



F EM - -B l Vi — . 



The gravitomagnetic tidal tensor as seen by the moving 
gyroscope, H a( 9, can analogously be obtained it terms of 
the tidal tensors (E u ) a p, (M u ) a p measured by the static 
observers, using the decomposition (the dual of (74) be- 
low): 



*R 



a(3 



75 



4e^ WAU ^(E«)/l-2e r a /(E«)/] 



+4(H") [ [V] Mo 



C \T 1 SfiU /TUTU 



(M u )„ T u x u v 



The force exerted on the gyroscope reads, to linear order 
in the fields, and also to first order in v: 



G 



-B. ji Sj ~ -(T)' < S j +2{W 1 ) l{ * l e>\ l v k S ] , (44) 



F° ~ —W°Si ~ - QEPy'vjSi 



(45) 



These expressions can also be obtained by setting (to- 
gether with the above substitutions for the electromag- 
netic and gravitational potentials) u a — > U a ~ (1, t7) in 
the expressions for the corresponding tidal tensors Eqs. 
(12), (17)-(20) of [B]. For this application however the 
expressions in the form given above will suit us most. 

We note that, to this accuracy, the spatial part of the 
forces, apart from global signs and a factor of two in 
the second term of (44) compared to (42), differ essen- 



tially in the fact that the former is symmetrized, whereas 
the latter is not. The dynamical consequences of this 
fact are exemplified in Sec. |IV| below. Also the time 



components are different: F^ M = — (-Fem)o = (this 
is an exact result), telling us that the total work done 
on the dipole, as measured by the static observers, is 
zero; but (i*b)o 7^ 0, telling us that there is a non- 
vanishing (and non-negligible in the linear regime) work 
being done on the gyroscope. Also one may check that 
whereas F£U a = (i.e., Fq is a spatial force), F£ M has a 
non- vanishing time projection in the CM frame which, to 
first order in v, reads: F^ M U a ~ (B u y i fijVi. The phys- 
ical significance of the time components of the forces is 
discussed in Sec. [V] 



7 Since the gravitational field we are dealing with is the linearized 
Kerr solution, which does not embody a central spinning mass, on 
more general grounds one would describe the gyroscope motion 
by saying it moves with velocity v relative to the static observers, 
i.e., observers of 4-velocity u a = u°(l,0), tangent to d/dt, in 
Boyer Lindquist coordinates. 
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It should also be noted that the forces above do 
not translate in a trivial fashion into accelerations, cf. 
Sec. |IIF[ F a in general is not even parallel to a a , as the 
test particles possess hidden momentum. Assuming the 
proportionality fjf* — aS a , from Eq. ( 39 ) we have in the 
electromagnetic case 



In the gravitational system we have, from Eq. (39): 



mo a 



F 



EM 



2B W- 



B %3 

u\lj i 



(46) 



The last two terms of ( 39 1 are herein neglected. As for the 
term ^F^Dyfi /dr, for (i a = crS a , it follows from Eq. S 
that it is of second order in the fields; since in the gravita- 
tional case we are dealing with linearized theory, we must 
coherently consider electromagnetic fields weak enough 
so that only their linear contributions are relevant to the 
dynamics. The term DPg^/dr = e a 0lS U s D{S> 3 a?) / dr is 
also negligible in this approximation if, among the many 
possible solutions allowed by the condition S a/3 Up = 
(see e.g. [2S])> we choose the "non-helical" one; ac- 
tually imposing -P^ idT ~ amounts, in this application 
(not in general! 8 ), to ensure that we pick such solution. 
This is justified as follows. 

For a non-helical representation, the acceleration 
comes in first approximation from the force F£ M 



subtracted _by DP^ idEM /dr; i.e. 



u\lj. 



kept in (46). Since |(S u ) y j u j 



from the terms we 
~ crS'/i s /r 4 , and 



D(S xua) 
~dr 



fc ;iy £ /Ltj| ~ aSQv/r i 1 we have 



III 



m 



where we used: D(S Xjja)/dT S Xu Da/dr, since 
a XjjDS/dr involves terms quadratic in the fields by 



virtue of Eqs. (23) and (46); dr/dt = v-f/r; S ~ mv Iot R, 
and we allow u rot (the rotation velocity of a point at the 
body's surface) to take its maximum value v ro t = 1; R 
is the radius of the spinning test body. The second term 
of the second expression above is neglected since it is 
quadratic in v. The order of magnitu de of the contri- 
bution DPhidi/dr to the acceleration (46), corresponds 
thus to a multiplication of the magnitude of first term 
of Eq. ( [46] ) by a small factor R/r, being r ^> R en- 
sured by the far field assumptions. It is worth noting 
also that, independently of any weak field or slow mo- 
tion consideration, if no assumptions are made about the 
multipole structure of the body, the contribution P£ id i to 
P a could not even be kept in a consistent pole-dipole ap- 
proximation, as it is of the typical order of magnitude of 
the quadrupole contributions to the hidden momentum, 
which are to be neglected; this is discussed in detail in 
Appendix [C2| 



In general one cannot a priori neglect the term S X u a based on 
the arguments above; for instance if the particle is charged and 
suffers a Lorentz force, cf. first term, of Eq. | |39| l; that yields a 
contribution to the acceleration of zeroth order in S. 



F 



G 



-(E u ) ji Sj + 2(E u ) l( - l e 3 \ l v k S J (47) 



Again the last term of (39) is negligible if we assume 



the non-helical representation, as by arguments entirely 
analogous to the ones above we show it is of the order 
~ \(M. u y i Sj\R/r. But note that in the purely gravita- 
tional case (by contrast with electromagnetism in gen- 
eral, see Footnote [8]), we can actually impose (to this 
degree of accuracy) the condition Phidi ~ as a means 
of picking the non-helical representation. Note also that 
(since there is no other type of hidden momentum) this 
effectively amounts to impose the momentum velocity re- 
lation P a ~ mU a . 

We shall now compare these results with the ones 
recently presented in 29 . We note that Eq. (46) 
agrees with Eq. (Ill) of making q — 0, neglecting 
quadrupole terms in the latter, and identifying B 1,3 = 
(B u yi and = (B u )^ , being E and B, in the nota- 
tion of |29| , the electric and magnetic fields measured by 
the static observers. It then reads in our notation: 

ma 1 = (fl^.y - e jkl v k (E u ) H ^ + ^ (£")^V (48) 

And Eq. Q agrees with Eq. (113) of 29 , provided 
that the quadrupole terms and terms quadratic in v in 
the latter are neglected, and one identifies £? ! J = W 3 /2, 
S l 'i = E 1 - 7 , where £ is the gravitoelectric field {£ = — G 
in the notation of 34 J, B is the gravitomagnetic field 
(B = — H/A in the notation of [31]), it reads in our nota- 
tion 



ma 1 = -(IHn/S-' +2e 3kl v\E G ) l ^Si 



\ {E G ) l . M v k S> 
(49) 



The pair of Eqs. (48)-(49), is equivalent to (46)-(47l, 



only written in a different form. The main difference is 
that the former makes linearized gravity look very simi- 
lar to electromagnetism; whereas the latter makes more 
transparent key differences between the two interactions. 
Namely the crucial ones which occur when the test par- 
ticle moves (even in the slow motion, weak field regime) 
and are manifest in the fact that the tensor in the second 
term of ( |47[) is symmetrized, by contrast with (46). In 
the form (|48|)-( 49 ) of [29], by their turn, for the problem 



at hand, they formally differ only by a factor of 2 second 
term of (49) (the term d l [S ■ (v x Eq)] in the notation 
of |29j). comparing with the corresponding term in (48). 
This factor alone accounts for the aforementioned differ- 
ences in this formulation, which tends to mask them. The 



examples in Sees. IV A and IV exemplify how dramatic 
they can be. 

Thus, to conclude: in the literature concerning the 
weak field gravito-electromagnetic analogy (e.g. [2]) the 
gravitational force acting on a gyroscope is commonly 
cast in the form F G = V(H ■ S)/2 (being H the 
"gravitomagnetic field", defined in those treatments as 
H l = e!-fc<7ofe,j)' similar to its electromagnetic counter- 
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part Fem = V(.B • ft) (the latter being an exact ex- 
pression, holding in the inertial frame momentarily co- 
moving with the CM, where B and [2 are measured, c.f 
Sec. II B I, seemingly implying a similarity between the 
two interactions. We emphasize that such expressions 
are not suited to describe the dynamics of a gyroscope, 
as they can hold only if the gyroscope's center of mass is 
at rest in a stationary field. This is usually overlooked 
in the literature, despite in pQ, where this analogy was 
originally presented, being mentioned that it was derived 
precisely under those conditions. Otherwise (if the fields 
are time dependent and/or the test particle moves) the 
gravitational tidal tensors M a p involved are very differ- 
ent from their electromagnetic counterparts B a j3, even in 
the slow motion regime, as one can see with full gener- 
ality substituting the expressions for tidal tensors (12) 
and (17) - (20) of 6] in Eqs |!Jl). Even if one ana- 
lyzes the system from the test particle's frame (which is 
where its electromagnetic counterpart exactly holds), it 
is clear that /•/,. = —W l Sj ^ V 
more terms), see Eq. (20) of [B]. 

In this application we have seen that the two forces in- 
deed are similar when the test particles are at rest with 
respect to the sources, which allowed us to cast gravi- 
tational the force for this case as a spin-spin force, in 
analogy with electromagnetism. However in the general 
dynamical case the test particle is allowed to move, the 
two forces differ significantly even to first order in the ve- 
locity (and in the fields), cf. Eqs. Q-p)]). That is, the 



j ^ V l (i? • S)/2 (it contains 



"spin-orbit" interaction of gravity and electromagnetism 
are very different. 



Sec. \II1\ in brief 



• In the stationary, weak field regime, when the 
particles are at rest with respect to the sources, 
the interactions are very similar, and there is 
a gravitational spin-spin force analogous to the 
electromagnetic one. 

• When the particles move, differences imposed 
by the symmetries of the tidal tensors are im- 
portant and non-negligible in any slow motion 
approximation. 



IV. DYNAMICAL MANIFESTATIONS OF THE 
SYMMETRIES OF THE MAGNETIC TIDAL 
TENSORS 

A. Radial motion in Schwarzschild spacetime 

According to Eqs. |T}l) of Table |l| it is the magnetic 
tidal tensor, as seen by the test particle of 4-velocity 
jja _ Jj0(i^tf^ thai determines the force exerted upon 
it. Consider first a magnetic dipole in the field of a static 




Figure 1: An illustration of the physical consequences of 
the different symmetries of the tidal tensors. A gyroscope 
dropped from rest in Schwarzschild spacetime will move ra- 
dially along a geodesic towards the source, with no force ex- 
erted on it. A magnetic dipole in (initially) radial motion 
in a Coulomb field, by contrast, suffers a force. Due to the 
hidden momentum, the force is approximately opposite to the 
acceleration! 



point charge Q, and with a purely radial initial veloc- 
ity. Unlike what one might naively expect, a force will 
be exerted on the dipole; the radially moving dipole in- 
deed sees a vanishing magnetic field B a , but a nonvanish- 
ing magnetic tidal tensor B a p, which can be seen from 
the expressions ((57|) below. Taking the perspective of 



the frame comoving with the particle, this is explained 
through the laws of electromagnetic induction: the mov- 
ing dipole "sees" a time-varying a electric field; by virtue 
of Eq. (Jl]3a) (which in vacuum is a covariant form for 
VxB = dE/dt), that will induce a curl in B, i.e., an anti- 
symmetric part in the magnetic tidal tensor B a p . Which 
implies that the tensor B a p is itself non- vanishing: 



Bt 



a/3] 



and therefore (except for the special case v \\ jl) a net 
force will be exerted on the dipole; computing it explicitly 
we obtain: 



f em — 0; -Fem — B 



W= -^-(vxfty . (50) 



Note this important point: for this configuration, actu- 
ally B a p = B[ a p], i.e., the force Fg M comes entirely from 
the antisymmetric part of B a p. Let us now consider 
the analogous gravitational setup, a gyroscope in radial 
motion in Schwarzschild spacetime. Since the gravit- 
omagnetic tidal tensor M. a p is symmetric (in vacuum): 
H[ Q( 8] = 0, we expect, in the spirit of the analogy, H a/ 3 as 
"seen" by the gyroscope to vanish under these conditions, 
and therefore, by virtue of Eq. (JTJ lb) of Table [IJ 

Fg = -W pa S = 

Indeed this is the case, as one may check comput- 
ing explicitly M^Sp = -kBP^UJJvSp, with U a = 
(JJ , U r , 0, 0), for this metric. This means that a gyro- 
scope in radial motion moves along a geodesic (it is the 
trivial solution of the equations of motion with Pirani 
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condition, see Appendix CI), with no forces applied on 
it. For instance, a gyroscope dropped from rest will fall 
into the singularity moving in a straight line. 

Finally, note that the force in the electromagnetic 
problem, due the hidden momentum, does not translate 
in a simple fashion into an acceleration. This is manifest 
in Eq. (39); for flat spacetime, and a particle whose only 
electromagnetic moment is fi a (q = 0), it reads: 



m a 



D 



{S p a?)U s - *F, 



dr 



The last term vanishes if one assumes fx = aS, since: 
B a (U) = for radial motion; thus, from Eq. (p3|, 



Dfi^/dr = S v a v U^ and *F p a D^/dr = -B{U) a S v a v = 
0. The second term can also be taken to a good approx- 
imation as being zero, which, as explained in Sec. |III| 
to an accuracy of order 0(S 2 ), amounts to say that we 
pick the "non-helical" solution allowed by the Mathisson- 
Pirani condition. Therefore, since, in this application, 
B( a /3) — 0, we are led to the conclusion that m a a rj 
_ p a I 

EM ' 



B. Equatorial motion in Kerr and Kerr-de-Sitter 
spacetimes 

In this application we compare equatorial motions of 
gyroscopes in the Kerr and Kerr-de-Sitter spacetimes to 
magnetic dipoles in the field of a spinning charge, to point 
out another important physical consequence of the dif- 
ferent symmetries of the tidal tensors: whereas in the 
electromagnetic case B a p can never vanish for a mov- 
ing dipole, as a consequence of the laws of electromag- 
netic induction, in the gravitational case there is a ve- 
locity field for which H Q/ 3 = 0, i.e., a gyroscopes moving 
with such velocities feel no force. This opens the possi- 
bility of circular geodesies for spinning particles, which 
we show to exist in Kerr-de-Sitter spacetimes. In or- 
der to obtain those results we make use of another ex- 
act gravito-electromagnetic analogy [T1TH3] (other than 
the tidal tensor and the inertial force analogies of Sees. 
II B and II C), this time a purely formal one (cf. [2]), 
between the decomposition of the Maxwell tensor F al3 
in electric and magnetic fields and the decomposition of 
the Riemann tensor (in vacuum) in electric and magnetic 
tidal tensors. In particular, from the formal analogy be- 
tween the scalar invariants of F a @ and the second order 
invariants of R a /3 7 s (which will be studied in more detail 
in the forthcoming paper |15j ) we note that in the equa- 
torial plane of Kerr spacetime the gravito-magnetic tidal 
tensor vanishes for some velocity field, in analogy with 
the vanishing of the magnetic field (not the magnetic 
tidal tensor!) in the corresponding electromagnetic sys- 
tem. The later means that a magnetic dipole with those 
velocities does not precess (apart from Thomas preces- 
sion); for completeness we also investigate the analogous 
gravitational problem — check if there is a velocity for 
gyroscopes such that they do not "precess" relative to 



the distant stars; the answer is affirmative, and turns out 
to asymptotically match its electromagnetic counterpart. 
This result will be of importance in the context of the 
physical significance of curvature the invariants, and the 
clarification of their relationship with gravitomagnetism, 
to be discussed in [151 . 



1. A magnetic dipole in the field of a spinning charge. 

Velocity field for which B a — 0. — We start by the 
electromagnetic system which will serve as a guide for 
the gravitational case. With respect to a unit timelike 
4- vector u a , the Maxwell tensor (6 independent com- 
ponents) splits irreducibly into the two spatial vectors 
(3 independent components each) (Eff) a = F a ^up and 
(B u ) a = *F al3 u/3, cf. decomposition In), which are a co- 
variant definition for, respectively, the electric and mag- 
netic fields as measured by an observer of 4- velocity u a . 
Despite both (E u ) a and {B u ) a depending on the observer 
4- velocity u a , combining them one can construct the two 
second order scalar invariants (e.g. [TT1 H"2l 170] ): 



1 



F a pF° 



(51) 



E a B a = --F a p*F a ? (52) 



These are actually the only two independent scalar invari- 
ants one can construct from the Maxwell tensor. They 
have the following physical interpretation [T51 [701 ITT] : 
i) if E a B a ^ then the electric E a and magnetic 
B a fields are both non-vanishing for all observers; ii) if 
E a E a - B a B a > (< 0) and E a B a = 0, there are ob- 
servers u a for which the magnetic field (B u ) a (the electric 
field {E u ) a ) is zero. 

The electromagnetic field produced by a spinning 
charge (magnetic moment fL s ) is described by the 4- 
potential A a = (<p,A): 



Q 1 ft, s x r _ 1 /i,, sin I 



(53) 



The invariant structure for this electromagnetic field is: 



E 2 - B 2 = 



^(5 + 3 cos 20) 
2r 6 



> 



p> 3 2u s Qcos9 

E ■ B = = (= m the equatorial plane) 

(54) 

which tells us that in the equatorial plane 8 = it/ 2 there 
are observers who measure B a to be locally zero (since 
E ■ B = and E 2 — B 2 > therein) . It is straightforward 
to obtain the 4-velocity of such observers. The magnetic 
field (B u ) a = -kF a ^up seen by an arbitrary observer of 
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4-velocity u a — (it*, u r , u , u^) is given by: 



2// s cos0 t 

D = 5 U 



B" 



Bf = 



sin 6 



B*=^ (2u r cos6 + ru e sint 



Thus, the condition B r = implies 8 = tt/2 (i.e., equato- 
rial plane, as expected); in the equatorial plane, B l = 
implies u e = 0, and B e = implies u^/u 1 = /i s /(Qr 2 ). 
If we additionally assume that the charge and mass are 
identically distributed in the body, its gyromagnetic ra- 
tio is jU s /J = Q/2M, and we therefore conclude that 
observers with angular velocity 



J 



2Mr 2 



v 



(B=0) 



(55) 



see a vanishing magnetic field in the equatorial plane. 
No restriction is imposed on the radial component of the 
velocity, apart from the normalization condition u a u a — 
— 1; that is, the velocity for which the magnetic field van- 
ishes at a given point is not unique — it comprehends a 
class of observers [57]. The velocity field corresponding 
to the case u r — is plotted in Fig. [2^). The vanishing 
of B a for these observers comes from an exact cancella- 
tion between the magnetic field generated by the relative 
translation of the source and the one produced by its 
rotation. And it means that a magnetic dipole possess- 
ing a velocity of the form ( 55 I does not undergo Larmor 



precession, since the second term of Eq. ( 23 ) vanishes 



Hence dipoles moving in circular lines tangent to the 
velocity field of Eq. (55) undergo Fermi- Walker trans- 



port. Note that this does not mean they do not "precess" 
relative to the distant stars (in this case, a global inertial 
frame); they do, since such worldlines are accelerated, 
due to Thomas precession. The latter, from the point of 
view of the inertial frame, is encoded in the first term of of 



Eq. ( 23 1 ; and from the point of view of the frame comov- 



ing with with the dipole, where Eq. DS l /dr = holds, 
and with axes fixed to the distant stars, it is encoded in 

ft*, 



the the connection r|, 

0] 



cf. Eq. (24). Let us 



denote hereafter by ft* the angular velocity of rotation, 
relative to Fermi-Walker transport along the particle's 
worldline, of a comoving spatial triad fixed to the dis- 
tant stars; and by — ft* to t the total precession rate of the 



particle's spin relative to this frame. Thus, cf. Eq. (24): 



dS 
It 



ft* 



x S: ft* 



-ft* - <jB 



(56) 



In flat spacetime, ft* equals minus the Thomas preces- 
sion; it can be shown (e.g. [53] p. 564) that ft* = 
7« x a/(7 2 + 1), where v is the velocity of the dipole 
relative to the distant stars, and a the spatial part of 
the acceleration a a = DU a /dr. Hence, ft* = only 
for a magnetic dipole in straightline motion, momentar- 
ily moving with a velocity obeying (55). This observation 



is important for the comparison with the gravitational 



analogue in Sec. IV B 2 



B a p never vanishes. — The force (JTJla) exerted on the 
dipole, however, does not vanish, as it is only the field, 
not its derivatives (or its tidal tensor), that vanishes. Ac- 
tually, the magnetic tidal tensor B a p, as measured by 
a moving dipole, never vanishes. As seen by a generic 
observer, it has following components in the equatorial 
plane: 

B r g = a(r 2 Qu -3^ s u*); B 6r = a(2r 2 Qu <t ' - 3/x s u'); 



B, 



r(j> 



-aQr 2 u u \ B, 



r = —2aQr 2 u e 
Be<j, = aQr 2 u r ; B^g = -aQr 2 u r 
B tr = ia^ s u e ] B t e = 3afi s u r 



(57) 



with a = l/?' 3 . Thus we see that in order to make B^ a ^ 
vanish, we must have u e = u r = and 



Mr 2 



J {B {af>) =0) 



(58) 



(differs from a factor of 2 from the velocity (55) which 
makes B a vanish; the second equality again assumes 
that the charge and mass are identically distributed in 
the source). However, B[ a p] only vanishes if v = 0; 
hence it is not possible to find an observer for which 



B a /3 — B 



(a/3) 



B[a/3] — 0. 



Again, the fact that B[ a p] 
cannot vanish for a moving observer is a direct conse- 
quence of Maxwell's equations/the laws of electromag- 
netic induction: a dipole moving relative to the spin- 
ning charge always sees a varying electromagnetic field; 
that endows B a p with an antisymmetric part, by virtue 
(from the point of view of the MCRF) of vacuum equa- 
tion VxJJ = dE/dt, or covariantly, by Eq. |l]3a). Note 
that this is true even if one considers a dipole in cir- 
cular equatorial trajectories around the central source: 
D * F a p/dT — 2_B[ Q(3 ] ^ along such worldline, which is 
due to the variation of the electric field along the curve 
(it is constant in magnitude, but varying in direction). 



2. A gyroscope in Kerr spacetime 

Velocity field for which the gyroscope does not "pre- 
cess". — Drawing a parallelism with the previous section, 
we start by posing the question: is there a gravitational 
analogue to the velocity field ( [55| (i.e., a cancellation be- 
tween translational and rotational gravitomagnetism)? 

Before answering this question one must first note some 
subtleties. What we did in the previous section was find- 
ing a velocity field, given by (55), for which B a = 0; so 



to study the analogous gravitational problem, the first 
difference to note is that herein the problem does not 
translate into the vanishing of physical field such as the 
magnetic field B a . As explained in detail in the compan- 
ion paper {34j, the gravitomagnetic field H, which comes 
at first glance as the most natural candidate if one is to 
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discussed in Sec. II C The only part involved in it is tt; 



Figure 2: a) Equatorial plane of a spinning charge. Black 
arrows: velocity field £?(b=q) , which makes the magnetic field 
B a vanish; magnetic dipoles in straightline motion, momen- 
tarily with such velocities, do not precess relative to the dis- 
tant stars, b) Equatorial plane of Kerr spacetime. Black ar- 
rows: velocity field i?(n*=o) f° r which a gyroscope in straight- 
line motion momentarily does not "precess" (with respect to 
the distant stars); asymptotically it matches its electromag- 
netic counterpart. Red arrows: velocity field Vm=o) which 
makes the gravito-magnetic tidal tensor W a p vanish; means 
that gyroscopes moving along trajectories tangent to V(a=o) 
feel no force, Fq =0. If A > (Kerr-dS spacetime) circular 
geodesies for gyroscopes do even exist (Sec. IV B 3). V(w=o) 
has no electromagnetic analogue: for a moving dipole always 
B a p 0, by virtue of Eq. B^p] = -kDF a p/dr, generically 
implying Ff M / . 



draw a parallelism with electromagnetism, is not a phys- 
ical field, it is a mere artifact of the reference frame; so 
obviously it is not reflected in any field invariants, and 
therefore an analysis in the spirit of the one in the previ- 
ous section is immediately ruled out (an analogous invari- 
ant analysis can be used instead to solve for the vanishing 
of the gravitomagnetic tidal tensor H Q| a, as we shall see 
below). Moreover, unlike B a , which depends only on a 
given observer's 4- velocity (not on its derivatives), H is 
defined with respect to a congruence of observers, and 
consists on the sum of the vorticity uj of the observer 
congruence plus the rotation (l of the local tetrad rela- 
tive to Fermi- Walker transport 34.. Thus in the most 
general case H can be anything, depending on the refer- 
ence frame one wishes to set up. Hence, more than some 
naive notion of gravitomagnetic field, what makes sense 
here to compare the are the physical effects. We know 
that a gyroscope at finite r in Kerr spacetime "precesses" 
with respect to a frame anchored to the distant stars 



(how such a frame is set up is explained in Sees. II C and 
in [101 This reflects an effect — frame dragging - 
which is physical (even though it has no local meaning, 
and can only be measured by locking to the distant stars 
by means of a telescope), and analogous to the preces- 
sion of a magnetic dipole subject to a magnetic field, as 



which (in this case that the reference tetrad is locked to 
the distant stars) gains thereby a physical meaning, that 
we can compare with electromagnetism. 

The total precession of the test particle relative to the 
distant stars ( 56 ) is thus a quantity one can define both 



for the magnetic dipole of the previous section and for 
the gyroscope in Kerr spacetime. Hence it will serve as 
a guide to find the gravitational analogue of the veloc- 
ity field (55). Firstly we note that, as discussed in the 



previous section, the velocity field ( 55 ) means not that 
f2*tot vanishes for a dipole moving in circular lines tan- 
gent to U(b=o) 7^ in that case due to the Thomas 

precession); but that Q*tot — f° r a dipole carried by 
an observer moving in a straight line with a velocity that 
momentarily coincides with V(b=o)- So the analogous 
problem in gravity seems to be finding the velocities at 
which observers in straightline motion carrying spatial 
triads with axes fixed to the distant stars must travel, 
in order for such triads to be Fermi- Walker transported 
(i.e., il* = f2*tot = 0, so that a gyroscope it carries 
does not "precess" relative to the distant stars). How- 
ever a difference must be noted: in the electromagnetic 
case the vanishing of f2*tot comes from the fact that by 
moving in straight-line, the Thomas precession vanishes, 
and by having a velocity momentary obeying (551 the 
Larmor precession —o~B also vanishes; the two appear- 
ing in Eq. (24) as clearly separate terms. In the grav- 
itational case, by contrast, we only have tt± — — f2*totj 
which is itself the total rate of rotation of the reference 
tetrad relative to Fermi- Walker transport, and embod- 
ies the Thomas precession. The latter is not zero for 
straightline motion in curved spacetime 9 , since these are 
accelerated observers in general. Hence one might won- 
der, which one is the problem analogous to the vanishing 



of B a for the velocity field (55) of the previous section: 
the velocity field such that SI* = for a gyroscope mov- 
ing with zero acceleration, or moving in "straight line"? 
(The result will depend on this choice, due to the Thomas 
precession, that contributes in the latter case). We ar- 
gue it is the latter, since, as will be shown below, it is 
the one that asymptotically matches the electromagnetic 
counterpart, and coincides with vanishing of the Post- 
Newtonian (PN) gravitomagnetic field Hp^, as defined 
in the comoving PN frame (note that -ffpN is meaning- 
ful 10 in this context, since the PN frames are anchored to 



In fact it is not clear how one should define "straight line" in 
a curved spacetime, since light rays and spacelike geodesies are 
similarly deflected by the gravitational field. In the case of the 
equatorial plane for the Kerr solution, however, axisymmetry 
provides a good definition of "radial direction" and "tangential 
direction", which we can explore to define straight line, see the 
discussion following Eqs. ( |li4| . 

One might wonder whether the exact problem at hand could be 
solved using the gravitomagnetic field H, i.e., instead of basing 
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the distant stars; note also that the Thomas precession is 
cast as a part of Hpn |Z2])- It can thus be interpreted as 
a cancellation between translational and rotational grav- 
itomagnetism, just like the vanishing of B a for observers 
moving with velocities V(b=o) arises from an exact can- 
cellation between the magnetic fields generated by the 
rotation and relative translational motion of the source. 

We shall now compute such velocity field. The Boyer- 
Lindquist line element for 8 — tt/2 can be written as 



To compute the Fermi- Walker transport law we need 
the connection coefficients TJ?-, defined, as usual, by 

VEp.Ec, = T? Ea- They can be codified in the Cartan 
connection forms 



, ,a pa , .7 

Id 3 — 1 - aid , 

7/3 ' 



(60) 



which form the unique solution to the Cartan structure 
equations 



-V dt 



ds 2 



where V = 1 



2Ma 
rV 



dci> 



W^dr 2 +r 2 WV~ 



2M 



and W 



1 



2M 



This 



form of the metric suggests the orthonormal coframe 



= V?dt 



2Ma 

TV 2 



= W~'dr; 



= rW2V-*d(t>, (59) 



whose dual frame {-E5, Ef, E^} can be defined geomet- 
rically. Indeed, E§ is simply the unit vector tangent to 
d/dt, which is the the unique (up to a constant) Killing 
vector field which is timelike at infinity. Thus the con- 
gruence with four-velocity Eq corresponds to the only 
observers who can claim to be "at rest" with respect to 
the distant stars — the static observers, cf. point [8] of 
Sec. ID The vectors Ef and Ei can be geometrically 
defined by the condition of being respectively orthogonal 
and tangent to the surfaces of constant r. These sur- 
faces have geometric meaning, as they are the orbits of 
the equatorial plane's isometry group (generated by the 
Killing vector fields d/dt and d/d<fi). 



duj & + Ld & p A uo = 0. 



A straightforward calculation yields 

J, = u * A = ¥ r v- 1 wiJ-^v- 1 J 



v 

Ma, 



(61) 



(62) 



Id 



w 6 



1 - 



Ma, 



-V~ w + 
AM + AM 2 Ma 2 



(63) 



We now have to decide how to define the concept of 
"pointing in the same direction" for spacelike vectors or- 
thogonal to Eq. Clearly Ef and Ei do not point in 
the same direction throughout the equatorial plane, be- 
cause they fail to do so even in the flat Minkowski case 
M = a = 0. By analogy with what one would do in this 
case, we consider the frame 



El 



En 



(64) 



cos 4>Ef — sin 4>E^ 
E^ = sin <j)Ef + cos <f>Ei 



the analysis in the local tetrad associated to the single observer 
carrying the gyroscope, one could embed it (in some neighbor- 
hood) in a congruence, and then take the frame adapted to it 
(uj = f2), in order to define a gravitomagnetic field H = 2lj. 
That would seemingly allow us to take advantage of the analogy 
between the transformation laws for the gravitomagnetic and 
magnetic fields (see Eqs. (4.14) and (8.3) of |9|), in the spirit 
of what is done in 1721 ; and to solve this problem by simply de- 
manding H(U ) = 0, just like the in the electromagnetic case the 
solution (I55J1 amounts to solve for v the Eq. B(U) = 0. This 
scheme however is of no use in our exact problem, since the H 
computed by this procedure does not tell us about the preces- 
sion relative to the distant stars; what it measures is the rotation 
of a tetrad locked to the vorticity of the congruence relative to 
local Fermi- Walker transport. This is not in general a compar- 
ison of the latter with the "distant star" fixed tetrad (only in 
very special cases, see Sec. pi C\ and 34 ]), and in particular it 
is never so if the congruence has shear. One may check that the 
vorticity (thus the H of the adapted frame) of the congruence of 
"boosted" observers |66| does not match the rotation fl* (relative 
to Fermi- Walker transport), of the boosted tetrads fixed to 
the distant stars: uj = H /2 ^ H*. The scheme would only work 
in the weak field (Post Newtonian) regime, and for PN boosted 
frames (frames with constant coordinate velocity), since to PN 
accuracy these frames are fixed to the distant stars, and their 
shear is negligible (only the expansion is not), see 134) . 



This frame, which we will call the fiducial frame, will 
define what we mean by "pointing in the same direction 
with respect to the distant stars". We will say that a 
gyroscope does not precess with respect to the distant 
stars if the components of its spin vector in the fiducial 
frame remain constant 11 . Moreover, this frame will also 
define what is meant by "straight line motion". We will 
say that a particle is moving in straight line with constant 
velocity if the components of its velocity vector in the 
fiducial frame remain constant. From the Lie brackets 



M 



jV-'W-zEq; 



[Eo,E^}=0, 



we see that [E^,E^] and [E^,E^\ are proportional to -Eg, 
meaning that the distribution generated by E^ and any 
linear combination of E^ and Er with constant coeffi- 
cients is integrable. In other words, particles moving 



11 Setting up this frame experimentally is another matter; it does 
not suffice to point Ej, say, towards the apparent position of 
some distant star, because of gravitational lcnsing. 



21 



in straight line in the same direction but with different 
speeds meet the same static observers throughout their 
paths (albeit at different times). This shows that our 
definition of straight line motion is consistent. 

To compute the connection coefficients of the fiducial 
frame we notice that its dual coframe is obtained from 
Eqs. (64) above by replacing E& uj a , and hence for 
4> = (which can always be assumed due to the axial 
symmetry) we have dui° = du>°; 



duj = dio r 



dijj 1 = dujt 



>Aw r 



From Cartan's structure Eqs. (61) we obtain 
■ 5 = w f, UJ ^ = w 5 = 

1 



or 



ri ~ 2 

U f —d<f> = uA 



(65) 



Consider now the frame obtained from the fiducial frame 
( 64 ) by a global boost of velocity v in the direction, 



E 



Ei = JE V , 

vjE 6 + 



(66) 



E 2 



(where 7 = — v 2 ). Its dual coframe is w 1 = U) ; 



vjlu ; 



-vjoj + 7W 2 . 



It is well known that the matrices and Q of the connec- 
tion forms in the fiducial and boosted frames are related 
by 



where 



S 




(67) 



(68) 



is the change of basis matrix (hence dS = in this case). 
Therefore we have 



vfUJ ^ + 7w 1 2 . 



(69) 



The components of the spin vector of a gyroscope at rest 
in the boosted frame will remain constant if and only if 
Ei and E§ are Fermi- Walker transported along Eq, that 
is, if and only if the space components of V^.-Ej and 
Vjj. i?2 vanish. This is equivalent to having Tq 2 = 0, 



that is, to u> \ having no component in ur. Using (62), 



( 65 ) and ( 69 ) we can write this condition as 
2 r/ 5M 



v 



r 

Ma 
+ 1 = 



r 



6M 2 

r2 



2Ma 2 

r3 



- T/3 



(70) 



(assuming a 7^ 0; if a = 0, i.e. in the Schwarzschild case, 
this condition is equivalent to v = 0). This equation is 
quadratic in v, and in principle has two solutions; how- 
ever, since the product of these solutions is 1, only one 
of them will correspond to the velocity of the for which 
f2* = 0. In the far field region, where M/r -C 1, (70) 
approaches 



9 2r 

v 2 v + 1 =0, 

a 

whose physically acceptable solution approaches 

r— > 00 d 



2r' 



(71) 



(72) 



Note the striking similarity with the velocity field (55) 



that makes B a vanish in the analogous electromagnetic 
problem: both velocities depend only on r and on the ra- 
tio a = J/M (if the charge and mass are identically dis- 
tributed, in the electromagnetic case), and asymptotically 
they match. This weak field limit is the one one would 
obtain using a boosted Post-Newtonian metric (e.g. [73]), 
and solving for -ffp N = 0, using the (most common) naive 
notion of gravito-magnetic field defined as the curl of the 
off-diagonal part of the metric: iJpN — -Vxi, Ai = goi, 
since for an orthonormal tetrad (attached to the back- 
ground inertial frame, which in this case means the dis- 
tant starts), the expression dHp N /dr = (S x 
always holds, as shown in [57] , 

Close to the ergosphere r = 2M, on the other hand, 



(70) approaches 



2v + 1 = W(n t = 



=0) 



(73) 



and hence the velocity for which there is no precession 
approaches the speed of light. These results will be of im- 
portance for the companion paper |15j (namely to clarify 
misconceptions on the relationship between gravitomag- 
netism and the curvature invariants discussed below). 

As discussed above, this is a cancellation of trans- 
lational (which comes as combination of Thomas and 
geodetic precessions) and rotational gravitomagnetism 
(which contributes as the Lense-Thirring precession). 
Another problem which in this spirit would also arise 
as a natural candidate to be the analogue of the elec- 
tromagnetic problem of the previous section would be 
considering gyroscopes in geodesic motion, and look 
for the velocity field for which (momentarily) they 
do not precess; this eliminates Thomas precession, 
and we would be effectively looking for a cancellation 
between the geodetic and Lcnse Thirring precessions. 
This is an equally interesting problem, but not the 
one that asymptotically matches the electromagnetic 
counterpart. Let us denote such velocity field by v'; 
in the PN regime it can be computed from e.g. Eqs. 
(3.4.38) of [3], or (40.33b) of [3]], setting therein a = 
and the PPN parameters 7 = fi = 1 (note also that — fi* 
in our notation corresponds to fl in |3 , and Q in |31|). 
Solving for v — v^e^ the equation i\ = 0, we obtain an 
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asymptotic velocity v' ~ 2a/3r — (4/3)v(b- 
differs by a factor from ( 55 ) and ([72| . 



o) 



which 



Velocity for which M a p = 0. — We have seen in 
Sec. |IVB1| that, by virtue of Eq. (|T}3a) of Table [T] (i.e, 
the laws of electromagnetic induction), for a moving ob- 
server, Bt a m ^ 0; and thus B a p ^ 0, always. But we 
have seen also that, in the equatorial plane, there is a 
velocity field for which -B( Q; a) = 0. Since the gravitomag- 
netic tidal tensor is symmetric in vacuum: = H( Q( 3) , 
we expect, in the spirit of the analogy, that H Q ^ = in 
the analogous gravitational setup. We shall now find the 
velocity fields for which H Q( g = 0. It turns out that this 
can be done by a procedure very similar to the problem 
of finding the velocity field W(b=o) I0r which the mag- 
netic field B a (not the tidal tensor B a ^ vanishes, due 
to the formal analogy between the decompositions of the 
Weyl and Maxwell tensors in electric and magnetic parts 
pn~HT5] . and the quadratic invariants they form. 

In vacuum the Riemann tensor becomes the Weyl (10 
independent components), which can be irreducibly de- 
composed (see e.g. |13|). with respect to a unit timelike 
4- vector u a , in terms of two spatial tensors, the gravito- 
electric E Q/ 3 and gravito-magnetic H Q/ g tidal tensors: 

R$ = 4{2« [a ^+ 5[ ^}(E«)^ (74) 
+2 {ea^CffytW + e^ v {m u )^u a] u v ) 

In vacuum, E a/ g and M. a p are both symmetric and trace- 
less, possessing 5 independent components each, thus en- 
coding the 10 independent components of R a ^s- 

The decomposition ( 74 1 is formally analogous to the 
decomposition (JTJ of F^in terms of E a and B a . More- 
over, the gravitational tidal tensors form the second order 
scalar invariants formally (see discussion in Sec. 4 of | 14j) 
analogous to the electromagnetic invariants ( 51 )-( 52 1 : 



E Q7 E« 7 



^ 7 H Q7 = \R aM sR aMS = ^R R (75) 



E Q7 I 



' ' ^ a ^ s = — *R R 



— Y^Raf3~tS * R 



16 



(76) 



which are usually combined into I = (R-R+i*R-R)/8. 
However, these invariants are not sufficient 12 for an anal- 
ysis like the one we did for the electromagnetic case. 



Whereas the invariants ( 51 )-( 52 ) are the only two inde- 
pendent scalar invariants of F a ^ , in the case of R a /3-yS 
(for vacuum) there also two independent cubic invari- 
ants, given by A = R^R^ pa R pa af} /l6 and B = 
R a P^ R ^ pa * R pa alj /16, and combined into J = A — iB. 
It turns out (cf. [TH [75l [76]) that one obtains formally 
equivalents statements to (i)-(h) of the electromagnetic 



12 We thank L. Wylleman for his input on this issue. 



case, replacing F by R, provided that the condition M = 
I 3 / J 2 — 6 > (real or infinite) is added to ii); that is: i) 
*R • R =/= =>■ E Q7 and H Q7 are both non- vanishing for 
all observers; ii) *R • R = 0, R • R > 0, with M > => 
there are observers for which H Q7 vanishes ("Purely Elec- 
tric" spacetime) 13 . Further details and comments on this 
classification shall be given in [T5] . 

For the case of Kerr spacetime, which is of Petrov type 
D, the condition M > (real) is satisfied, since I 3 = 6 J 2 
(see e.g. [76]). Thus one only has to worry about the 
invariants ((75j)-([76|), which have the structure (the exact 



expressions are presented in |15j): 



E Eq, 7 - H 7 H, 



6M 2 



> o 



E Q7 H a7 w — — — ^ (3S _^ ( = o in the equatorial plane) 
formally analogous to its electromagnetic counterpart 



(54). Note in particular that the result E a7 H Q7 = 
for the equatorial plane (8 = it/2) is exact. This means 
that in the equatorial plane , there are observers for which 
the magnetic tidal tensor vanishes. It is straightforward 
to determine the 4-velocity of such observers. In the 
equatorial plane, the non-zero components of the gravito- 

-kRaiipvU^u" seen by an ar- 
are 



magnetic tidal tensor H f 
bitrary observer of 4-velocity u a — (it* 
given exactly by (a = 3M/r 3 ): 



[2a z 

a (a 2 1 
aa (au 



r 2 ) u^v} 



2 ) (< 



- a (a 2 



u 1 ) u e 



(77) 



te0 = eta [(a 2 + r 2 ) — au~\ W 
l et = -a[(a 2 +r 2 )ii -au*]u r 



-2aau r u 



r„,e 



Hi ! 



a (2a 2 + r 2 ) u r u e 



Then it is straightforward to show that in order to 
make all components M. a p vanish, we must have u 6 = 
(i.e., the observer must move in the equatorial plane, as 
expected and in analogy with the electromagnetic exam- 
ple above) and 



J (H=0) 



(78) 



=o) 



Thus, observers with angular velocity = 
vanishing gravito-magnetic tidal tensor in the equatorial 
plane. Again, no restriction in imposed on u r , apart from 
the normalization condition u a u a = — 1. The velocity 



13 The case *R. ■ R. = 0, *R ■ R. < would mean that there would 
be observers for which E a7 vanishes (but no "Purely Magnetic" 
vacuum solutions are known, and it has been conjectured that 
do not exist , see e.g. [771 178] ). 
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field corresponding to the case u r = is plotted in Fig. 
[2}d. It is interesting to note that, asymptotically, this 
field matches the velocity (581 for which the symmetric 
part of the magnetic tidal tensor B a p vanishes in the 
electromagnetic analogue (and, up to a factor of 2, the 
velocity (72) for which the precession vanishes). 

As discussed above, the velocity field ( 78 ) has no elec- 
tromagnetic analogue. The magnetic tidal tensor B a p 
can never vanish for a moving observer, due to Maxwell's 
(vacuum) equation VxS = dE/dt; or, in other words, 
due to the laws of electromagnetic induction. The fact 
that in the gravitational case one can find a velocity field 
for which H ai g = is an illustration of the physical con- 
sequences of the different symmetries of H Q( g compared 
to B a p, and means that there is no gravitational ana- 
logue to the curl of B induced by a time- varying field. In 
other words, it signals the absence of electromagnetic-like 
induction effects in gravity. 

Note that these fundamental differences between the 
two interactions are manifest even in the weak field and 
slow motion regime, for taking the field to be weak, either 
by going far away from the source, or by taking a to 
be small, that only amounts to make the velocity for 
which Fq vanishes (by contrast with F£ M ) smaller, since 
\v\ s» a/r. That illustrates how misleading can be the 
usual treatments in the literature (e.g. [2]) concerning 
the topic of "gravitoelectromagnetism" in the framework 
of the linearized theory, which naively cast the force on 
a gyroscope in an expression of the type Fq — — V(5 • 
Bq), identical to the electromagnetic force on a magnetic 
dipole. As discussed in Sec. Ill the similarity between 
the two forces occurs only in the very special case that 
the gyroscope's center of mass is at rest relative to the 
source. 

No circular geodesies for spinning material particles in 
Kerr spacetime — the vanishing of Fq for gyroscopes 



carried by a congruence of observers of the form (78) 
makes us wonder if a spinning particle can move in cir- 
cular geodesies around a Kerr Black Hole, which we will 
now check. Since the observers ( 78 1 have prograde motion 



with angular velocity , we are only interested in check- 
ing prograde circular geodesies, whose angular velocity 
reads (e.g. [79]): 



1 



(79) 



Equating this expression to ( 78 ) , we obtain 
a \[M 



a 2 + r A 



r 3 / 2 + aVM 



& r 



a 
M 



Thus, for given M and J, at r = a? jM there would be 
circular orbits along which H Q( g = 0. However such orbits 
are not possible since r does not lie outside the horizon. 

2 , thus the condition r > r + implies 



r+ = M 
M 



>M + My/1- a 2 /M 2 & 1 - A 2 < 



Vi - a 2 



where we defined the dimensionless parameter A = a/M. 
A = 1 is the extreme Kerr case; A > 1 is a naked sin- 
gularity; thus, excluding the naked singularity scenario, 
the circular orbit would exist only in the extreme case, it 
would be precisely at the horizon, and thus (apart from 
being obviously unstable), it would be a null geodesic. 
Finally, it is interesting to note that the velocity field 



( 78 1 appeared independently, in a completely different 
context (not related to tidal effects), in a recent paper, 
see Eq. (30) of [80 . Therein it is shown that the Kerr- 
Newman metric can be obtained by a simple rescaling 
of an orthonormal tetrad field in Minkowski space, con- 
structed from spheroidal coordinates in differential rota- 
tion, each spheroidal shell r = constant rotating rigidly. 
It turns out that the angular velocity of the shells is 
precisely = v^ m=0 y both in Minkowski space and in 
the Kerr-Newman metric (with respect to the static ob- 
servers); that is, in the light of our discussion above, the 
gyroscopes comoving with the shells measure vanishing 
magnetic curvature M. a p and feel no force: Fq = 0. 



3. Circular geodesies in Kerr de Sitter Spacetimes 

The failure to obtain circular geodesies for material 
spinning particles in the previous section was due to the 
fact that the angular velocity of circular geodesies in the 
Kerr spacetime dies off as r 2 / 3 , whereas the angular ve- 
locity for which H Q( 3 = dies off as r~ 2 ; in other words, 
"geodesies are too fast". But they should be possible 
in other spacetimes; in this spirit, Kerr de Sitter come 
as natural candidate, since a repulsive A should "slow 
down" the circular geodesies. Next we will show that 
this is indeed the case. In Boyer-Lindquist coordinates, 
the Kerr-de Sitter metric takes the form 



A, 



,2 E 



ds 2 = ^- (dt - asm 2 6d6) + —dr 2 + —d9 2 

A fl sin 2 i 



— - [adt - (a 2 + r 2 < ' 



X 2 S 



(80) 



where 



A r EE 


r 2 


- 2Mr + 






A 2 


A e = 


1-1 


- —a cos 




1-f 


O 

A 2 


X 


3° 


v = 


r 2 


-1- a 2 cos 2 



2 A 2/ 2 , 2\ 



(81) 



First one notices that since A ^ =>• ^ the classifi- 
cation based on scalar invariants used in the previous sec- 
tion does not apply to herein the Riemann tensor. Such 
classification still holds however generically for the Weyl 
tensor C afi p u , see e.g. [7H]- And since A does not generate 
mass currents, we have H Q ^ = M.r a p) = rl a /3, denoting 
rl a p ee *C a ^f} u U^U u the magnetic part of the Weyl ten- 
sor. Thus, solving for W a p = amounts to check the 
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conditions under which H a p vanishes, which reduces to 
the same procedure of the previous section, but this time 
using the invariants of the Weyl tensor. The invariants 
have a similar structure, similarly leading to the conclu- 
sion that in the equatorial plane, there are observers for 
which H a f3 = H Q( g = 0. Actually, the magnetic tidal 
tensor H a/ g for the metric ( 80 1 is generically obtained by 
simply multiplying expressions (77) by 9/(3 + a 2 A) 3 : 



(H 



Kerr-dSja^ 



(3 + a 2 A) 3 v 

Thus, the velocity field for which H a/ 3 



same Eq. ( 78 ) 



UKerrJ a/3 ■ 

= is given by the 



Now we need to check if this velocity field can corre- 
spond to a family of circular geodesies. We can very eas- 
ily derive the geodesic equations from the Euler-Lagrange 
equations 



dX (dii* 



dC 







with Lagrangian 



(82) 



(83) 



To compute the circular geodesies we only need the in- 
equation, 



dX 



(fjr 



1 



which for circular trajectories (9 = 0, r 
to 



(84) 

= f ) reduces 



ff00,r(Vgeo) 2 + 2g t!j>: rVf eo + g t t,r = 0, 



(85) 



where v. 

Solving this equation in the equatorial plane 6 = ^ gives 



; is the angular velocity of the geodesies. 



-Ma + far 3 ± ^/Mr 3 



gco 



l± = 



At-6 
3 ' 



- 3 - a 2 M + 4a 2 r 3 



(86) 



which reduces to standard Kerr case, Eq. ( 79 ), when A = 
0. 

There are two things we need to check, though. The 
first is that the geodesies lie outside the black hole event 
horizon (and inside the cosmological horizon), the sec- 
ond is that the geodesies are time-like. The horizons are 
located at the real roots of A r = 0, which gives us the 
equation 



2Mr 



A 



r 2 (r 2 + a 2 ) = 0. 



(87) 



To guarantee that the geodesic is time-like we obviously 
have to make sure that ds 2 < 0, which for our case gives 
the following condition 

-t 2 + 2Mr(l - v% Q a) 2 + r 2 (v% Q ) 2 (a 2 + r 2 ) 
+ ^r 2 (r 2 +a 2 )(l-< o a) 2 <0. 



To find spinning particles that follow circular 
geodesies, we have only to equate the prograde solutions 
of Eq. pBb to d78l|, 



-Ma 



far 3 



Mr 3 



A r 6 
3 ' 



r 3 -a 2 M+±a 2 r 3 



(89) 



Unfortunately, we can no longer analytically solve this 
equation for r in general, but for our present purposes it 
suffices to numerically show that such an r exists for some 
particular cases of a and A. Consider for example the case 
a/M = 0.8, AM 2 = 0.001. Solving equation Q for r, 
we find, as the only acceptable solution r ~ 14.2025Af 
(the other roots are either complex or fall within the hori- 
zon ). We can easily see from equations ([87| and ([88]) 



that this geodesic is both time-like and outside the event 
horizon. Obviously, several other solutions of (89) for 
different values of a and A are possible. We generically 
find that, for fixed a/M, decreasing values of AM 2 cor- 



respond to solutions of (89 1 with increasing values of r. 



These solutions seem to be all unstable, however 



Sec. \IV\ in brief — the physical consequences of 
the different symmetries of B a p and M. a f3 

• In clectromagnetism, due to vacuum Eq. 

B [af}] = *F al3 . 7 liy2, a force F% M = B a ? ^ 
is exerted on the dipole whenever it moves in 
an inhomogeneous field (except for very special 
orientations of ft). 

• In gravity, M.\ a m — 0, and there are velocity 
fields for which H Q( g = 0, i.e., for which gyro- 
scopes feel no force: 

— in the examples studied, they correspond 
to the situations where B a p = Bt a m in 
the electromagnetic analogue; 

— there are even geodesic motions for 
spinning particles: radial geodesies in 
Schwarzschild, circular geodesies in Kcrr- 
dS spacetimes. 

• Formal analogy between the quadratic scalar 
invariants of R a p 7 s an d F a p are useful to ob- 
tain velocities for which H Q ^ = 0. 



V. TIME PROJECTIONS OF THE FORCES, 
AND WORK DONE ON THE TEST PARTICLE 

A fundamental difference between the gravitational 
and electromagnetic interactions concern the time pro- 
jections of the forces Fq and F^ M in the different frames, 
which we shall explore in this section. We start by ex- 
plaining the meaning of the time projection of a force in 
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a given frame, and its relation with the work done by it 
and the particle's energy. 

Consider a congruence of observers 0(u) with 4- 
velocity u a , and let U a denote the 4- velocity of a test 
particle. The following relation generically holds [5]: 



-u a U a 



(90) 



where v a — U a /j — u a is the velocity of the test particle 
relative to the observers 0(u); in the frame u l = 0, v % is 
the ordinary 3-velocity. The energy of the test particle 
relative to 0(u) is E = —P a u a , and its rate of change 
(the "power equation") 



dE 



-F a u a - P a u a . p U f) , 



(91) 



where F a = DP a /dr denotes the 4-force. Thus we see 
that the variation energy of the particle relative to 0(u) 
consists of two terms: the time projection of F a along 
u a , plus a term depending on the variation of of u a along 
the test particle worldline. The first term is interpreted 
as the rate of work, as measured by 0(u), done by the 
force on the test particle. Using u a — U a /j — v a , we can 
write it as 



F a U 

F a u a = --^+F a v a . 
7 



(92) 



In order to better explain the meaning of this expres- 
sion, let us start by the simplest case, a point particle 
with no internal structure (in our expansion, a monopolc 
particle). In this case the momentum is parallel to 
the 4- velocity, P a — mU a , and its mass is a constant, 
m = mo; hence the force is parallel to the accelera- 
tion: F a = DP a /dr ~ m a a . Since U a is always or- 
thogonal to a a (this is a general statement, which fol- 
lows from the normalization condition U a U a = —1), 
that implies F a U a = 0; i.e., the time projection of F a , 
in the particle's frame, vanishes. Such force is said to 
be spatial with respect to U a , just like the accelera- 
tion. That leads to the expression: —F a u a = F a v a , 
telling us that the time-projection of the 4-D force F a 
is the familiar power (i.e., the rate of work per unit 
of proper time r) transferred to the particle by the 3- 
D force (h u )*F» (see e.g. [35l [49]). Consider also for 
simplicity an inertial frame, so that the second term of 
(91) vanishes; then —F a u a = dE/dr = mod'y/dr, i.e., 



F a v a — mod'y / dr is the rate of variation of kinetic en- 
ergy of translation of the particle's CM. It is clear in 
particular that dE/dr — F a v a — in a frame comoving 
with the test particle, since therein v a = 0. This is the 
type of force we are more familiar with; an example is the 
Lorentz force on a charged particle, DP a /dr — qF a ^Up, 
whose projection along u a reads —u a DP a /dr = r )qv a E a , 
yielding the rate of work (per unit proper time) done by 
the electric field on the particle moving with velocity v a 
relative to 0(u). 

However, if the particle has internal structure, its inter- 
nal degrees of freedom may store energy, which in general 



will be exchanged with the energy of the external fields 
and the kinetic energy of the center of mass. Therefore 
the proper mass of the particle m = —P a U a no longer 
has to be a constant, cf. Sec. HE (For the case of a 



quasi-rigid spinning particle, its internal energy is essen- 
tially kinetic energy of rotation about the center of mass, 
which varies in the presence of an electromagnetic field, 
as explained in Sec. VI below, see also [55TI67] .) Also 



the momentum will not be parallel to U a , as the parti- 
cle in general possesses hidden momentum, cf. Sec. |IID| 
These, together (as we shall see next), endow F a with a 
nonvanishing time projection: F a U a ^ 0. Also the con- 
tribution F a v a to (92) is no longer solely a variation of 



translational kinetic energy, see e.g. Eq. ( 103 ) below 



Let us turn our attention now to the second term of 
Eq. ([91]). Decomposing (e.g. PIUEl]) 



^q;/3 = -~a(u) a Ui3 + LO aj3 + 6 aj3 



(93) 



where a(u) a — u a pu" is the acceleration of 0{u) {not 
the particle's!), uj a p = {h u )a{h u YpU^-. v } the vorticity, 
and 9 a p = {h u )a(h u )pU^) the total shear of the con- 
gruence (9 a p = a a p + 9(h u ) a p/3, being a a /3 as usual 
the traceless shear and 9 the expansion scalar). Let us 
denote by G(u) a — —a(u) a the "gravitoelectric field" 
[9l [34] measured in the frame u l = 0. Decomposing 

pa 



the second term of Eq. (91 1 becomes 



P& d , cf._Eq. (|26|), and U a using (|90|, 



-P a u a ; P U? = mi 2 [G(u) a - 9 a pv?]v a 

+ 7 P h tt id [G(u) a - (cj a0 + 9 a p) V ] 



(94) 



This part of dE/dr depends only on the kinematical 
quantities of the congruence. That is, unlike the term 



(92 ), which arises from the 4-force F a , the term (94) does 



not depend on any physical quantity one can locally mea- 
sure; it is locally an artifact of the reference frame, which 
can always be made to vanish by choosing a locally iner- 
tial one. Its importance (in a non-local sense) should not 
however be overlooked. To understand it, consider a sim- 
ple example, a monopole particle in Kerr spacetime, from 
the point of view of the congruence of static observers 
(i.e., u a parallel to the time-like Killing vector field 
£ = d/dt, in Boyer Lindquist coordinates). Since the con- 
gruence is rigid, 9 a p = 0; also, for a monopole particle, 
P£ id = 0, and, in a gravitational field, F a = (the parti- 
cle moves along a geodesic). Therefore, the energy vari- 
ation reduces to dE/dr = —P a u a .^U^ = m7 2 G(u) Q u Q ; 
which is the rate of "work" (per unit proper time r) done 
by the gravitoelectric "force" 0HGH] m>y 2 G(u) a . (In the 
Newtonian limit, reduces to the work of the Newtonian 
force mG.) Hence we see that (94) is the part of (91) 
that encodes the change in translational kinetic energy 
of a particle (relative to static observers) which occurs 
due to the gravitational field, without the action of any 
(physical, covariant) force, i.e., for particles in geodesic 
motion. 



2G 



Substituting Eqs. (92 1 and (94) in (91), we obtain a 



generalization, for the case of test particles with varying 
m and hidden momentum, of the "power" equation (6.12) 
of [5] (the latter applying to monopole particles only). 



A. Time components in test particle's frame 

One fundamental difference between the tensorial 
structure of H a ^ and B a p is that whereas the former is a 



spatial tensor in both indices: M.(u) a ff' 



,P = 



l{u) a pu a 



(this follows from the symmetries of the Riemann ten- 
sor), the latter is not: B a p(u)u a = but B a p(u)v,P = 
*Faj-.pu 1 uP in general. That means that whereas 
Fq is a spatial force, F^ M has a non vanishing time pro- 
jection in the particle's proper frame, —F£ M U a . Let us 
discuss its physical meaning. First note, from Eq. §4a), 
that 



Fg M U a = B' 3a U a ^ = e 0SliV U s E^fx^ ; 



(95) 



where we see that it consists on a coupling between \i a 
and the space projection of the antisymmetric part of the 
electric tidal tensor E a @ m easur ed in the particle's frame, 
which, as discussed in Sec. II B (see also [SI 02]), encodes 



Faraday's law of induction. Indeed, for an arbitrarily ac- 
celerated frame (non-rotating and non-shearing) we may 
replace E a p by the covariant derivative of the electric 
field E a -p, cf. Eq. (20), and Eq. (95) becomes, in vector 
notation, F^ M U a = — (V xE)-fl. Its significance becomes 
clear if one thinks about a magnetic dipole as a small cur- 
rent loop of area A and magnetic moment /2 = riAI, see 
Fig. [3k. It then follows: 



■F£ M U a = {VxE)-nAI = I dE = 1 E = V\ 



s 1 



ind 



loop 



(96) 

where dE is the 3-D 2-form dE = -Ej-^dx 1 A dx 3 , and E 
is the induced el ectr ic field, coming from the generalized 
induction 14 law (21 ). In the second equality we used the 



fact that the loop is (by definition) infinitesimal, the third 
equality is the application of Stokes theorem in the 3-D 
rest space of the dipole. 

Thus —F^ M U a = Vi n d is the rate of work transferred 
to the dipole by Faraday's law of induction. Using Eqs. 



(31) and (32), we see that it consists in the variation of 



the proper mass m, minus the projection along U a of the 
derivative of the hidden momentum (to which only the 





Figure 3: A magnetic dipole (depicted as a current loop) in 
the inhomogeneous magnetic field of a strong magnet, from 
the point of view of two different frames: a) the particle's rest 
frame; b) the frame of the static observers (at rest relative to 
the strong magnet). /2 = I An (magnetic moment); A = area 
of the loop; I = current through the loop; n = unit vector 
normal to the plane of the loop; E = induced electric field. 
In the dipole's frame there is a non-vanishing work done on 
it by E, at a rate Vind = —F£ M U a , which is reflected in a 
variation of proper mass m. From the point of view of static 
observers u a , the work is zero (—F£ M u a = Vind + Ptrans = 0), 
manifesting that a stationary magnetic field does no work. 
That may be regarded as an exact cancellation between Vind 

and Ptrans- 



"electromagnetic" hidden momentum contributes): 

drn _ DPSd Ua = * 
dr dr dr dr 



Eind = 



n p a 

hidEM tt 



(97) 



Note also, from Eq. (31), that (since m has the inter- 
pretation of the energy of the dipole as measured in its 
proper frame), Vi n d is the variation of the dipole's en- 
ergy E = — Pq as measured in a momentarily comoving 
inertial frame. It should be mentioned that the failure 
to notice that F£ M has a time projection, and that m 
varies, has led to the difficulties in [151 E] to covariantly 
describe the force on a magnetic dipole (namely to the 
claim that no covariant description of such forces is this 
scheme is possible). 



The induction phenomenon in ( 96 ) has no counterpart 
in gravity. Since M a p is a spatial tensor, we always have 



F£U a 







(98) 



which means that the energy of the gyroscope, as mea- 
sured in a momentarily comoving inertial frame, is con- 
stant. The proper mass m is also constant, since, in the 
gravitational case, P a a a = —U a DP{f id /dT is always zero 
(see also Sec. HE). Hence, we see that the spatial char- 



acter of gravitational tidal tensors precludes induction 
effects analogous to the electromagnetic ones. 



This generalization of Maxwell-Faraday law for accelerated 
frames is needed if one is to deal with the electric and magnetic 
fields measured in the test particle's frame, which in general ac- 
celerates. One could instead base the analysis in the inertial 
frame momentarily comoving with it, as done in Sec. V of [6], 
where dB/dr = —V X E holds; the two treatments are equiva- 
lent. 



B. Time components as measured by static 
observers 

With respect to an arbitrary congruence of observers 
u a , the time projection of the force exerted on a magnetic 
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dipole is, cf. Eq. (92 1: 



F EM U /3 
7 



ind 

7 



fEM"« (99) 



where, in accordance with discussion above, we identify 
Vi n d = — as * ne P ower transferred to the dipole 

by Faraday's induction, and F^ M v a is the power trans- 
ferred by the 3-force (h )^F^ M exerted on it. Consider 
now a congruence of observers u a for which the fields are 
covariantly constant F a ^ 1 u 1 = 0, which we dub in this 
context " static observers" 15 . The time projection of the 
force in this frame vanishes: 

D P 

~ (F E M) a u a = -^u a = *F a0;1 U^ a u"< = 0. (100) 

UT 

That tells us that the total work done on the dipole, as 
measured in this frame, is zero. The scalar E = —P a u a 
(the energy of the particle in the u l = frame) is thus a 
conserved quantity, and using Eq. ( 26 ) , we can write it 
in the form 



E 



T + E\ 



hid 



constant, 



(101) 



where we dub E^id = ~Phid u a as the "hidden energy" 
(i.e., the time component of the hidden momentum in 
the u l = frame), and T = (7 — 1)to is the kinetic en- 
ergy of translation of the center of mass (cf. e.g. [H] p. 
70), as measured in this frame. The reason for the later 
denomination is seen taking the Newtonian regime, where 
T m mv 2 /2. In the (very few, to the authors' knowledge) 
existing literature addressing this problem, the fact that 
the total work done on the particle is zero is in [39] justi- 
fied with an exact cancellation of dT/dr and dm/dr; or 
in [55HBT] . for the case of a spherical spinning charged 
body, between the variations of T and kinetic energy of 
rotation about the CM (as shown in Sec. VI A 2 for such 
body and if PuL = 0, dm/dr is a variation of kinetic en- 
ergy of rotation). Eq. (101) shows however that, in the 



general case, the energy exchange occurs between three 
parts, Ehid also having a role. A suggestive example are 
the bobbings of a particle with magnetic dipole moment 
orbiting a cylindrical charge considered in Sec. III.B.l of 
[29] (see also [25]). 

In this work we are especially interested in the case: 
Phid = -^i«d = 0); m this case, m + T = constant, 
i.e., the energy exchange, due to the action of the force 



15 The reason for this denomination is the fact that the condition 
F a ^ 1 u 1 = corresponds to the observers at rest with respect 
to the sources in the electromagnetic applications herein (the 
coulo mb charge in Sec. |IVX| the spinning charge of Sees. |III| and 
|IV B| or the magnet in Fig. |3b)). Note that for e.g. observers u ,OL 
in circular orbits aro und a Coulomb charge we have F a ?yu'"> ^ 
(see expressions | |57| and discussion therein), even though u' a is 
in that case a symmetry of F a p: C u iF a p = 0, and the latter is 
explicitly time-independent in a rotating frame. 



Fg M , occurs only between proper mass (i.e., internal en- 
ergy) and translational kinetic energy. Also, "Pmd = 
dm/r, and for static observers in flat spacctimc u a .p = 0, 
so that the second term of (91) vanishes. We can thus 
write: 



dE 
dr~ 



'FftM u a 



where 



dT 
dr 



1 

1 

7 



dm 



o 



F, 



EM u a 



(102) 



(103) 



is the rate of variation of translational kinetic energy, and 
the the third equality may be obtained by differentiating 
( 101 ) and comparing with (99 1 . An example is the prob- 
lem depicted in Fig. [3|t 



a magnetic dipole falling along 
the symmetry axis of the field generated by a strong mag- 
net. Pj" d = for this configuration. That this is a solu- 
tion (a non-helical one) of the equations of motion sup- 
plemented with Mathisson-Pirani condition is shown by 
arguments analogous to the ones given in Appendix |C 1| 
for the gravitational analogue depicted in Fig. [J below. 
From the point of view of the static observers, E = and 
only magnetic field B is present; since the latter does no 
work in any charge/current distribution, no work can be 
done on the dip ole; thus naturally F£ M u a = 0. Accord- 
ing to Eq. (102), this arises from an exact cancellation 



between Ptrans and Vi n d- on the one hand there is (in the 
u % = frame) an attractive spatial force F^ M causing the 
dipole to gain translational kinetic energy; on the other 
hand there is a variation is of its internal energy (proper 
mass to) by induction, which allows for the total work to 
vanish (this is in agreement with rationale in [33], p. 21). 
See also in this respect Sec. |VI C| and Appendix |B 4| 

In gravity, since those induction effects are absent, cf. 
Eq. (98), such cancellation does not occur. For arbitrary 
u a we have thus 



Fgu a = Fgv a ? 



(104) 



This implies in particular that, from the point of view 
of static observers, by contrast with its electromagnetic 
counterpart, Fq does work on the test particle (i.e., a 
stationary gravitomagnetic tidal field does work on mass 
currents). Let us discuss the implications. A conserved 
quantity for a spinning particle in a stationary spacetime 
is (e.g. [22J I2HH1 OS]) 



1 



a/3 



constant 



(105) 



where £ = d/dt is the time-like Killing vector field. Con- 
sider the congruence of the static observers 16 0{u), de- 



See point [8] of Sec. |l D| In stationary asymptotically flat space- 
times, such as the Kerr metric studied below, these are observers 
rigidly fixed to the asymptotic inertial rest frame of the source. 
They are thus the closest analogue of the flat spacetime notion 
of observers at rest relative to the source in the electromagnetic 
systems above. 
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fined as the unit-time like vectors u a tangent to we 
may write = £u a , with £ = \J — (the lapse, or 
rcdshift factor, of the static observers, see e.g [1]). The 
first term of (jl05h , ~P a £, a — Et; is the "Killing energy", 



a conserved quantity for the case of a non-spinning parti- 
cle (S a P = 0) in geodesic motion, which yields its energy 
E with respect to the static observer at infinity 17 . It can 
be interpreted as the "total energy" (rest mass + kinetic 
+ "Newtonian potential energy") of a monopole particle 
in the gravitational held (e.g. [82] ). The energy £/tot can 
likewise be interpreted as the energy at infinity for the 
case of a spinning p artic le. To see the interpretation of 
the second term in (|105|, V = £[ a;i g]5 Q,9 /2, consider the 



case that P{? id = 0. 







dE t , 



(It 



We have 

-Ea^a 

dV 
dr 7 



dV 

~dr 



■3) 



(106) 
(107) 
= 0. 



where we used the Killing equation £V 
—£,F£u a = £,FqV is the rate work (per unit of particle's 
proper time r) of Fq, as measured by the static observers 
at infinity, and thus V is the spin-curvature potential en- 
ergy associated with that work 18 . In order to compare 
with the electromagnetic Eq. (101 1, note that d^/dr = 
--fG(u) a v a , and that for Pg m = 0, E = 7m = m + T, 
thus we can expand dE to t/dT = d(£E + V)dr in the form 



'dr 



£mj 2 G a v c 



dV 
~dr 



= 0. 



(108) 



The second term accounts for the "power" of the grav- 
itoelectric "force" m7 2 G(u) (which, as explained above, 
is not a physical force, it arises from the acceleration of 
the observer congruence, being non-zero even for geodesic 
motion); in the weak held limit this is the variation of 
Newtonian potential energy. Eq. ( 108 ) tells us that the 



variation of translational kinetic energy T comes from 
the potential energy V, and the power of G(u) (being m 
constant); this contrasts with the case of the magnetic 
dipole discussed above, where (also for P£ id — 0) the ki- 
netic energy comes from the variation of proper mass m, 
no potential energy being involved. In terms of the work 
done on the particle, Fq is thus more similar to the elec- 
tromagnetic forces exerted on a monopole charge or on 
an electric dipole (for Dd a /dr — 0), see Appendix B4 



where the proper mass is as well constant and energy 
exchange is between T and potential energy. 



17 If the particle is in an orbit with always finite r, one can imagine 
this measurement process as follows: let -Etot(Ti) be the total 
energy of the particle at n; this means that if, at that instant, 
the particle was by some process converted into light, and sent 
to infinity, that radiation would reach infinity with an energy 
E = -u a P a = £tot(Ti). 

18 One may check explicitly that dV/dr = i a -p~ l S a ^U^ = £Fgv a , 
noting that DS a @ /dr = if P{f id = 0, and using the general 
relation for a Killing vector £ m; „a = Rxaav^" ■ 



Gravitational 




Figure 4: a) Gyroscope (small Kerr black hole) in the field of 
a large Kerr black hole; b) black hole merger. Evidence that, 
unlike its electromagnetic counterpart, gravitomagnetic tidal 
field does work: spin-dependent part of the energy released is 
the work (as measured by the static observers at infinity) of 
FS. 



There is a known consequence of the fact that Fq does 
work (and of the interaction energy V): the spin depen- 
dence of the upper bounds for the energy released by 
gravitational radiation (Fig. |4b) when two black holes 
collide, obtained by Hawking ]30j from the area law. 

In order to see that, consider the apparatus in Fig. [4] 
two Kerr black holes with spins aligned, a large one (mass 
M, spin J = aM) which is our source, and small one (4- 
vclocity U a , spin S = y / S a S a ) which we take to be the 
test particle, falling into the former along the symmetry 
axis (How this is setup with Mathisson-Pirani condition 
is discussed in Appendix CI). For axial fall, and given 

0. (Consider 

e7& 



Miidl 



that S is also along the axis, P3 d 
again the frame of the static observers 0(u), u a 
i.e. the observers with zero 3- velocity in Boyer-Lindquist 
coordinates). For this configuration, V is a pure spin- 
spin potential energy (that one might also conclude from 
Sec. IV A| where we have seen that for J — 0, Fq = 0); 
it reads: 



V(r) 



± 



2aMSr 
(a 2 + r 2 Y x 



-(r) 



ZFgUcdT 



the +/— sign applying to the case that S and J are par- 
allel/antiparallel. The second equality follows directly 
from Eq. (107), and can also be easily checked noting 
that, in Boyer-Lindquist coordinates, £,FqU = (Fg)o, 
and computing explicitly the time component (Fg)o for 
axial fall, Eq. (37) of PQ. Thus we see that V(r) is minus 
the work done by Fq as the particle goes from infinity 
to r. Let us comment on the presence of the lapse fac- 
tor £ in the integral above. Computing the work of Fq 
does not amount to integrate the power measured by the 



local static observers, Eq. (92): —Fgua = FqV (i.e., 
to sum up the work elements dW = FQV a dr), as that 
would mean summing up energies measured by different 
observers; but to integrate instead the quantity t;FQV a , 
which can be thought as summing work elements mea- 
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sured by static the observer at infinity. Let us now ana- 
lyze the problem of the black hole merger. The increase 
of translational kinetic energy of the small black hole dur- 



ing the fall is given by Eq. ( 108 ) above (or by the power 
Eq. (91); note in this case that 6 a p — 0, since 0(u) is a 
The second term of (108) is the gain 

and 



rigid congruence). 

in kinetic energy due to the "Newtonian" attraction, 



exists regardless of S a ; V by its turn is a spin-spin en- 
ergy; thus the kinetic energy of the particle (and therefore 
energy available to be released by gravitational radiation 
in the black hole merger) depends on S. Upper bounds 
for this energy which are, accordingly, spin dependent, 
were obtained in [3D] by a totally independent method. 
From these limits, and for the case of the setup in Fig. 
[4j Wald [T] obtained an expression (Eq. (35) therein) for 
the amount of energy AE S by which the upper bound is 
increased/reduced when S is parallel/antiparallel to J, 
comparing with the case S = (fall along a geodesic). 
This energy is precisely minus the value of V{r) at the 
horizon r + : AE S = — U(r + ); that is, it is the work done 
by Fq on the small black hole as it comes from infinity 

to the horizon: AE S = f^ r+) (-£Fgu a )d,T. 

We close this section with some remarks on the mean- 
ing of the work done by the gravitomagnetic tidal field. 
To the frame of the static observers in the Kerr space- 
time one ca n as sociate a gravitomagnetic "vector" field 
H (see Sec. Eel and e.g. [1 O |M] ; 



_JJ and e.g. [1 U M EH] 5 in the weak field 
regime this field is well known to be very similar to the 
electromagnetic counterpart, e.g. 2, 3, 83 ) which causes 
inertial (i.e., fictitious) accelerations in test particles of 
the type v x H, formally similar to the magnetic force 
v x B. Namely the "force" is orthogonal to the veloc- 
ity; hence this analogy might lead one to conclude that, 
similarly to its electromagnetic counterpart, the gravito- 
magnetic field can do no work on test particles, i) One 
must bear in mind that H (by contrast with B) has no 
local existence, it is a mere artifact of the reference frame; 
hence it would never be involved in a covariant quantity 
like the 4-force DP a /dT, or the work done by it; ii) both 
in electromagnetism and in gravity, are the tidal fields 
that yield the force; that is manifest in force Eqs. |ljl) of 
Table |U The elect romagnetic tidal tensors are essentially 
derivatives of the fields; for this reason we were able to 
reason in terms of the fields in the applications depicted 
in Fig. [3] (even though are their derivatives that show 
up in the equations). But the gravitational tidal tensors 
cannot be cast as derivatives of the GEM fields, even in 
the weak field regime, except for very special conditions 
(as shown in Sec III see also 134] [57]); the force Fg 
is thus in general very different from its electromagnetic 
counterpart. Namely it is so whenever the test particle 
moves relative to the source — so that the work of Fq 
can dramatically differ from that of F£ M , which is well 
exemplified by the contrast herein: in the test particle's 
frame, we have (Fem)o 7^ 0, {Fc)o = 0; in the frame of 
the static observers, we have precisely the opposite situ- 
ation: (Fem)o = 0, (F G ) ^ 0. 



Sec. in brief — the work done on the particle 
(magnetic dipole vs gyroscope) 

• The time projection of the force, —F a u a , is the 
rate at which it does work on the particle, as 
measured by the observer of 4- velocity u a . 

Time projections along the particle's worldline (U a ) 



• The electromagnetic one is non-vanishing 
F% M U a 7^ 0; it is the rate of work done 
by Faraday's induction, arising from E\ a fs\ (or 
equivalently, from B a pU°)] reflected in a vari- 
ation of m. 

• The gravitational one is zero, FqUo, = 0; the 
gyroscope's proper mass m is constant; no ana- 
logue of Faraday's induction (as H. a pljP = 0). 

Time projections relative to static observers (u a ) 

• The electromagnetic one is zero, F^ M u a = 0; 
a stationary electromagnetic field does no work 
on magnetic dipoles. 

• The gravitational one is non-zero, F^u a 7^ 0; 
gravitomagnetic tidal field does work — there 
is a spin-curvature potential energy; explains 
Hawking- Wald spin-spin interaction energy. 



VI. BEYOND POLE-DIPOLE — THE TORQUE 
ON THE SPINNING PARTICLE 



In our pole-dipole approximation, as we have seen in 
Sec. IIC it follows, from Eqs. ^ or (23), that for purely 



magnetic dipoles (d a = 0) if ft = aS, S is a constant 
of the motion. This might be somewhat surprising. If 
one imagines the magnetic dipole as a spinning charged 
body, one would expect that, in a time- varying magnetic 
field, the induced electric field will in general exert a net 
torque on it, which will accelerate 19 the rotation of the 
body. But in Eq. (23) we only find the term ft, x B, 



coupling the field to the dipole moment (which is there 
in any case, even if the field is constant), and no term 
coupling to the derivatives of the electromagnetic fields; 
i.e., no trace of induction phenomena. Moreover we have 
seen in Sees. Ill El and [V] that the induced electric field 
does work on the spinning particle, causing a variation 
dm/dr = —fl-DB/dr of its proper mass m; that variation 
is shown in the non-relativistic treatments in [65 , 67J (for 



19 Unlike the torque due to the magnetic field, the torque due the 
induced electric field will not be orthogonal to S, and hence will 
in general change its magnitude. For instance, in the application 
in this section, i?i n( j has circular lines around 5, so that Tj n( j || S. 
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the case of a "rigid", spherical body; it is not cast therein 
as a variation of proper mass m, but of the Hamiltonian 
contribution —fi-B) to be but a variation of its rotational 
kinetic energy. Thus we expect it as well to be associated 
to a variation of the spinning angular velocity, and hence 
of S. 

As we shall see below, this apparent inconsistency is an 
artifact inherent to the pole-dipole approximation, where 
terms 0(R 2 ) (being R the size of the particle), which are 
of quadrupole type, are neglected; indeed, whereas the 
contribution of induction to the body's energy is of the 
type p ■ B, i.e., of dipole order, the associated torque in- 
volves the trace of the quadrupole moment of the charge 
distribution. And there is no analogous torque in the 
gravitational case, confirming the absence of an analo- 
gous gravitational effect. 

For clarity, we will treat the two interactions (elec- 
tromagnetic and gravitational) separately, and consider 
first, in Minkowski spacetime, a spinning charged body 
subject to an external electromagnetic field, and then in 
a gravitational field. 



A. Electromagnetic torque 

We will start by the electromagnetic case in flat space- 
time. The equation for the spin evolution of an extended 
spinning charged body is, up to quadrupole order [21] : 



D(S C 



\a/3 



eh 



(109) 



where (P D ix) a and {S clin ) a P are defined by Eqs. (A4) 



and ( A5 ) and consist on the sum of the physical mo- 
menta P a , S a ^ plus electromagnetic terms P' a , S a ^; 
see Appendix |a} (S can ) al3 and (Poix) Q + qA a are shown 
in [S3] to be the canonical momenta associated to the 
Lagrangian of the system. Q a @ is the electromagnetic 
dipole moment defined in m 01 ^ 1 = AQ^ 01 ^ 1 /3, being 
Q a P~t an electromagnetic quadrupole moment defined as: 



Q 



a/37 



J' 



a[/3 7 ] 



(110) 



where J a ^ v and q av are, respectively, the current and 
charge "quadrupole" 20 moments (cf. Eqs. (3.8)-(3.9) of 



20 Following the convention in some literature 29 , 85], we are de- 
nominating expressions of the type | |111| -| [TT2] | as "quadrupole 
moments". Note however that it is also customary in the litera- 
ture to denominate by charge quadrupole moment the traceless 
part of q a P . The former includes the trace and is non-zero even 
for a uniform chimerical body; the latter definition, more consis- 
tent with the actual picture of a quadrupole of charges, measures 
a certain type of deviation from spherical symmetry (and is zero 
for an uniform spherical body). q a @ may also be called the sec- 
ond moment of the charge, according to nomenclature in e.g. 
[31], p. 977. 



E(r,C7) 



J° 



(111) 
(112) 



S(r,C/) 



The vector u; 7 , defined as a vector such that displacement 
of every point by w 1 dT mapsS(r) into E(r+dr), is given 
in flat spacetime by Eq. (All; it can be taken in def. 
(112) above as w 1 ~ J7 7 , since in this approximation we 
are neglecting terms up to r 3 . Hence we may write: 



(113) 



Using Eq. ( A8 i) we may re- write Eq. ( |109 ) explicitly 
in terms of the physical angular momentum S a/3 : 



DS afi 
dr 



£>(<Scan)^ 

dr 



DS' afi 



S' af) = F^q^ (114) 



Note that DS' a/3 /dr is a quadrupole type contribution. 

We are interested in the torques t°, i.e., the vector 
that measures the rate of deviation of the spin vector 
from Fermi-Walker transport: 



D F S a 
dr 



1 



t af) 



Us 



DS af> 
dr 



thus using Eqs. (109 1, (114) it follows: 
D F S a 



dr 

'DEM 
'QEM 
'QEMcan 



' DEM ' ' QEM 



'QEMcan 



1 



A 



2 a 
1 



EM (q\ - 
a Dq*r 



dr 



(115) 

(116) 

(117) 
(118) 
(119) 

(120) 



t dem i s the dipole torque already present in Eq. J9|, 
i.e., just a covariant form for the familiar torque r = 
jl x B + d x E, where no trace of electromagnetic in- 
duction is found. We split the quadrupole torque Tq EM 
in two parts: 7Q EMcan , that we may dub the "canonical 
EM quadrupole torque", is the torque 21 coming from the 
third term of (109) (i..e. the electromagnetic quadrupole 
contribution to D(S can ) al3 /dr). The other part is t" 7 = 



21 Note that Tq EMcan is actually what in most literature is por- 
trayed as the quadrupole torque, see e.g. 1291 . It is clear however 
from Eq. ( | 1 16 [ ) above (written in terms of the physical angular 
momentum S a ^ , S a ), that it is not the total quadrupole torque 
r QEM» ant ^ * ne discussion herein emphasizes how crucial this dis- 
tinction is. 
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he a n 5 U$DS Q/3 /dr, and plays a crucial role in this dis- 



cussion, as the first term of (120) is the torque due to 



the electric field induced in the CM frame, by the gener- 
alized Maxwell- Faraday law (21). This is what we shall 
now see. For that it is convenient to work in the Lorentz 
frame momentarily comoving with the particle's CM. In 
this frame the torque (about the center of mass) due 
to the induced electric field is f| n( j = J p c r x E[ n< id 3 x, 
where p c = ~j a U a is the charge density in the CM 
frame. Let us expand £ in a Taylor series around the 



CM: E l = E l cu + (E, 



^CM) l jr 3 + ... (for the integral above, 
and to quadrupole order only terms to linear order in 
r are to be kept in this expansion), which, splitting 



E 1 * 



E l = E, 



+ (EcM^'^'fj, we may write as 



cm " 2^ x ( V x £ )cm] 



CM 



,(i>i), 



The first term is the part of E that has a curl, and is thus 
the induced electric field: Ei n d(r) ~ —fx (V x _E)cm/2 
(the second part may be cast as a gradient of some scalar 
function and thus not related with induction). Therefore, 
recalling the definition of q a/3 , Eq. (Ill) above: 



1 



(V x E, 



CM 



f Pc [r 



ay 



d 6 x 



l 



(121) 



which, by virtue of relations (20), is a non-covariant form 
for 



'ind 



1 e< 
2 



'^E^ [ 9 « _ sw j (122) 



i.e., the first term of (120). Note that, using Eqs. (|Tj4a) 
of Table |TJ we can also write it as 



2 B ? 



(123) 



It is interesting to consider the special case of a spin- 
ning body with uniform mass and charge density; we may 
write 2<jl a p = (q^(h u )'p —q a p), where I a p is the moment 
of inertia tensor (see e.g. |31j). which we write covari- 
antly as: 

1% = [ J S [^r y {h u )^ - r a r p ] d^ s , (124) 

JY,(t,U) 

being J a — —T a ^Up the mass/energy current (with re- 
spect to the observer U a comoving with the CM), and 
a = q/2m is the classical gyromagnetic ratio. In this case 
we may re- write ( 122 ) as 



'ind 



Rigid spinning charged body 



(125) 



Consider the case that the test particle is a charged, 
"quasi-rigid" body E3] , such that its charge current 



density j, as seen in the CM frame, is j(r) = p c il x r, 
where the vector f2 Q , spatial with respect to U a , denotes 
the angular velocity of the body. tt a is defined as follows: 
let A a be some vector (also spatial wrt U a ) with origin 
at the CM, and co-rotating with the body; then 



D F A a /dr = SlfAP; n aP = ep^M^U" . 
From definitions |7|) and (11), we have: 



(126) 



0jS 



r^j^U a di: 



which, for the rigid body, in the CM frame, becomes: 



or, covariantly, 



(127) 



The rate of work done on this body by the induction 
V 



torque T; 



(122) 



ind' 



T a O 
'ind""' 



is, thus, from Eq. (123), or 



'ind 



n a = -B a u?p a = 



(128) 



That is, we obtained precisely the work Vind of Sec jV A| 
Eqs. (|95h-( 96 1. This is the result we seek: we have just 



proved that the work transferred to the body by Fara- 
day's law of induction, which, to pole-dipole order, is 
manifest in the projection along U a of the dipole force 
Fg M , and on the variation of the proper mass dm/dr, 
is indeed associated to an induction torque, which causes 
S 2 to vary as expected (since r^ d is not orthogonal to S a 
in general). This torque was known to exist from some 
non-rclativistic treatments in the literature [65-67 , deal- 
ing with the special case of spinning spherical charged 
bodies. Only that torque is not manifest to pole-dipole 
order as it involves the second moment of the charge q a p , 
which is of quadrupole order. Note as well that in order 
to assign a spinning particle a moment of inertia and an 
angular velocity, one must go beyond dipole order, as I a p, 
Eq. ( |124[) , involves the "quadrupole" of the mass (mg) a ^, 
Eq. JT49J). But the rate of work it does, T^f^, by con- 
trast, is manifest to dipole order, since, for a rigid body, 
q a p and Q a combine into the magnetic dipole moment 
p a , by virtue of Eq. (fl27|. 



2. Torque on spherical charged body 

In this context, and in view of a comparison with the 
gravitational problem, it is interesting to consider (in 
Minkowski spacetime) the case of an uniform, spherical 
charged body, whose quadrupole moments of j a reduce 
to the trace of q a p, and we expect the total quadrupole 
torque on the particle r a to come essentially from T-^ d . 
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Figure 5: Current density j for a spherical body rotating 
about the z axis, as seen in the CM frame. In this frame 
the spatial components of J afiv vanish: J ijk = 0. All the 
components involving z axe zero: J %3Z = 0, since j z = 0; and 
jizj _ q f rom t ne symmetry with respect to the equatorial 



plane. It is easy to see that: J x 



Jx 2 fd 3 x 







J 1 



fyxj y d 3 x = = J mvx =J V 



J" 



So firstly let us explicitly compute the quadrupole mo- 
ments for this type of body. It is clear that the charge 
quadrupole Eq. (Ill I is such that, in an orthonormal 



frame comoving with the center of mass, its time compo- 
nents are zero: q 00 = 0, and its spatial part reduces to 
its trace: q i] = 5 ij q k k /3. 

Such tensor is covariantly written as: 



(129) 



As for the tensor j7" a ^ 7 , its spatial components J % ^ k 
vanish in the CM frame, see Fig. [5] The only non- 



vanishing components are = Jr l r : >j d 3 x = q % K 

Hence J a ^{h u )» = 0, thus, by virtue of Eq. (ml, 
japu = q a0jju ^ and therefore: 



(130) 



Substituting (129) and (130) in (110 1 



Q a ^ = ^q\{h u ) a ^U^ = ^q T T g a ^U^ . (131) 

Now let us compute the quadrupole torque exerted on 
the test body. From E qs. (|A8[ i) and (A9) of Appendix 
A 1 together with Eq. (|129|, we have 



Bi 



Qr 



(132) 



Substituting fl32| and ( pL3T| ) in Eq. pl9| ), we obtain 

TQEMcan = U^ a F \e a f U 5 = , (133) 



the second equality holding in vacuum (which is the prob- 
lem at hand) by virtue of Maxwell's equations F a ^ = 
47rj Q . This means that Tq EM = r' a . Now in order to 
compute the torque Eq. (120), we must give a law of 
evolution for q a p. Eq. (129) ensures that the body is 



spherical; we also demand dq^/dr — 0, so that it has 
constant size; together these relations imply that q a p is 
Fermi- Walker transported: D F q a p/dT — 0, i.e., it has 
constant components in an orthonormal tetrad comoving 
with the particle's center of mass, as expected. Hence, 
from Eqs. ( 120 ) and ( 122 ), it follows that the quadrupole 



torque reduces to 



'QEM 



'ind 



a/3 1 



"25V 7 : 



'ind 



a/3 



*/3] 



i.e., up to acceleration dependent terms, r, 



QEM 



(134) 
is the 



torque due to the induced electric field. 

For a body with uniform charge and mass densities, it 
follows from Eqs. ( [124] and ( fTTT] ) that qf/3 = = 
a I. where / 



I, 



I yy denotes the moment of 



inertia of the sphere with respect to any axis of rotation 
passing through its center. Thus we may put = 
-ale^E^ = -aIB a p UP. We shall write, in the CM 
frame, and in vector notation, the equation for the spin 



evolution (115) for such body: 



DS 



— T DEM + find ' 



G 



-a x E 



^ - DB o-I _ -> . , 

= uxB-al— axE. 135) 

p dr 2 y J 

to point out that this is the relativistic generalization of 
Eq. (1) of [55], or Eq. (6) of ESI; those non-relativistic 



results therein follow from Eq. ( 135 ) above by replacing 



t — > t, and neglecting the acceleration dependent term. 
Let us now compute the work, rT^, done by the total 



torque, t c 



'DEM 



'QEM 



on the particle. First note 



that, for a "quasi-rigid" body, the relation S a = / Q/3 r^ 
holds [22] ; which, for a uniform spherical body, cf. defi- 



nition (124), becomes: 



s a = m c 



(136) 



Hence, assuming the proportionality fi a = o~S a , it fol- 
lows from Eq. (117) (with d a = 0, which is the prob- 



lem at hand) that the work of the dipole torque is zero: 
rg EM O a = 0. Thus, T a n a = T« EM n a . From Eqs. p27b 



and (129), we have 



1 



(137) 



and therefore, from Eqs. (134) and (28), 
1 



^aB^E /V — T ind^o- — o^hidEM ™- 



Now consider a setup for which the electromagnetic 
hidden momentum is zero: fy^EM = 0j as i s the case of 
the application in Fig . [6k . In this case, T a D, a = TjJ d f2 CT . 
From Eqs. (juHJ) and ( fl36[ ), we have also r CT = IDQf/dr; 
therefore: 



dr 



ljdn 2 

2 ~d~r~ 



'ind 



(138) 
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where in the last equality we used (1281. Itt 2 /2 is the 



body's kinetic ene rgy o f rotation about its CM, see e.g. 
[351133]; hence Eq. (jl3Sh above tells us that for this setup, 



the rate of variation of the body's kinetic energy of ro- 
tation equals the rate of work, as measured in the CM 
frame, done by the dipolc force F^ M on the particle (that 
is, its projection —F£ M U a along the particle's worldlinc). 

implies also, cf. Eq. (28), 
DB a /dr; therefore 



The condition -P^ idEM 



that B^U? 



1 



'hid 



q\B a U^ 



1 



DB° 

~d~T~ 



and we may write, using (137) 

1 dn 2 



-- _ a _ 



(139) 



(140) 



Using D^/dr = aDS a /dr, together with Eqs. ( |116| , 
(fTlTb, ((1341) and ((1391), we can rewrite Eq. (fl40l) as: 



WO 2 

2 ~dT 



d{B a fi a ) crq\ ( dB 



dr 



6 



dr 



a/3 



a a E l3 B r 



(141) 

which is the relativistic generalization of Eq. (10) of 
(the acceleration dependent term, which arises from our 
demand that q a p be Fermi- Walker transported, is absent 
therein). From this equation we see that, for this setup (a 
quasi-rigid, uniform, spherical charged body), the vary- 
ing part of the mass, —B ■ jl, present in the dipole ap- 
proximation, cf. Eq. ( 36 I , is kinetic energy of rotation 



(not potential energy, as claimed in some literature, e.g. 
[2"2l |3"6"I l4"6"p — this establishes, in a relativistic covari- 
ant formulation, and in the context of Dixon's multipole 
approach, the claims in [64H67] . The remaining terms in 



Eq. (141), involving the charge quadrupole q a p, a re not 
manifest in the (dipole order) mass equation ( 35 ) , since 
to that accuracy B a Dfj, a /dT = 0, by virtue of Eq. (23). 



and the distance between them given by y/r a r a . In the 
general case of a curved spacetime, the distance between 
two points is the length of the geodesic connecting them; 
r a is not even (exactly) a vector, and the closest notion to 
a separation vector is the bitensor —a K (x, z), which is the 
vector tangent to the geodesic at z whose length equals 
that of the geodesic connecting the two points [2H [27] . 
Also one must distinguish a coordinate system at x from 
the one at z, as the basis vectors change from point to 
point in a curved spacetime. 

In the previous sections, dealing with pole-dipole ap- 
proximations, we did not face this problem in defining 
the moments because cutting off the multipole expansion 
in the dipole term amounts to keep only terms linear in 
r; since, to first order, r a is a vector (see e.g. Eq. (7) of 
[86)). and \Jr a r a the distance between two points, the use 
of bitensors is by definition redundant. We can always 
define the moments by the same expressions as in flat 
spacetime. And in an approximation up to quadrupole 
order, albeit not true generically, if the gravitational field 
is not too strong, we can still, to a good approximation 
(cf. [22] p. 513) skip the bitensors and define the mo- 
ments by expressions identical to the ones in flat space- 
time (where a K = -r K , 6 kA ^ = g K ^g^ x and H a x ~ S x , 
see [521 155]). For that, and following [1 [SS], one sets 
up a locally nearly Lorentz Frame e a with origin at the 
reference worldline, momentarily comoving with the test 
particle (i.e., eo = U a and the spatial triad spans the 
instantaneous rest space of the center of mass) . In such 
frame r a — x a yields a good approximation to the dis- 
placement vector between a point x a and the reference 
worldline z q (t). 

From Eqs. (5.29), (5.9) and (4.20) of [H], we have, for 
the case F afi =0: 



J° 



B. Gravitational torque 



x/3-yS _ 1 ^ t -,[al3]S _ t <5[a/3] 7 ^ _ Tj[a p P] [-yS] _ jjhpS] [a/3] 

(144) 

where the moments t a ^ lS and p Q/37 , in the frame e a , are 
defined as: 



The equation for the spin evolution of an extended 
spinning body in a gravitational field is, up to quadrupole 
order El M 



DS K> 
dr 



= 2P [r U x] 



(142) 



leading to 
D F S a 



dr 



r QG T QG 



G 



J^R x lffU s . (143) 



The general definition of the gravitational quadrupole 
moment J a ^ s , given by Eqs. (9.12), (9.11) and (9.4) 
of [3S], involves bitensors (er K , Q kX ^ l ' and H a p in the 
notation of [38] ) . Let us briefly comment on that. Only in 
flat spacetime (and Lorentz coordinates) the separation 
vector between two points x a and z a is r a = x a — z a 



r & r^T^w a dT, a , 



S(t,U) 



(145) 



(146) 



Noting that T^ s dZ s = -T^ v U v U s dH d = J^U & dT, s , 
where J a denotes the mass/energy current, we may 
rewrite p"^ 1 as: 



d/37 



r " r Pj*fTj 5 dE s 



£(t,(7) 



(147) 



which has thus the interpretation of the quadrupole mo- 
ment of the mass current, analogous to the electromag- 
netic quadrupole moment of the charge current J a $~< , 
Eq. (|112[). Again we take w a w U a in (|145|) since, to 



quadrupole order, only terms up to r are being kept in 
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the integrands. Therefore t a ^ a (T u ) s a = p a ^U s , and we 
may decompose t al3lS as follows 



p a P~f may also be decomposed as 



6 

(148) 



P 



a/37 



where 



K) 



a/3 



r a rPj~<d^ 



E(r,t/) 



(149) 



is the mass quadrupole (or second moment of the mass, 
see [25] EU 155] a nd Footnote 20 ) , analogous to charge 
quadrupole (111). 



1. Torque on spinning "spherical" body 

Our goal in this section is to consider the gravitational 
analogue of the problem in Sec. |VI A 2 Therein we con- 
sidered a spherical charged body in flat spacetime, whose 
charge quadrupole moment was shown to reduce to its 
trace: q a g — q T T (h u ) a g/3, and the quadrupole of the cur- 
rent to J al3 » = q a a U' y (h u ) al3 /3. We prescribe the analo- 
gous test body for the gravitational problem by demand- 
ing it to have an analogous multipole structure (i.e., its 
"gravitational skeleton" [35 ) to its electromagnetic coun- 
terpart (rather than demanding its shape to be "spher- 
ical", which in a general curved spacetime is not a well 
defined notion. A body with such multipole structure 
will of course be a sphere in the case of flat spacetime; 
and otherwise may be thought as one if the field is not 
too strong) As shown above, the quadrupole moment 



p°<0i j Eq (147), has an analogous definition to J a $~i , 



Eq. (112), only with J a in the place of j a ; hence its 



structure must be (analogously to Eq. (130)) 



; {m Q y T {hF)^W 



(150) 



In a local orthonormal tetrad e a such that <=o = U a (i.e., 
the triad e-i spans the rest space of the center of mass), 
this tensor has similar components to J a ^ as written in 
the CM frame (only with J a in the place of j a ). From 
Eq. (1148 1) we have 



t a ^ s = ^(m Q ) r r (// / ) Q ' 3 t/ 7 [/ 5 + t af)Xa {h u )l{h u )i . 

(151) 



The second term of ( 151 1 is the quadrupole moment of 
T y , which has no electromagnetic analogue (the time 
projection of T a @: —T al3 Up — J a has an electromagnetic 
analogue, which is the current density j Q ; the space part 
of T a/3 has not). For a "quasi-rigid" spinning body, we 

have in the CM frame (e.g. [57]) T~ ij = pu\ ot u{ ot + 
where u" ot denotes the 4-velocity of the (rotating) mass 



element of the body and s lJ are the stresses. For non- 
relativistic rotation speeds v Iot -C 1, s lJ can be neglected 



[87] [88], and thus T ij w pv l Iot vi 



this regime, the second term of ( 151 1 is negligible com- 



therefore, in 



pared to the first one (involving T 00 = p) , It then follows 



ja/3 7 <5 _ -{m Q y r U[ a g^U 5 \ 



(152) 



(in agreement with Eq. (7.31) of [25]). Substituting in 
Eq. (143), we obtain the gravitational torque: 



T QG 







the second equality holding for vacuum [B} IV = 0), 
which (as in the electromagnetic case) is the problem 
at hand. Thus, no gravitational torque is exerted, up 
to quadrupole order, in a spinning spherical body. This 
means that there is no gravitational counterpart to the 
electromagnetic torq ue t ; ° h o f the analogous electromag- 
netic problem in Sec. VI A 2 that is generated (from the 
viewpoint of the particle's frame) by the induced electric 
field. This is the result we expected from the discussion 



in Sec. II B and in [6]. The electromagnetic torque (122) 
comes from the antisymmetric part of E a p; or, equiva- 
lently, from the time projection of B a p, cf. Eq. (123), 
which encode Faraday's law of induction, cf. Eqs. (\20\ 
and (|l]4a). The absence of an analogous effect in gravity 
is thus natural, since the gravito-electric tidal tensor E Q( g 
is symmetric, and the gravitomagnetic tidal tensor H Q ^ 
is spatial, meaning that the dynamical effects which in 
electromagnetism are caused by the curl of the electric 
field E, have no counterpart in the physical gravitational 
forces and torques. 



C. Summarizing with a simple realization 



VI A| and |VIB entirely corroborate 



The results in Sees 
the discussion in Sec. V] (and also Sec. HE); namely the 
manifestation of electromagnetic induction and the ab- 
sence of an analogous phenomenon in the physical grav- 
itational forces and torques. In this context it is inter- 
esting to consider the analogous problems in Fig. [6} a 
spinning spherical charge moving in the field of a strong 
magnet (or another spinning charged body), and a spin- 
ning "spherical" mass moving in Kerr spacetime. Let 
us start by the electromagnetic case. A force F£ M , 
Eq. (|l]la) will be exerted on the particle, causing it to 
move (thereby gaining translational kinetic energy, at a 
rate 7-trans, Eq. (103)). And as it moves in an inhomo- 
geneous magnetic field, a torque r^ d is exerted upon it, 
due, from the the viewpoint of the observer comoving 
with the particle, to the electric field induced by the 
time-varying magnetic field. That torque will cause a 
variation of angular momentum S a , and therefore of the 
angular velocity il a — S a /I of the particle (measured 
with respect to the comoving Fermi- Walker transported 
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Figure 6: a) A spinning, positively charged spherical body 
being pulled by a strong magnet; _Ei n( j = electric field induced 
in the body's CM frame, b) A spinning spherical body falling 
into a Kerr Black hole. As the spinning charge move s in the 
inhomogeneous magnetic field B, a torque r^ d , Eq. (1341, is 
exerted on it due to -Emd, i.e. to the antisymmetric part E\ a m , 
or equivalently, to the time projection B^pU^. That causes 
S = %/ S a S a , and the body's angular velocity Q = S/I, [cf. 
Eqs. (|136[)-|T26|], to vary. T^ d does work at a rate T^" d f2 a = 



■pind, which exactly matches time projection of the dipole force 
Fg M , cf. Eq. ( 128 1. That causes the kinetic energy of rotation 



about the CM to decrease, manifest in a variation of proper 
mass dm/dr, and cancels out the gain in kinetic energy Ptrans, 
so that the total work transfer is zero (cf. Sec.|VB}. In the 
gravitational case no analogous induction effects occur (as 
expected since Er a| gi = H a ^(7' 3 = 0): no torque is exerted on 
the spinning particle, its S, angular velocity f2, and proper 
mass to, are constant; and there is a net work done on it by 
Fq at a rate Vtot = —F§u a , corresponding to an increase of 
translational kinetic energy. 



tet rad, cf. Eq. (1261). Clearly, as we see from Eqs. (135) 
or (134), S 2 is not conserved (as would be the case in a 



pole-dipole approximation, see Sec. II C). The variation 



of the magnitude f2 of the angular velocity also implies a 
variation of rotational kinetic energy of the particle; that 
is reflected in a variation of proper mass dm/dr, as shown 
by Eq. ( 128 1, which is precisely the one we found to exist 
in Sees. II E| and |Vj encoded in the time projection of 
As discussed in Sec. 



the force Ff M 



VB 



if we choose 

a frame corresponding to the "static observers" u a (cf. 
point [8] of Sec. ID and Footnote 15) then F£ M u a = 0, 
and the total energy of the particle, E = —P a ii a , is con- 
served. That means that the rate of variation in kinetic 
translational energy Ptrans (due to the work of F^ M ) is 
exactly canceled out by the variation of rotational ki- 
netic energy Vi n d (due to the work of Ti nc j), ensuring that 
a static magnetic field does no work. 

In the gravitational case, there is also a net force Fq 
on the body, cf. Eq. (|l]lb) of Table [TJ causing it to gain 
kinetic energy at a rate Ptrans = F G Vi. But no torque is 
exerted on it; up to quadrupole order we have: 



(i.e., the spin vector of the spinning spherical mass is 
Fermi- Walker transported) implying also SI = constant. 
This is consistent with the constancy of the proper mass 
(and the fact that Fq is spatial, meaning that no work 
done by induction), because since there is no torque, the 
kinetic energy of rotation is constant. Thus also the gain 
in translational kinetic energy is not canceled out by a 
variation of rotational kinetic energy, and therefore a sta- 
tionary gravitomagnetic field will do a net rate of work 
"Ptrans on the particle. 

We close this section with a few remarks. 1) The dis- 
cussion herein explains the reason why, in the framework 
of the inertial fields from the 1+3 formalism, the analogy 
in Eq. ( 24 ) holds even when the fields are time- varying 



D F S a 
dr 



0; S = constant, 



(while we expected it to break down due to the phe- 
nomenon of electromagnetic induction) — it is due to the 
fact that, by neglecting moments beyond dipole, we are 
in fact neglecting the torque T^ d due to Faraday's induc- 
tion, since the latter involves the "quadrupole" moment 
q a fi of the test particle (unlike the work "Pi n d = f^aP-a 
done by it, and the corresponding variation of proper 
mass, which is also due to electromagnetic induction, but 
already manifest in the pole-dipole approximation, as it 
comes in terms of the magnetic dipole moment \i a ; this is 
the reason why m varies whereas S 2 is a constant in the 
dipole approximation). 2) The application in Fig. |6]illus- 
trates an important aspect of the physical effect of frame 
dragging, and the electromagnetic analogue. For clarity, 
let us consider the case that initially the test bodies have 
no rotation. In the electromagnetic case, Fig. [6^,, as the 
ball falls, it gains an angular velocity (relative to the FW 
transported tetrad) due to the torque r; n d; in the gravi- 
tational case, Fig. [HJd, no such rotation arises. However, 
from the point of view of a frame anchored to the distant 
stars, the angular velocity of the spinning mass indeed in- 
creases as it approaches the black hole. This can be easily 
seen as follows. In the linear limit, and in the frame of 
the static observers, it is well known [2] 131 l34 l 157 1 l89 l fTT2] 
that gravity becomes very similar to electromagnetism. 
In electromagnetism, from the point of view of static ob- 
servers, the torque fj n d does not come form the induced 
electric field (which is zero in that frame, as the fields are 
static therein) , but from the overall effect of the Lorentz 
force dqv x B applied to each charge element dq of the 
spinning ball. The gravitomagnetic "force" v x H leads 
to an entirely analogous "torque". These, however, are 
pure artifacts of the reference frame, with no local phys- 
ical existence; they can only be measured by locking the 
frame to the distant stars (by means of, e.g., a telescope). 
The ball in Fig. [6]d, if initially H a = 0, will never gain 
any rotation relative to the local compass of inertia; in 
fact, an observer sitting firmly with his tetrad on top of 
the ball will not detect any sign of rotation: he will not 
measure any Coriolis forces acting on the test particles 
he launches, and sees gyroscopes axes fixed. That, of 
course, is consistent with the constancy of S a . In fact, 
it is the frame of the static observers that rotates rela- 
tive to the compass of inertia, which is manifest in the 



3G 



fact that they measure a non-zero gravitomagnetic field 
H (i.e., they see test particles in geodesic motion being 
deflected by the Coriolis "force" v X H , and gyroscopes 
"precessing"). That is, we have a gravito-electromagnetic 
analogy based on GEM "acceleration fields" (namely H) 
in the frame frame anchored to the distant stars; but 
locally it vanishes, and this is the essence of the frame 
dragging effect. The electromagnetic counterparts, by 
contrast, and electromagnetic induction, are locally mea- 
surable quantities, which manifest themselves covariantly 
in the physical tidal forces and torques. 



ticle 



Sec. | V7| m brief — the torque on the spinning par- 
manifestation of Faraday's induction) 



• The torque due to Faraday's law of induction, 
i ind , comes from a coupling of Ei a m to q a p (the 



re- 



charge "quadrupole") : 

— dipole approximation ignores q a p, as well 
as the moment of inertia; that explains 
why Tf nd is not manifest to dipole order; 

— but the work it does, T" nd fl a , is mani- 
fest to dipole order! It equals the pro- 
jection —F£ M U a of the dipole force along 
its worldline. 



T ind nas n0 gravitational analogue (consistent 
with E Q(3 = 0). 

A time- varying electromagnetic field torques a 
spherical charged body; the gravitational field 
never torques a "spherical" body; m (includes 
kinetic energy of rotation) , 5 1 , and angular ve- 
locity relative to FW transport, are constant. 



VII. CONCLUSION 

In this work we studied the dynamics of spinning test 
particles in general relativity, in the framework of exact 
gravito-electromagnetic analogies. A detailed summary 
of the main results and realizations is given in Sec. [TJ so 
we conclude with some additional remarks. 

Both the equations of motion — force and spin evo- 
lution — of a spinning particle in a gravitational field 
are related to their electromagnetic counterparts by ex- 
act analogies, valid for generic fields. And expanding the 
force equation explicitly as an equation for the center of 
mass acceleration, another exact analogy emerges, for the 
"hidden momentum". The first one (the analogy based 
on tidal tensors for the force) was introduced by one of 
the authors in jS]; the second one (the analogy based on 
the 1+3 incrtial fields for the spin equation) is also still 
relatively little known (to this degree of generality can 
be found only in [HI 02], to the authors' knowledge); the 



hidden momentum analogy, Eq. (29 1, is presented here, 
in its exact form 22 , for the first time. They are all herein 
extracted from the rigorous equations of motion for pole- 
dipole particles (e.g. J2S [23 HH1 ESI [3HI I3U H31 1101 IHI] ) , 
from which they emerge if the Mathisson-Pirani spin 
condition is employed. 

These analogies are useful essentially for two reasons: 
practical purposes, as they provide a familiar formal- 
ism to treat exactly otherwise complicated gravitational 
problems; and theoretical purposes, as they provide suit- 
able tools for a transparent comparison of the gravita- 
tional and electromagnetic interactions (this is especially 
the case of the tidal tensor analogy, which based on physi- 
cal forces, common to both theories). It is the main point 
of this work that there is a lot to be learned (about both 
of them) from a comparative study of the gravitational 
and electromagnetic interactions. 

The first remark we want to make, is that it is impor- 
tant to realize that the existence of these exact analogies 
does not mean that the interactions are similar. The 
GEM analogies are best known in the context of lin- 
earized theory, e.g. [THS1 [52]; not that well known are 
the exact ones, as explained above. These however have a 
different status from the analogies drawn in the linearized 
theory, in the sense that the latter, in order to hold, re- 
quire a degree of similarity between the interactions to 
which the exact ones mentioned above are not bound. 
These are functional analogies: the magnetic tidal ten- 
sor of electromagnetism B ai plays in the equation (jlji) 
for the force exerted on a magnetic dipole the same role 
as the gravito-magnetic tidal tensor H Q7 in the exact 



Eq. |I|b) for the gravitational force exerted on a gyro- 
scope: both forces are determined by a contraction of the 
spin/magnetic dipole 4- vector with a magnetic- type tidal 
tensor as measured by the test particle of 4- velocity U a . 
And the functional analogy extends to Maxwell and Ein- 
stein field equations, as manifest Table [TJ Also, in the ap- 
propriate frame (unlike the former, which is covariant and 
fully general, this analogy is less general, not covariant, 
attached to a special reference frame, as explained in de- 
tail in [33]) the gravitomagnetic field H plays in the "pre- 
cession" of the gyroscope an analogous role to B in the 



precession of a magnetic dipole, cf. Eq. (24 1. The analogy 



extends, under special conditions in stationary fields, to 
the equations for the geodesies, for the force on the test 
particle, and in the field equations, see e.g. [5[ |H1 133] • 
But the analogies do not imply, even in seemingly anal- 
ogous setups (for instance, the electromagnetic field of 
spinning charge, and the gravitational field of a spinning 



22 In an approximate form it was previously obta ined in 
the spin condition S a ^ Pp 



C3b 



[29J using 
Note also 



0; see Appendix 
that the formal similarity of the hidden momentum term S Xi! 
with the electromagnetic field momentum — ft X E of Eq. (D5i 



(opposite to the hidden momentum of electromagnetic origin) 
was already noticed in [5]. 
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mass, considered in Sees. Ill and IV B I that these ob- 



jects are similar. Firstly, the gravitomagnetic field H 
and tidal field H Q ^, unlike their electromagnetic coun- 
terparts, are non-linear in the potentials (except for very 
special cases, such as ultrastationary spacetimes, studied 
in e.g. Sec. IV of 6]). But that is not all; even if one 
goes to the weak field regime (where the non-linearities 
of the gravitational filed can be neglected), the symme- 
tries and the time projections of the tidal tensors con- 
tinue to differ crucially. This shows how misleading can 
be the apparent similarity suggested by the usual linear 
approaches in the literature, which is detailed in Sec. |III| 
Indeed we have seen that the electromagnetic and grav- 
itational effects can be completely opposite: as shown in 
Sec|V]for an observer comoving with the test particle, the 
time component of Fq is zero, whereas F£ M is non-zero, 
and corresponds to the power transferred to the dipole 
by induction; for static observers u a , F-£ M u a = (conse- 
quence of the fact that stationary electromagnetic fields 
can do no work on a magnetic dipole) whereas FqUq ^ 0, 
meaning that a gravitomagnetic tidal field does work on 
mass currents. 

The exact analogies are suited instead for a compar- 
ison between the interactions: as it amounts to com- 
pare mathematical objects that play analogous dynam- 
ical roles in both theories. Such comparison unveils 
suggestive similarities, useful in terms of the intuition 
they provide. But, and especially in the case of the 
tidal tensor formalism, it is actually through the dif- 
ferences it makes transparent that it proves illuminat- 
ing. The differences between the tensorial structure of 
the gravitational and electromagnetic tidal tensors are 
related to fundamental differences between the interac- 
tions, namely the phenomenon of electromagnetic induc- 
tion, and the way it manifests itself in the electromagnetic 
tidal forces and torques, which has no analogue in gravity. 
The results in Sec. |Vj concerning the time components 
of the force, and in Sec. |VI| concerning the torque ex- 
erted on the spinning particle, are manifestations of the 
absence of a gravitational counterpart to the antisym- 
metric part of E a p (or, equivalently, to the projection of 
B a p along U a ); Ei a m encodes the Maxwell- Faraday law 

V x E = —dB/dt; the gravito-electric tidal tensor by con- 
trast is symmetric: E[ a| gj = 0, translating in an absence of 
analogous induction effects in the physical gravitational 
forces and torques. And the results in Sees. |IV A| and 
IV B| showing that in a non-homogeneous gravitational 
field there are moving observers for which M a fj — 0, so 
that gyroscopes can actually move along radial or circu- 
lar geodesies (in Schwarzschild and Kerr-dS spacetimes, 
respectively) , manifest that there is no gravitational ana- 
logue to Maxwell Eq. V x B = dE/dt; as in the electro- 
magnetic systems, due to this law (more precisely, in co- 
variant form 2_B[ Q ^] = *i< 1 a ^ ;7 {7 7 ), B a p is non-vanishing 
whenever the dipole "sees" a varying field (as is the case 
when the particle moves in a non-homogeneous field) , and 
therefore (except for some special orientations of p) an 
electromagnetic force F$ M is exerted on it. 



This contrast is interesting, of course, from the theoret- 
ical point of view, as the deep connection between the two 
theories is still very much an open question; but it is also 
valuable in terms of the intuition it provides to practical 
applications. Acknowledging the differences in the sym- 
metries of the tidal tensors, namely that there is no grav- 
itational counterpart to the electromagnetic antisymmet- 
ric parts, allowed us to predict, in those conditions where 
the electromagnetic force Eq. (JTJa) comes exclusively from 
B[ a /3] > that the force (jljb) would be zero in the analogous 
gravitational setup. In this framework we predicted (Sec. 
IV B 3[ ) the existence of circular geodesies for spinning 
particles in Kerr-dS (which, to our knowledge, were not 
known in the literature 23 ), or that a gyroscope dropped 
from rest in Schwarzschild spacetime would move along 
a geodesic towards the singularity (Sec. IV B). Note that 



even the latter, which is the simplest application in this 
work, could be a complex problem (involving possibly 
complicated descriptions, and difficulties in setting up its 



initial conditions, see Appendix C 1 ) outside the tidal ten- 
sor formalism and the Mathisson-Pirani condition. We 
have also understood (Sees. VA and HE I the conser- 



vation of proper mass of a gyroscope from the spatial 
character of HP^, as in the electromagnetic analogue the 
variation of proper mass arises from the work Pmd of the 
induced electric field, encoded in the projection of B a p 
along the particle's 4-velocity U a (or equivalently, from 
the anti-symmetric part of E a p), which has no gravita- 
tional analogue. In the same spirit we physically inter- 
preted the spin-dependence of Hawking's upper bound 
[5U] for the energy released when two black holes col- 
lide: if one considers (from the point of view of static 
observers) a magnetic dipole falling into a strong magnet 
[as illustrated in Fig. [3]d)], there is no net gain in energy 
of the dipole, as the stationary magnetic field can do no 
work on it. That comes from an exact cancellation (since 
the hidden momentum term vanishes in this case) be- 
tween the loss in proper mass T^nd = dm/dr, which con- 
sist on the rate of work done by the electric field induced 
in the CM frame, and the gain in translational kinetic 
energy "Ptrans- Since V[ n d has no gravitational counter- 
part, such cancellation cannot occur, and therefore we 
expect a stationary gravitomagnetic tidal field to do a 
net work —F^Ua = F^Va on a "gravitomagnetic dipole" 
(a gyroscope); which turns out to be the case: it leads 
the to the Hawking- Wald spin interaction energy [T]. In 
other words: we show that the gravitational spin interac- 
tion energy, and the spin dependence of Hawking's upper 
bound for the collision energy, are justified with the fact 
that, unlike its electromagnetic counterpart, a stationary 
gravitomagnetic tidal field does work. 

Its is worth noting also that, on some conditions where 



23 In the context of spinning particles in Kerr-dS spacetimes, it 
is worth mentioning the work [94] . where another problem was 
addressed: equilibrium positions. 
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the forces differ significantly (for instance when the parti- 
cle moves relative to the source), certain qualitative sim- 
ilarities remain, from which one can still get useful intu- 
ition to visualize gravitational effects. An example is the 
deviation from geodesic motion of a gyroscope (in non- 
radial motion) in Schwarzschild spacetime 24 . Eqs. (|ljl) 
show that in both cases it is the magnetic tidal tensor as 
seen by the test particle that determines the force. Hence 
the gyroscope deviates from geodesic motion by the same 
reason that a magnetic dipole suffers a force even in the 
coulomb field of a point charge: since it is moving in the 
electric field, in its "rest" frame, there is a non-vanishing 
magnetic tidal tensor (one must however bear in mind 
that in this case the two tidal tensors can be very dif- 
ferent; and that whereas in the case of a moving dipole 
one has always B a p ^ 0, the same does not necessarily 
happens in the gravitational case; see Sec. IV A I. 



netic parts, and their scalar invariants) we use in 



As for the other physical analogies, based on "iner- 
tial" fields — the analogy for spin precession of Sec. |IIC[ 
and the hidden momentum analogy of Sec. IIP the first, 
though not extensively used in this work, still plays an 
important role, as one of the results in Sec. |VI| is the 
clarification of the reason why it holds for arbitrary fields 
(while one might expect it to break down for time depen- 
dent fields, due to electromagnetic induction, see Sees. 



IIC and VIC). It is also involved in the discussion in 



Sec. IV B 2 of the gravitational analogue to the velocity 
field for which a magnetic dipole does not precess. This 
analogy, in the context of weak, stationary fields, is one 
of the best known (yet its precise mathematical meaning 
not often well understood; see [33]) and widely used, as 
it can be a applied to effects such as Thomas-Precession 
EE CTJ or the Lense-Thirring effect IS EE ISO [53] . 
The hidden momentum analogy, by its turn, is used in 
the companion paper [25] to explain the dynamics of the 
famous Mathisson helical motions, and in [29j to explain 
the bobbings observed in numerical simulations of binary 
systems (though therein not using the exact analogy in 
Eq. ( 29 ) , but an approximate one obtained by employing 
Tulczyjew-Dixon spin condition instead). In Appendix 
|D 1| we interpret the formal similarity between the "iner- 
tial" S x G and electromagnetic fix E hidden momenta, 
through an analogue model. 

The purely formal analogy between the Maxwell and 
Weyl tensors (their decompositions in electric and mag- 



24 Note that this application is beyond the scope of previous ap- 
proaches to the Papapetrou equation in the context terms of an 
analogy with electromagnetism. As we have seen in Sec. |III| 
the expression Fq = — V(£? • Bq) of the usual linearized the- 
ory (e.g. [2]) is valid only when the gyroscope is at rest (in a 
stationary spacetime). Hence, this expression can only account 
for the coupling between the spin of the gyroscope and the spin 
of the source (that is zero for Schwarzschild black hole). The 
exact equation (5) of [5], derived in the framework of the Quasi- 
Maxwell formalism is valid only for stationary spacetimes, and 
when the gyroscope's worldline is tangent to a time-like Killing 
vector. 



Sec. IV B differs from the physical analogies discussed 
above in the sense that the parallelism drawn is not be- 
tween objects playing analogous dynamical roles (see Sec. 
4.3 of [H] and [33]) in the two theories. Yet it was useful 
in this work as it provided familiar, intuitive tools to solve 
the problem in Sec. |IVB"} From the structure of the scalar 
invariants of F a p for the electromagnetic field of a spin- 
ning charge, we know that in the equatorial plane there 
are velocity fields U a for which the magnetic field B a 
vanishes, thus for which a magnetic the dipole does not 
precess. Given the formal similarity with the quadratic 
invariants of R a p~tS, we anticipated that in the analo- 
gous gravitational setup (Kerr spacetime, which asymp- 
totically, describes the field of a spinning mass) there is 
a velocity field for which H a/ 3 = 0. One must always 
however bear in mind the non-physical character of the 
analogy; as M a p = means not that the gyroscope does 
not precess (which it does, for this velocity field, as shown 
in [15]) but that it feels no force. (Again this exemplifies, 
now in the more obvious context of a non-physical, purely 
formal one, that the existence of an analogy does not im- 
ply any similarity; as the effects at stake are in this case 
actually totally opposite. This observation is crucial for 
the understanding of the physical meaning of the scalar 
invariants of the Riemann tensor, which is discussed in 
the companion paper |15j.) 

In the course of this work a clarification of a number 
of issues concerning the dynamics of spinning particles 
in general relativity was needed. Firstly the problem of 
the covariant, relativistic equations of motion for parti- 
cles with pole-dipole gravitational and electromagnetic 
moments; the gravitational part is well established, but 
difficulties exist in the electromagnetic part, there be- 
ing different versions of the equations in the literature 
and inconsistencies in their physical interpretation. Their 
clarification is the purpose of Appendix |A 1| Then in 
Appendix [B] we dissect the differences between the dy- 
namics of electric and magnetic dipoles (which may be 
cast as the contrast of the two dipole models: a current 
loop in the magnetic case, a pair of opposite monopoles 
in the electric case), which is reflected in the spin vec- 
tor evolution, hidden momentum, proper mass (and its 
variation), and the force exerted on the particle (both 
its time and space projections). In particular, the time 
projections of these forces, their physical content, and 
relationship with the mass of the particle and the work 
done by the fields, is ignored in most literature, or mis- 
understood (e.g. [2U[3ni[3S])- (Failure to notice that the 
force is not spatial, and the proper mass not constant, 
has even led some authors [TB] [TJJ to believe that a co- 
variant treatment in the dipole scheme is not possible.) 
These are rigorously discussed, for the case of the mag- 
netic dipole, in Sees. [Vjand II E[ and in Appendix [B] for 
the electric dipole. In Sec. V A we clarify the meaning of 
the projection of F£ M parallel to U a , as being the rate 
work transferred to the dipole by Faraday's induction, 
and its relation with the variation of proper mass m; the 
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later being further discussed in Sees. HE and VIA In 



particular we clarify that, contrary to usual claims in the 
literature, e.g. [351 123 EE], but in agreement with the 
point of view in |39| . the contribution —ft - B to m is not 
potential energy (as B does no work), but internal energy 
of the particle, which, for a spherical quasi-rigid spinning 
charged body, is shown in Sec. |VI A~2] to be but kinetic en- 
ergy of rotation about the center of mass [cf. Eq. (141)]. 
That is in agreement with the results previously obtained 
in the non-relativistic treatments [64H67] ; as for these, 
our treatment in Sees. |VI A 2| fll E| shows how they fit in 
the relativistic multipole scheme (namely that the kinetic 
energy of rotation is encoded in m) . These are aspects of 
the dynamics of magnetic dipoles (and comparison with 
the gravitational case) that are central in this work; they 
involve however many subtleties which, we believe, can 
only be fully assimilated by contrasting with the case of 
the electric dipole, which we do in Appendix |B) 

Another important clarification regarding the electro- 
magnetic equations was made in Sec. VI A[ where we 
studied Dixon's [29( 138] e qua tions to quadrupole order. 
The dipole torque in Eq. (22) (assuming /2 = crS), only 



causes the magnetic dipole to precess (since f _L S); its 
spin vector would have constant magnitude S = \/S a S a 
accordingly, regardless of the time dependence of the 
field. This seems puzzling, as if one thinks about the 
magnetic dipole as arising from a spinning charged body, 
placed in a time varying magnetic field, a torque find (see 
Fig. |6| due to the induced electric field would be exerted 
on it, changing its angular momentum. This torque is 
known to exist as it has been computed (for some special 
cases) in some non-relativistic treatments [65l467| . The 
constancy of S would also apparently be in contradiction 
with the varying proper mass, and our assessments in 



Sees. V A HE that that variation reflects the work done 



on the dipole by the induced electric field. In Sec. |VI A| 
we clarify these apparent contradictions, and how find fits 
in Dixon's multipole formalism, which indeed involves 
some subtleties: i) the torque ( |122 1 is not manifest to 
dipole order [where only the multipoles of first order in 
r are kept, cf. integral s ([5]) , ([7|] because it involves the 
"quadrupole" moment ( |111| (i.e., the second moment of 
the charge, of order r 2 ); however the work done by it, 
and the corresponding variation of m, are manifest in 
dipole order, as they can be written in terms of fl (that 
is: the variation of angular momentum S caused by T? d is 
not manifest to dipole order, whereas the corresponding 
variation in kinetic energy of rotation is!), ii) The exis- 
tence of this torque is overlooked in the literature [25J [35] 
concerning Dixon's quadrupole order equations (which is 
down to the unclear physical significance of the quan- 
tities involved in the original derivation); in their usual 
form they are equation s for the "canonical" angular mo- 
mentum (S'can)"' 9 , Eq. (A3 1; not for the physical angular 
momentum S a ^ . iii) The confusion between canonical 
and physical angular momentum, and subsequent con- 
flicts from the apparent absence of ~r\ n d , are common also 



in the traditional treatments of classical rotors, see [53] 
for details. This clarification is of special relevance in the 
context of this work as the non-existence of an analogous 
gravitational torque is at the root of fundamental differ- 
ences between the two interactions emphasized herein. 

In Appendix |C 1| we compare the description given 
by different spin conditions, in particular the Dixon- 
Tulczyjew and the Mathisson-Pirani conditions, clear- 
ing up some common misconceptions in the literature 
(Sec. C 2 ). As for the latter condition, on which the exact 



analogies in this work rely, it has been subject of skepti- 
cism due to the exotic helical solutions it admits (in addi- 
tion to more intuitive non- helical one), which have been 
deemed unphysical in some literature [57J IM1 HS1 HZ] ; this 
is dissected in the companion paper [25j , where these as- 
sessments are shown to be unfounded, and this spin con- 
dition to be perfectly valid. And herein (Sec. C 1 ) we 



exemplify, with the applications in this work, that it can 
be the most suitable one. 

As future direction, we plan an investigation of the 
gravito-electromagnetic analogies in the equations of mo- 
tion for spinning particles to quadrupole and higher or- 
ders in the multipole expansion. 



Appendix A: The equations of motion for spinning 
test particles in the literature 

It is perhaps surprising that the problem of the covari- 
ant equations describing the motion of spinning parti- 
cles subject to gravitational and electromagnetic fields is 
still not generally well understood, with different meth- 
ods and derivations leading to different versions of the 
equations, and the relation between them not being clear. 
And that it is the electromagnetic (not the gravitational) 
field that has been posing more problems. In some works, 
e.g. [T8l - f23] . equations of motion been presented which 
are symmetric in terms of electric and magnetic dipoles. 
If not properly interpreted, this would lead to physically 
inconsistent predictions, given the different nature of the 
two dipole models (see discussion below and Appendix 
[B| . Some authors [T51 E] have even concluded that such 
covariant description is not possible. 

The equations of motion for spinning polc-dipolc par- 
ticles are derived in rigorous and unambiguous forms 
in [39], for special relativity, and in [27j in the con- 
text of general relativity; rigorous derivations are also 
in [3T] [5H HH1 |35]; m this case, however, one must be 
aware of the subtleties involving the interpretation of the 
quantities showing up in the equations (namely, the ten- 
sors P a and 5 Q/3 therein). Herein we will explain how 
equations (|8|-([9|) compare to the literature. 

Eqs. (|8|-([9| correspond to equations (6.31)-(6.32) 
of [27] (cf. also (3.1)-(3.2) of [35]), with the following 
simplifications in the definitions of P a and S al3 (= J a ^ 
in 27\): 1) the bitensor — ct ;q , which generalizes the con- 
cept of separation vector of special relativity r a = x a — z a 
(see also e.g. Appendix A.l in [52]), is approximated by 
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r a . Clearly this is not exactly a vector in general rel- 
ativity; but it is a vector to first order (for an explicit 
demonstration, see e.g. Eq. (7) of [55]), whose magnitude 
yields also a first order estimate of the distance between 
two points. And since we are cutting off the expansion 
in the dipole term, which is linear in x a , indeed, to this 
degree of accuracy, a' ,a — —r a . For the same reasons, 
the bitensor of geodesic displacement g K a becomes ~ 6* 
to this degree of accuracy. In this way, Q and |5]) , agree 
with the definitions for P a and S a ^ given in e.g. [H 181 j - 
2) The vector w 1 in the definition of Q a /3 (more precisely 
in the the moments in [5TJ [57J 133), which is a vec- 
tor such that displacement of every point by w 1 d,T maps 
£(t) into E(r + dr) is, to pole-dipole order, to be taken 
as w 1 ~ U 1 . The justification is more easily seen consid- 
ering the case of flat spacetime 25 [2H IH3 , where we have, 
for £ orthogonal to U a : 

^ = ^(l-^^)= £ n ( l + r„a«). (Al) 

This leads to a quadrupole-type contribution to j a/3 [21] 
if a a ^ 0: 

j af} = [ r Q /u> 7 dS 7 

Js(t,U) 

r a fU^dE^ - a a [ r° 'r a fW 'd£ 7 , 

S(t,!7) JT,(s,u) 

which is negligible in this scheme. 3) The moments P a , 
S a P and Q al3 are defined relative to an hypersurface of in- 
tegration E(r, U) which is normal to U a , as done in [21], 
whereas in e.g. [551157] hypersurfaces S(r, P) orthogonal 
to P a are used. That does not change the shape of the 
equations to dipole order though; in the purely gravita- 
tional case (F a P — 0) we know it 26 from the conservation 
equations T a ^ = 0. In the electromagnetic case that can 
be checked comparing with [21] (identifying the appropri- 
ate quantities, see Sec. Al below), or with independent 
derivation in 



[39] . to pole-dipole order the form of the 
equations remains the same. 



1. Dixon's equations 

In later works by Dixon [5TJ [55] [37] EE] these equations 
are presented in a different form, e.g. Eqs. (1.33)-(1.34) 



25 It suffices for this purpose to work in flat spacetime; a gener- 
alization of w a to curved spacetime only amounts to general 
relativistic corrections to something which is already negligible 
in Special Relativity. 

26 In the pole-dipole approximation only terms linear in r are kept; 
since to first order the spacetime can always be taken as flat, 
then by the same flat spacetime arguments (see e.g. _ | 31p one can 
show that T a ^ = implies that the integrals defined at 
a point z(t). over an hypersurface £(t, u), are independent of 



of [35], which are symmetric with respect to the electric 
and magnetic dipoles. Taking into account the different 
signature and conventions, they read: 



1 



1 



D{Poix)a = qF a ' 3 U p + -F^ a Q^--R a ^S^ , (A2) 



dr 



D(S can ) a0 
dr 



2(P D ix) [ "(7 /31 + 



F cc] 



(A3) 



It must be noted that (P Dix ) Q and (S ca „) a/3 (P a , S af} 
in the notation of [5TJ [55] [37] [35]) are not the physical 
momentum and angular momentum given by Eqs. Q 
and ([5| above, but instead contain additional electro- 
magnetic terms, cf. definitions (5.1) and (5.2) of [55]. As 
shown in [S3], {P^) a +qA a = (P can )« and (S can ) Q/3 have 
the interpretation of "canonical momenta" associated to 
the Lagrangian of the system 27 . (P C an) Q is the quan- 
tity conserved in collisions [53], and its time part is the 



scalar conserved in time-constant fields, cf . Eq. ( B8 ) be- 
low (see also Eq. (5.3) of [22], Eq. (29) of [29]). (S can ) afl , 
as explained below, generalizes the canonical angular mo- 
mentum of some non-relativistic treatments [5SH57] . To 
pole-dipole order we may skip the bitensors in these def- 
initions, to obtain: 



(p Di x) a = r + p°, p a 



9 a j*dEp, (A4) 



E(z,C7) 



(Scan)"' 3 = S af} + S' af} , S' af} = 2 f r [a ^ ] pdE^ , 

(A5) 



with 



* a = - [ F a p(z + ur)r du , 
Jo 



(A6) 



¥ 3 = - [ rf^z + ur)r du . (A7) 
Jo 

Note that in the case flat spacetime, these expressions 
are exact J5JJ. The lowest order approximation to these 
integrals is to take only the zeroth order term in the 
expansion of -F^around z, i.e., to take F a @ w constant 
along the particle's size; this is the problem at hand, 
as higher terms in the expansion of F a @ would lead to 
contributions of higher multipole moments to P™ and 
S ,otf) . Thus we obtain: 

P' a = -F a j( P (i) , S' aP = F^ a q^ a (it) , (A8) 

where q a @ is "the charge quadrupole moment", defined 
by (3.9) of [21] . As such S' a/3 is negligible to pole-dipole 



We thank A. Harte for discussions on this point. 
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order, but it is of crucial importance in Sec. |VI| where 
terms up to quad r upol e order are kept. Substituting P' a 



in equations (A2)-(A3) (and taking S af) ~ (Scan)"' 9 ), we 



immediately obtain equations (|Sj)- ([9]) 

We note also that (S can ) al3 is a generalization of the 
concept of canonical angular momentum which arose 
in classical (non-covariant) treatments on electromag- 
netism 65-67 J. In |67j . a canonical angular momen- 
tum vector, Eq. (31) therein, is obtained differentiat- 
ing dC/duj (C = Lagrangian of the system, uj = an- 
gular velocity of the body). Such 3- vector is but a 
non-covariant form for the spatial 28 vector (S can p — 

S^ + S'^, with 



[/ M /2, as can be easily shown. (Scan) 7 = 



5 / 7 



1 

— e 
2 



7 QlaPTTH 



2 l*"!* 



B a 
~2~ 



(A9) 



where we used the orthogonality conditions q a ^U a = 
q a PUfj = 0. If the body has uniform mass and energy 
density, we can write S 7 = (q/2m)B a I^ r , where is 
the moment of inertia defined in Eq. ( 124| ). In this way, 
in the center of mass frame {U l = 0), (Scan) 7 = (0, Scan) 
which reduces to Eq. (31) of [67] . 

The distinction between P^ ix , S™£ n and the ph ysica l 



momenta P a , S a ^ is crucial; overlooking it in Eqs. (A2) 



( A3 ) will then lead one to believe that electric and mag- 



netic dipoles are dynamically similar, since those Eqs. 
are symmetric with respect to them. That, given the dif- 
ferent nature of the two dipole models — the current loop 
(magnetic dipole) and the two opposite charges (electric 
dipole) — leads to physically inconsistent predictions: i) 
the electric dipole would have a hidden momentum B xd, 
just like a magnetic dipole has a hidden momentum fix E; 
this would violate the conservation equations, because in 
the electric dipole model the integral of the electromag- 
netic pointing vector vanishes, ii) a static electric field 
would do no work on a moving electric dipole, nor in a 
rotating one, which we know to be false; iii) the force 
on an electric dipole would have a time projection in the 
CM frame exactly analogous to its magnetic counterpart 
(whose the time projection is the rate of work done on the 
current loop by the induced electric field), which would 
seem to imply that the electric dipole is a current loop 
of magnetic monopoles. iv) also the spatial part of the 
force would not be consistent with the results known from 
previous works on classical electromagnetism. The prop- 
erties of the electric dipoles, with an account of these 



problems and a comparison with the magnetic dipole are 
presented in detail in the next section. 



Appendix B: The electric dipole 

In order to better understand some key issues in this 
work: the physical meaning of the time projection of the 
force on a magnetic dipole, the variation of its proper 
mass, the work done on it by the external fields and their 
relationship with Faraday's law of induction, as well as 
the hidden momentum, it is useful to contrast with the 
case of an electric dipole. 

It is clear from Eqs. (|8])-(|9| that both the force and 
the equation describing the evolution of S a ^ are differ- 
ent for electric and magnetic dipoles. These differences 
arise from the very intrinsic differences between the two 
dipole models: a system of two close monopoles of oppo- 
site charges, which is the case of the electric dipole, and 
the current loop model, which is the case of the magnetic 
dipole. For an electric dipole (fi al3 = 0, q = 0) in flat 
spacetime Eqs. ([8|-(|9| read: 



DP a 
dr 



DS al: 
dr 



E a pdP 



DdP_ 



2P [a U p] 



2d\ a F f) \U 1 , 



(Bl) 



(B2) 



where E a p = F a ~.^U' 1 is the electric tidal tensor [6]. 
Firstly we note that, unlike its magnetic counterpart 
Eq. ((ijla) of Tablc[TJ the force on an electric dipole is not 
generically given by a contraction of a rank 2 tidal tensor 
with the dipole moment vector (only if Dd a /dr = 0). In- 
deed, the force on an electric dipole is not entirely a tidal 
effect: there is the extra term F a pDd@ / dr (overlooked in 
most literature) which does not involve derivatives of the 
electromagnetic field. This term is physically interpreted 
as follows. From Eq. (3.23a) of [3T] we have: 



Dd~< 



J 1 - U^q 



where q is the total charge and Dixon's moment J a , de- 
fined as J a = Jwj; s \ j a w 1 dY>~ f , is the integral of the cur- 
rent density j a in the 3-space orthogonal to U a . Thus, 
for the electric dipole (q = 0), the extra term reduces to 
FZfJ 1 and Eq. (Bl ) may be re- written as: 



F a B jt 



(B3) 



28 Note that definition of (5 can ) 7 is not a dualization of (S l C an) a,3 i as 
neither (S C an) Q/3 or S' a ' 3 are spatial tensors if we take the Pirani 
condition S a f>U p = (not (5 ca n) Q/3 ^ = or {S^)^ Pp = 0, 
as done in 1211 or 1221 . respectively). Hence S la ^ and (Scan) 7 do 
not contain the same information as the corresponding two-forms 
(only their spatial part). 



Hence, the term F^J 1 has a straightforward physical 
interpretation: if the dipole vector d a varies along r (e.g., 
if the dipole rotates) it generates a net electric current 
in the CM frame; therefore, a magnetic force F^J 1 is 
exerted on the dipole, in addition to the tidal term E^d 13 . 
This magnetic force is not a tidal effect, it is simply the 
coupling of F a @ to the current j a . 
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Secondly we note that in the tidal term E a ^dp the 
indices of the tidal tensor are reversed compared to the 
force on the magnetic dipole Eq. |l]la). Since the tidal 
tensors E a p, B a p are spatial in the first, but not in the 
second index, this signals a fundamental difference, which 



is discussed in Sec. B 2 below: by contrast with F£ M — 
B$ a \xp, the tidal term E a ^dp has zero time projection in 
the CM Frame, which means that, as measured in that 
frame , it does no work on the dipole. 



1. No hidden momentum for electric dipole 

Unlike the current loop, a two-monopole type of dipole 
cannot store hidden momentum of electromagnetic ori- 
gin, see e.g. |59j . The expression for the momentum of 



an electric dipole is obtained contracting Eq. ( B2 ) with 
Up, leading to {U a d a = 0) : 



P<* = mU«-— [ -Up = P£ z 



pa 
MiidI 



(B4) 



showing that the only hidden momentum present is the 
pure gauge one -Pj°j dI arising from the spin condition 
(which exists regardless of the electromagnetic multi- 
pole structure of the particle). That is also to be ex- 
pected from conservation arguments, as unlike its mag- 
netic counterpart, the electric dipole does not generate 
electromagnetic momentum (cross momentum P", see 
Appendix]!)]) when placed in an electromagnetic field; the 
integral of the electromagnetic pointing vector vanishes 
[96] . Now consider, for simplicity, a stationary configu- 



ration; in this case the conservation Eqs. (T t( 



o 



imply P" ot — 0; if the dipole was to have any hidden 
momentum, it would not thus be canceled out by the 
electromagnetic field momentum, and violate the conser- 
vation equations. 

This gives a strong hint on how crucial it is to dis- 
tinguish between the physical momentum P a and the 
canonical momentum (Ppix)" = P a + P' a in Eqs. (A2)- 



( A3 ) . These equations are symmetric in terms of electric 



and magnetic dipoles; hence, as can be checked contract- 
ing (A3) with U a , or directly from (|A8[), (i-bix) Q includes 



a term e a 6fia d B^W 7 analogous to the hidden momentum 



hidEM 



9\- 



[i E U T of the magnetic dipole: 



p,a = pf* UaU a _ e a f)fi<7 d e B fJ -U' T 



2. Proper mass and time projection of the force in 
CM. frame 



We have seen in Sec. [VA] that the force on a magnetic 
dipole is not a spatial vector; it has a time projection in 
the particle's proper (CM) frame, which we showed to be 
the power transferred to it by Faraday's induction; and 
which is manifest in a variation of the dipole's proper 
mass, cf. also Sec. |IIE| It is thus important to contrast 
with the case of the electric dipole, where such induction 
phenomena cannot occur, due to the different nature of 
the dipole mod el. 

Contracting (Bl) with U a one obtains: 



Fc\U a — —E. 



Dd* 1 
^r~ 



(B5) 



where E a = F a/3 U/3 is the electric field as measured 
by the test particle. Hence, also in this case the time 
projection is generically non-vanishing. However, it is 
very different from the time projection of the force on 
a magnetic dipole, Eq. ( 96 1 ; the former has nothing to 



do with electromagnetic induction, as it arises from the 
variation of the dipole vector d a along the test particle's 
worldline, and not from the variation of the field as in 
the magnetic case. Indeed, since the electric dipole has 
no hidden momentum of electromagnetic origin and thus 
P a a a = =>■ dm/dr = — U a DP a /dr, we also have 



dm 
~dr~ 



~FS.Ua — E. 



Dd* 1 



(B6) 



(compare with Eq. (34)). In particular, we note that 



assuming the dipole vectors d a to be constant along the 
particle's worldline, we have F°yU a — dm/dr = 0, which 
contrasts with the situation for a magnetic dipole, where 
dm/dr is zero only if DB a /dr = 0, and not D[i a /dr = 0. 

The result ( |B5| ) makes sense: J 1 is the total current as 
measured in the dipole 's frame. When it is non- vanishing 
(for instance, due to a rotation of the dipole) a non- 
vanishing work, in this frame, is done on the dipole by 
the electric field. Note that the density of Lorentz force 
on a arbitrary current is f a = F a ^jp (see e.g. [70]); if 
the current originates from a moving charge, the force 
is spatial with respect to the charge's velocity U a , since 
j a = pU a (i.e., U a is parallel to the current) and thus 
f a U a = 0. But that is not the case when the current 
arises from rotation or oscillation of a dipole; in this case 
j a is not parallel to U a (and the 3-current j is non-zero 
in the test particle's frame) and f a is not spatial with 
respect to U a . 



Thus, if one overlooks the fact that (Pdix)" is n °t the 
physical momentum, one is led to believe that the elec- 
tric dipole has a hidden momentum just like a magnetic 
dipole; which not only would make no sense for a two- 
monopole model of dipole, as it would violate the conser- 
vation equations. 



3. Time component as measured by static 
observers and work done by the electric field 

As discussed in Sec. |Vj the projection of the force 
along a congruence of time-like vectors u a yields the rate 
of work done on the test particle as measured by the 
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observers 0(u) of 4-velocity u a . And we have seen in 
Sec. VB| for the case of a magnetic dipole, that if a con- 
gruence u a exists along which F a p is covariantly co nstan t 
F a ? u 1 = (which we dub static observers, cf. Sec. VB I, 



the projection of the force along it — F^ M u a (i.e., the time 
component as measured in the u l — frame) vanishes, 
meaning that no work is done on the dipole. And that 
was understood through an exact cancellation between 
the variation of translational kinetic energy of its center 
of mass, and the variation of its internal energy (the lat- 
ter being, from the point of view of the CM frame, the 
work done on the dipole by Faraday's induction, cf. Sec 
VA}. 



In the case of the electric dipole, since it is not a cur- 
rent loop, the situation is different, as expected. The 
projection of the force along a congruence u a for which 
F a ^ui = is not zero, by contrast with its magnetic 
counterpart. Let again U a be the 4-velocity of the test 
particle, related to u a by Eq. (90); we have: 



V 



-F> Q = ^E^^aVv^ + (E u ) a 



Dd a 



(B7) 



where (E u ) a = F al3 u fj and (E u ) Pl = Fp^y? are, re- 
spectively, the electric field and electric tidal tensor mea- 
sured by the observers 0(u). 

Thus generically —F^u a ^ 0. That is the field does 
work on the dipole, as measured by 0(u). The first term 
is a quite natural result: in a inhomogeneous electric field 
((E u ) a /3 0), a net force is in general exerted on a elec- 
tric dipole; if it is allowed to move (v a ^ 0) that force 
does work, which (if no other forces are present) trans- 
lates into a variation of translational kinetic energy. The 
second term contributes when Dd a /dr ^ 0, and is non- 
zero even if the fields are uniform. Again this result is 
familiar: it is the work done by the electric field when 
the dipole rotates or oscillates. 



4. Conserved quantities, proper mass and work 
done by the fields 

In order to better elucidate the relationship between 
the work done by the fields and the variation of proper 
mass, as well as the role of electromagnetic induction, 
we will compare, in a static electromagnetic field, three 
test particles: a point monopole charge, an electric dipole 
and a magnetic dipole. Let u a be the 4-velocity of the 
inertial frame relative to which the fields are static. Then 
u a preserves the electromagnetic field: L u F a ^ — 0, and, 
therefore, from the constancy of expressions (5.3) of [22] . 
or (29) of [29 , we have that: 

P£ ix u a + qA a u a = P a u a + (E u ) a d a — q<f> = constant 

(B8) 

where <fi = —A a u a denotes the electric potential mea- 



sured by the static observers u a . Using Eqs. (26)- (27 1, 
it is useful to re- write ( B8 1 as 



T + V + E] 



hid 



constant 



(B9) 



where V — ~(E u ) a d a + q<j) is the potential energy of 
the particle under the action of the field, T = (7 — l)m 
is the kinetic energy associated to the translation of the 
center of mass, and i?hid = ~Phid u a the "hidden energy" 
(i.e., the time component of the hidden momentum in 
the chosen frame) 



VB 



see Sec 

is the particle's 4-velocity, cf. Eq. (|90) 



7 



^-U a u a , and U°> 
In this section, 



for simplicity, we will ignore the inertial hidden momen- 
tum fjfjjT, as it is of little relevance compared to the 
other terms, and the motion effects induced by it pure 
gauge and confined to the disk of centroids, see [UJ HO] 
(a region within the convex hull of the body's worldtube, 
typically much smaller than it). Thus, we take herein 



P 



P 



hid * hidEM- 

Point monopole charge (d a 
case, condition (B9) means m+T+q4> 



r>* 
Miid 



= 0). — In this 
constant. There 
is no exchange of energy with the proper mass of the 
particle, which is a constant: 



dm 
~dr~ 



DP* 
dr 



-qF a0 U a U P =0. 



This only tells us that, in a time independent electro- 
magnetic field the familiar notion of "total mechanical 
energy" of the particle — kinetic energy T, plus the po- 
tential energy V = q<p — is a constant of the motion. 
Therefore, in the absence of other external forces, every 
gain in kinetic energy must come from the electric po- 
tential energy V, hence there is no doubt that the elec- 
tric field E(u) — ~V(f> is doing work. The rate of work 
V is given by the time projection of the Lorentz force 



Fg = qF^Up: 



V 



-F£u c 



qj(E u ) 



dV_ 
~~dr~ 



Electric dipole (q = P^ d = 0). 
means in this case m + T — (E u ) a d a 



Condition (B9j 
= constant; from 
Eq. (|B6|) the proper mass m is not constant; this means 



that energy is exchanged between the three forms: poten- 
tial energy V = —(E u ) a d a , translational kinetic energy 
T, and internal energy m. Two special sub-cases are par- 
ticularly enlightening: 

1. the case where the dipole vector is covariantly con- 
stant Dd a /dT = 0, implying dm/dr = 0; in this 
case the energy exchange is similar to the monopole 
charge: every gain in translational kinetic energy 
comes from the potential energy V. It is clear that 
the electric tidal field is doing work, at a rate (cf. 
Eq. dB7l): 



V 



7 (£ U W 7 ^ 



dV 



2. The case where the dipole's CM is at rest (v = 0, 
r = t), i.e., T = 0; in this case, m — E-d = constant 
(note that E(u) = E if v — 0), and the energy 
exchange occurs between the potential energy V — 
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—E-d and proper mass to (which includes rotational 
kinetic energy of the particle). A familiar case is 
the quasi-rigid dipole (d a d a = constant) which is 
allowed to rotate (Dd/dr = ilx d ). In this case, cf. 
Eq. (B7|, the electric field does work at a rate: 



V 



-F>o 



(nxd) -e 



dm 
~dt 



dV 
~dt 



This is simply the rate of work V = f ■ ft done 
by the torque f — d x E which is exerted on the 
dipole by virtue of Eq. ( |B2[ ); that can be seen us- 
ing the vector identity (Q X d) ■ E = (d x E) ■ Q. 
The torque f causes the dipole to rotate and gain 
kinetic energy at a rate V . Again we see how i mpor - 



tant it is not to confuse -Pg ix , S"f a in Eqs. (A2) 



( A3 ) with the physical momenta. Overlooking the 



distinction would lead to the conclusion that, just 
like for a magnetic dipole, the electromagnetic field 
does no work on a moving electric dipole, nor in 
a rotating one — that would be true for a mag- 
netic dipole only, where, for fl = aS, the torque 
t = DS/dr = fix B, cf. Eq. (24), is orthogonal 
to S, and therefore only causes the dipole to pre- 
cess, implying no change in S 2 or rotational kinetic 
energy. 



Magnetic dipole (q 
in this case to + T - 



-- d a = 



0). 



Condition (B9) means 
constant; if we take fl = aS, 
from Eq. (36) we have to = too — ^ a Bp, and thus the 
condition becomes T — fJ. a Bp + -Ehid = constant. The 
energy exchange is between translational kinetic energy, 
proper mass and hidden energy. There is no potential 
energy involved (cf. [6"1H6"T]), which is consistent with the 
fact that the static field does no work on the magnetic 
dipole: 



V 



, 



cf. Eq. (100). A case of interest in the context of this 
work is the one depicted in Fig. ^jp, a magnetic dipole 
falling into a magnet along the field's axis of symmetry. 
In this case P£ id — E^id = 0; implying T — [i a Bp — 
constant. The energy exchange is only between trans- 
lational kinetic energy and proper mass; every gain in 
the former must come at the expense of latter (which, 
if the magnetic dipole is a spinning particle, is but rota- 
tional kinetic energy, cf. Sec. 
that there is a force driving t 



VI A 2 and 



). Note 

he dipole into motion; but 
what the field does is to interconvert translational kinetic 
energy into rotational kinetic energy |65[ 166] (or other 
forms of internal energy, reflected in the proper mass of 
the dipole). The variation of proper mass is the one that, 
when analyzed from the CM frame of a dipole moving in 
an inhomogeneous field, is the work transferred to the 
dipole by the electric field induced by the varying (along 
the test particle's worldline) magnetic field, see Fig. [3^,. 



Appendix C: Comparison of the Mathisson-Pirani 
and the Tulczyjew-Dixon spin conditions 



As is clear from the discussion in Sec. |II A[ the two con- 
ditions arc different, but equivalent, ways of describing 
the motion of a spinning body; they differ in the choice 
of a representative point in the body; hence we can say 
they are a matter of gauge choice. A supplementary con- 
dition is needed to close the equations of motion ([8|-(|9|; 
the condition S a "up = for some unit time-like vector 
u a ensures that we are taking as reference worldline the 
center of mass as measured by some observer of 4- velocity 
u a . Whether u a is itself the 4- velocity of the center of 
mass (Pirani condition), or it is parallel to P a (Dixon 
condition), or it corresponds to the static observers in 
Schwarzschild spacetime (Papapetrou-Corinaldesi condi- 
tion), or any other type of observer (there is an infinite 
number of possibilities), the option should be based on 
convenience. 

For a free particle in flat spacetime, clearly the Dixon 
condition provides the simplest description for the center 
of mass motion, which is straightlinc motion, as a non- 
degenerate solution. The Pirani condition includes in this 
case Dixon's solution as a special case, but also the helical 
solutions which are more complicated descriptions. 

However, in the presence of gravitational and electro- 
magnetic fields, Dixon's solution, in general, no longer 
coincides with any of Pirani's; and it turns out that the 
latter may actually provide the simplest and more in- 
tuitive description. This is actually always always the 
case for the equation of spin evolution: in the absence 
of electromagnetic field (or other external torques) S a 
is Fermi-Walker transported if S a ^Up = holds; i.e., 
the gyroscope's axis is fixed relative to a Fermi-Walker 
transported frame. This is the most natural description: 
in a gravitational field, and in the absence of other ex- 
ternal fields, no torque is exerted in a test particle with 
only pole-dipole moments (no quadrupole or higher mo- 
ments), as Eq. (22) shows; such particle is thus an ideal 
torque-free gyroscope. It is thus expected (being a gy- 
roscope an object that opposes to changes in direction) 
to have a rotation axis fixed with respect to the Fermi- 
Walker transported tetrad, since the latter is by defi- 
nition locally non-rotating. Actually the Fermi-Walker 
transported tetrad gets its physical significance as a lo- 
cally non-rotating frame precisely by being anchored to 
local guiding gyroscopes [3T| (defining the "local compass 
of inertia" [5]). Only if Mathisson-Pirani condition holds, 
is a gyroscope Fermi- Walker transported. Dixon condi- 
tion yields a different equation, Eq. (C16) below, mean- 
ing that S a undergoes transport orthogonal to P a ; there 
is no conflict because the transport is along a different 
worldline, but the description is not as sound. 

The Mathisson-Pirani condition may also provide the 
simplest description of the force / center of mass motion, 
and examples of that are the are the applications in Sees. 



IV B and VB which we shall discuss next. 
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1. Comparison of the spin conditions in the 
applications of Sec. Ill 

We start with the problem of fall along the axis in a 
Kerr spacetime discussed in Sec. |VB| 

Firstly, in order to compare the different spin condi- 
tions for the same problem, we must first establish how 
we ensure that we arc dealing with the same particle. A 
pole dipole particle is characterized by its two moments: 
P a and S a/3 ; these are defined, cf. Eqs. @-f5]), with re- 
spect to a point z a (r) and an hypersurface of integration 
£(it, r) (which may be thought as the instantaneous rest 
space of the observer measuring them); hence we pre- 
scribe the particle by giving these moments. Different 
representations of the same particle must yield the same 
moments with respect to the same point and observer. 
This is a crucial point made in |25j . 

We intend to describe what we naively expect, from 
symmetry arguments, to be a the fall along the z axis 
(8 = in Boyer-Lindquist coordinates) of a gyroscope 
whose initial position is some point along the z axis, 
and its initial velocity U and spin S are also along the 
axis of symmetry. It turns out, however, that such naive 
prescription of initial conditions not only does not com- 
pletely determine the problem, nor does it ensure its axial 
symmetry. The reason is that in Relativity, and in con- 
trast with Newtonian mechanics, the motion of the CM 
of an extended body is not fully determined by the force 
laws given its initial position and velocity; one needs a 
third condition, which is to prescribe the field of unit 
time-like vectors u a relative to which the CM is com- 
puted (i.e., the field such that S a Pup = 0), and which 
for an arbitrary choice breaks the axial symmetry. The 
momentum-velocity relation also depends on this choice, 
cf. Eq. (30), implying (hidden momentum) that U will 
not in general be parallel to P, so that they do not both 
lie along the z axis. Note that, as explained in Sec. |II D 



the acceleration of the CM does not originate solely from 
the force, but also from the variation of field u a along 
the CM worldline. 

In order to prescribe an axisymmctric problem, we ten- 
tatively set up as initial conditions Ui n = Ue z (by e z we 
denote, in Boyer-Lindquist coordinates, the basis vector 
e r = d/dr along z), as initial position of the CM Xq M (u) 
some point in the z axis, and also some vector field u a 
such that u = ue z . In this case we expect the solution to 
be an axial fall. 

The Pirani condition S a0 Up = 0, i.e., u a = U a with 
this setup, initially satisfies the symmetry requirements, 
so let us start with it. In this case the momentum be- 



comes (cf. Eq. (25)) 

P a = mU a 



MS 



(CI) 



uation of motion F, 



and the spatial part of the ec 
DP> 3 /dr = -W^Sa (cf. Eq. §lb)) follows as 







Initially, with Uj n = U°e + U z e z , one obtains Fein = 
— H za S a e z (it is straightforward check that along the axis 
we have W a = if i ^ z); thus the force is along z as 
expected, which could have been inferred from symmetry 
arguments, given the axial symmetry of the initial setup 
and the fact that H. a p = R a ^vU li \J 1 ' depends only on 
U a . It is clear from the equation above that the most 
natural result, which is motion along the z axis, with 
the body accelerating in the same direction of the force 
(and of S): a = ae z S Xjj a = 0, implying P a — 
mU a , and Fq = ma a , is one 29 of the possible solutions of 
(C2). It is a "non- helical" solution since it is straightlinc, 
thus the description we seek. Hence we have solved our 
problem, the momentum- velocity relation relation having 
been naturally established in the course of the analysis. 

Now let us compare with the equivalent descriptions 
for this problem given by other spin conditions. For a 
generic field u a with u not lying along the z axis, we no 
longer have axial symmetry and therefore we should not 
expect an axial fall as a solution. We know however from 
the analysis above with Pirani condition, that the axial 
fall indeed occurs given the appropriate initial conditions; 
so what we can hope for in the framework of a different 
spin condition, in general, is it to give a different (possi- 
bly exotic) but equivalent description of the same motion. 
However we would not know how to set it up, i.e., how to 
prescribe its initial conditions, without using the infor- 
mation extracted from the analysis above (it would not 
be clear that such motion was possible in the first place) . 
Note that if one naively sets up Xq M {u) lying on the z 
axis, and then P or U (there is also an ambiguity on this 
choice, as they are not parallel, cf. Eq. (30)) along z, the 



ma + 



D(Sxua) 



Fr 



-W a S a ei 



(C2) 



solution of the equation of motion in general not will not 
be an axial fall, as it will not even be a solution equivalent 
to it up to a gauge transformation. We ensure that we are 
dealing with the same spinning particle, by imposing, as 
discussed above, P a , and also S a ^ about the same point 
Xq M (P), to be the same. Hence P — Pe z . In general, 
the condition S a °up = 0, for some family of unit time- 
like vectors u a defined along the reference worldline, will 
define a centroid Xq M (u) at a different point compared 
to Xq M (U), not lying on the z axis, again manifesting the 
fact that the problem is no longer axisymmetric. Since 
U a will not in general be parallel to P a , cf . Eq. ( 30 1 , the 
centroid also does not move parallel to the axis. Writing 
S a p = £a/3/i^('S'")' i w ly , where (S u ) a denotes the new spin 
vector, the force now reads 



F§ = -^uX^S^ = -*B" x U»u»(S u ) X , (C3) 



29 Other solutions were possible, this is due to the incompleteness 
of the gauge fixing of the Pirani condition. Note however that 
these are not equivalent solutions (not just a matter of choice 
between helical and non- helical representations), because, since 
J/j" is fixed, they correspond to different P^'s. They correspond 
to helical representations but of different solutions. 
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which depends both on U a and u a , and in general will 
not be also along the z axis . This clearly leads to a more 
complicated description of the same problem. 

The case of the Dixon condition exemplifies some of 
these difficulties. Firstly we face the complicated equa- 
tion relating P a and U a [231157] 



U a 



M 2 \ 4M 2 + R^sS^S^ J { ' 



which in general are not parallel; so if one wishes to study 
the problem of fall along z, we face the difficulties setting 
up initial conditions discussed above. For instance the 
ambiguity of which vector, U or P (if any) should be ini- 
tially chosen parallel to z. Without using the knowledge 
extracted above with Pirani condition, and with Eq. ( C4 ) 



coupled to ( C3 ) it would be difficult to conclude that the 



axial fall occurs, and to what kind of solution it corre- 
sponds in this gauge — namely if it is or not an axisym- 
metric problem in this gauge. In order to represent the 
same spinning particle as above (which we know to repre- 
sent the desired problem) it is P that we must set along z, 
as discussed above. To find the force, given by Eq. ( C3 ) 
above with u a = P a /M, one needs to know both 30 U a 
and P a (not only U a , as with Pirani condition). It can 
then eventually be shown from (C4|, as done in e.g. [55] . 



even though not that easily, that for this setup, using 
Dixon condition we have P a = mll a , and therefore the 
solution coincides with the one obtained above with Pi- 
rani condition. Hence we end up precisely with the same 
solution, but taking a more complicated path. 

In Sec. |IVB"3] we concluded that in the equatorial plane 
of Kerr-dS, for suitable r and ir, spinning particles move 
in prograde circular geodesies. And we were able to do 
it only because we used Pirani condition. When this 
condition holds, the force is given by a contraction of the 
gravitomagnetic tidal tensor H Q/ 3 with S a , cf. Eq. (JTJ 1 b); 
and from the analysis of the curvature scalar invariants 
we deduced that in the equatorial plane there is a velocity 
field Eq. (78 1 for which H Q( g = 0; and that at some r = 
r gC o, solution of Eq. (89), that velocity corresponds to the 
velocity of the circular geodesic of a monopole particle. 
Along that circle, the equation of motion reduces to 



DP a 
dr 



= <^ ma a - e' 



U' 



D(S'V) 



= 



(C5) 



admitting a a = as trivial solution, which is obviously a 
"non-helical" one; in that case the spinning particle will 
move along the circular geodesic. We would not be able 
to reach this conclusion using other spin conditions: for 



30 It is for this reason that, although we can write, from Eq. j30| , 
the expression P a = M 2 U a /m + S al3 (F G ) l 3/m analogous to 
| |C1| , we cannot apply the same reasoning as with Pirani condi- 
tion, because we only know in which direction the force points 
after solving 1C41. 



a different condition S a ^up = 0, ^ U a the force is no 
longer given by a contraction of the spin vector with the 
magnetic part of the Riemann tensor H Q ^ (but instead 
by a rank two tensor TA. a p = -kR afJ- p u U ll u 1 ' involving both 
vP and UP , cf. Eq. (C3)), and therefore a similar anal- 
ysis in terms of scalar invariants of the Riemann tensor 
is not possible. In particular, starting with Dixon condi- 
tion S a PPp = 0, it would be virtually impossible, dealing 
the system formed by Eqs. (C9) an d ([9]), coupled with 
the momentum- velocity relation ( |C4[ ) , to ever notice this 
effect. 

As for the application in Sec. IV A the motion in 
Schwarzschild spacetime of a particle with radial ini- 
tial velocity, firstly notice that, contrary to what one 
might naively expect, we cannot extract from symme- 
try arguments much information about the force on a 
particle with generic spin S a , as, due to it, the prob- 
lem does not have spherical symmetry. Indeed, a force 
orthogonal to e r arises in the analogous electromag- 
netic setup, cf. Eq. ( |50| . Using Pirani condition, with 
U = U a eQ + U r e r , we have, as discussed in Sec. IV A 






i aj 3 = =>■ DP a /dr = 0. Hence again we have (|C5|) as 
equation of motion, with trivial solution a a 



pa 



mil", i.e., the gyroscope radially moves along a geod esic . 
In the case of Dixon condition again we face Eq. (C4) 
relating P a with U a , not being clear which one to set 
initially along e r ; using the information extracted with 
Pirani condition, in order to repres ent t he same prob- 
lem, it must be P; then using Eq. (C14), and noticing 
that (H p ) tt/3 = 0, we see that, as expected P a \\ U a and 
DP a /dr = 0. 

Now let us check the case of the electromagnetic anal- 
ogous problem: magnetic dipole with initial radial veloc- 
ity in the Coulomb field. Firstly we note that due to the 
(gauge independent) hidden momentum P^dEM! m g en ~ 
eral P a cannot be parallel to U a ; furthermore, F£ M ^ 0, 
and a a ^ 0; therefore, for the case of the Pirani condition, 
for a generic S, and in the exact case, it is not trivial to 
prescribe the initial conditions for the non-helical motion 
we seek (although in this particular example, to first or- 
der in S, we can impose it by simply taking S a ^ap ps 0). 
And in the case of the Dixon condition we face again a 
complicated equation relating P a with U a (and therefore 
Fg M with a a ), which is Eq. (35) of [29]. Hence, in the case 
of this application (not to compare with the gravitational 
case, but solely to to solve the equations of motion), a 
condition analogous to the Corinaldesi-Papapetrou con- 
dition, S a Pv,f3 — 0, with u a corresponding to the static 
observers: u a = 6$ = d/dt, seems the most appro- 
priate: it makes the inertial hidden momentum vanish: 
S aP Du p /dT = 0, leading to P a = mU a + P£ dEM , which 
is the simplest momentum-velocity relation possible for 
this problem. 

The bottom line is that dwelling on the relation be- 
tween U a and P a is sometimes unnecessary work; in the 
absence of electromagnetic field, the difference between 
P a and mll a is pure gauge, as shown in Sec. |II D) The in- 
ertial hidden momentum vanishes with a suitable choice 
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of the field u a relative to which the CM is evaluated 
(i.e., choosing it parallel transported along the reference 
worldline, see Eq. (30)), leading to P a \\ U a ; this may 
actually be chosen as a spin supplementary condition, for 
a detailed discussion on it we refer to 40 . The spin con- 
dition is gauge freedom, and as such one should choose, 
in each application, the one that suits it most. 

In the case of the applications presented, on top of 
the arguments above, we add the most important one in 
the context of this work, which is the fact that it is the 
Pirani condition that allows for the exact gravitoelec- 
tromagnetic analogies, and that it was thanks to that 



analogy that 1) in the application of Sec. VB we under- 
stood the time component of the force, its fundamental 
difference with its electromagnetic counterpart, and were 
able to explain the spin-spin interaction energy with the 
absence of gravitational electromagnetic-like induction / 
the non-zero work done by the gravitomagnetic field; 2) 
in the case of the applications in Sees. |IV B| and |IV A| 
we anticipated the existence of the circular geodesies for 
gyroscopes in the first case, and of the radial geodesies 
in the second case, which are related with very intrinsic 
differences between the two interactions that are made 
transparent in the symmetries of the tidal tensors (a for- 
malism that relies on the Pirani condition). We believe 
this is enough to convince the reader that indeed, for 
these applications, the Pirani condition represents a clear 
advantage. 



2. The difference between the two conditions for 
classical particles 

We have seen in [35] that for a free spinning particle 
in flat spacetime the worldlines defined by Dixon condi- 
tion (corresponding in this case to Mathisson non-helical 
solution) and all the Mathisson helices remain within a 
worldtube whose radius corresponds to the minimum size 
[4Tj [49] a classical spinning particle can have (if it is to 
have finite spin without violating special relativity). Let 
S+P denote the angular momentum taken about the cen- 
ter of (P) as measured in the P % = frame, i.e., 
SfPp 



0; and let S" be the corresponding spin vec- 
tor: = e a ^ v S^P u /M. The radius of such worldtube 
is [25J SOI E5 -Rmax = S*/M. However, in the presence 
of strongly inhomogeneous external fields, the point we 
chose to represent the particle (i.e., where the fields are 
evaluated) makes a difference; moreover the force equa- 
tion itself changes as we have seen in the previous section, 
cf. Eq. (C3). It is shown in 00], for the gravitational case 



(Kerr spacetime) , that the solutions corresponding to dif- 
ferent spin conditions diverge. Hence, depending on the 
supplementary conditions, we seemingly get different re- 
sults. What this means however (this in agreement with 
the point of view in [10]) is the breakdown of the pole- 
dipole approximations (not that the spin conditions are 
not gauge after all!). The pole-dipole approximation is 
only acceptable when the choice of center of mass (and 



the spin condition) does not matter; i.e., when the scale 
of variation of the external field is much larger than S/M, 
and if we follow the particle for a period of time that is 
short enough so that the many solutions do not diverge 
outside the worldtube. 

Finally, there arc a number of misconceptions which 
are usual in the literature regarding the effective differ- 
ences between the two conditions, which we shall clar- 
ify. It is often asserted (in the absence of an electro- 
magnetic field), that the two conditions are equivalent 
if terms quadratic in the spin are neglected [OH [HU [2H] , 
on the basis that, according to Eq. (C4), the difference 
U a and Pp/M is of order 0(S 2 ), which it is claimed 
to lead a difference between the mass dipole moments 



corresponding to each of the conditions, Eq. (13 1, of 
S ap Pp/M - S ap Up ~ 0(S 3 ). This is not correct; the 
analysis in [25], for the special case of flat spacetime, 
makes clear that the maximum distance between the 
worldline specified by S+^Pp = and the helical solu- 
tions allowed by S a ^U = is S*/M (i.e., \(S af} P /M - 
S a P Up)\ m ax = <?*)■ Also, the difference between mU a 
and P a is the hidden momentum; which, in the CM 
frame, reads Phid = S x a with Mathisson-Pirani con- 
dition, cf. Eq. (29). Thus in general it is not of order 
0(S 2 ) as stated m|81 ; and not even necessarily of order 
0(S 1 ): for an helical solution, since \a\ = us 2 R = Mv/S 
(see [25]), we have (in the CM frame) |P - mU\ = Mv 
which is of zeroth order in S. The origin of these mis- 
understandings is the fact that the momentum velocity 
relation depends on the spin condition, as is clear from 
Eq. ( |30[ ) ; and in particular Eq. ( C4 ) , on which the analy- 
sis in [69 , 81, 99 are based, already assumes the condition 
S^Pp = 0. 

We note also that it is not only when one compares an 
helical solution of S al3 Up — and the worldline fixed by 
S a "P/3 = that the difference disagrees with the one usu- 
ally mentioned in the literature; it also occurs between 
the latter and Mathisson-Pirani non-helical solution, if 
external fields are present. Let us see this with a con- 
crete example. Consider first a magnetic dipole at rest, 
in the absence of any external fields. In this case, the non- 
helical solution allowed by Pirani's condition S a ^Up = 
agrees with Dixon's center of mass worldline defined by 
S a PPp = 0, since in this case P Q = MU a , as shown 



by Eq. (C4). Thus, in this special case, both conditions 
pick the same reference worldline, i.e., a center of mass 
as determined in its own rest frame (a proper center of 
mass). Now bring a spherical charge close to the dipole, 
as depicted in Fig. [7j According to Pirani's condition, 
nothing changes, the center of (U) is still com- 

puted in the same frame. But it is not so with con- 
dition S al3 Pp — 0, due to the hidden momentum term 
-PhidEM = (1 XuE. This condition means that the center 
of (P) is computed in the zero 3-momentum 

frame (P = 0), which is a frame moving, relative to 
the rest frame of xg M (P) (the U = frame), in the di- 
rection along PhidEM- In order to see this, let us first 
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take the perspective of the P = frame; the vanish- rate to quadrupole order, embodying a quadrupole hid- 



ing of P in this frame implies, via Eqs. (26) and (28), 



U = — PhidEM/w; i.e., that Xq M (P) moves in this frame 
with 3-velocity v = U/~/, 7 = -U a P a /M. Therefore, 
the P — frame moves with velocity — v relative to the 
center of mass x^ M (P). This implies (see Fig. 1 of [2"5] 1 
that Xq M (P) is shifted to the left relative to Xq M (U) 
(which is where x^. M (P) was in the absence of the elec- 
tric field). Note that, since v is constant, it follows from 
Eq. (27) of [25] that the two centers of mass Xq M (U) 
and Xq M (P) are points at rest with respect to each other 
(and to the charge in the apparatus of Fig. [7j i.e., they 
are at rest with respect to the static observers therein, 
according to the definition given in point [8] of Sec. ID) 



The shift Ax a = x% M (P) - xg M (C7) is given by Eq7~[7) 



of [25] (setting therein u a = P a /M; and Eq. (8) therein 
yields — Ax a ). Thus, considering fi = aS, it is of order 
0(S 2 ). 



a. The 0(S 2 ) terms 

We will now clarify the above discussion around the 
0(S 2 ) terms, and the reason why this order of magnitude 
is relevant. It is often asserted in the literature, in diverse 
contexts, that the 0(S 2 ) contributions to P a are negli- 
gible, so that in the purely gravitational case one may 
take P a ~ mU a (which would imply, as explained above, 
that Mathisson non-helical solution would be identical 
to the worldline specified by S a/3 Pp = 0). Some authors 
[100) argue this is true in Post-Newtonian approximation 
(which is supported by our analysis in Sec. Ill); others 
argue it is so for "physical acceptable solutions" [69]; or 
that 0(S 2 ) terms are of the order of magnitude of the 
quadrupole terms, which are to be neglected in the pole- 
dipole approximation [26] , But usually without a full, 
satisfactory justification (to our knowledge, such justifi- 
cation has not been put forth) ; that is our purpose in this 
section. 

We will start by the gravitational case. We want to 
clarify if the second term in the expression for the mo- 



mentum ( 25 ) is negligible in a pole-dipole approximation. 



We have already seen that for a generic spin condition 
(that is, for an arbitrary choice of centroid) it is not; 
an example are, as discussed above, Mathisson's helical 
motions where this term is not negligible (nor of order 
0(S 2 )). In this analysis we will use the momentum ve- 
locity relation given by the Mathisson-Pirani condition, 
Eqs. ((26|-(27), and consider the non-helical solution (in 



the purely gravitational case it is close to the worldline 
specified by S al3 Pp and the following analysis is similar 
for both). To dipole order, the only form of hidden mo- 
mentum present in a purely gravitation al sy stem is the 
inertial one P hidI = e^pS^a 1 (see Sec. |lID[). However, 
quadrupole (and higher order) terms contribute also to 
the hidden momentum. Contracting Eq. (142) with U\, 
yields an expression for the momentum, Eq. (26), accu- 



den momentum term P^q'- 



pa _ a qP n 1 1 pa 
Mud — e 7^° a + MiidQi 



pS — 7d^[J n»] 77 
MiidQ — J " ' ' 



7a/3 



It is worth noting that P^q, unlike the first term P" idI , 
is not gauge (in general it cannot be made to vanish by a 
change of reference worldline) . Pj^ id Q can be decomposed 
using Eqs. (15b , leading to terms like (mo)' yS (E u ) ysU a 
where (toq)'*^ is the mass quadrupole defined in (149). 



Note also this important aspect: only the traceless part 
of (itlq) iS contributes in vacuum, since in that case 
(E u ) a a = (it is show using the above mentioned de- 
composition of J a p-yS that the contribution of (mg) 7 7 to 

P hidQ is zer0 )- 

We will exemplify the orders of magnitude of these 



hidden momentum terms using the application in Sec. Ill 



For a non-helical solution, the acceleration comes in first 
approximation from the force Fq, i.e., the terms kept 
in (47), which are of order O(S). We will focus on the 



second term of (47), since it is dominant for large r (the 



analysis of the first term is analogous); the contribution 
(i.e. S Xjy a) to the momentum is then typically: 



pa 
MiidI 



\S xu a\ 



-\(E U ) 



vm\(E u ) i:j \R 2 



(C6) 



where we used 5 ~ mt> rot P, and allowing u rot (the rota- 
tion velocity of a point at the body's surface) to take its 
maximum value v rot = 1; R is the radius of the spinning 
test body. Now, if no assumptions are made about the 
body's multipole structure (such as, for instance, a close 
to "spherical" shape), the traceless part of (toq) 7<5 (which, 
as explained above, is the only part that contributes) 
must be allowed to be of the order (itlq) 15 ~ mR 2 . That 
means that P^q involves terms of the same order of 
magnitude of P£? 



di- 



\(m Q r d (E u ) 7S v l \ ~ vm\(E l %\R 2 ~ \S x v a\ . 

This justifies the neglect of 0(S 2 ) hidden momentum 
terms in the pole-dipole scheme, and the claim in [26j . 
if no assumptions are made about the body's multipole 
structure (i.e., the only assumption made about the body 
is that it be small enough so that higher order multi- 
poles are negligible). On the other hand, if one assumes 
that the body is nearly "spherical", so that its traceless 
quadrupole is very small, which is a reasonable assump- 
tion, then it will make sense to keep these 0(S 2 ) hidden 
momentum terms in a pole-dipole approximation. This 
means that the treatments in [24] . studying the solu- 
tion given by Dixon- Tulczyjew condition without making 
P a « mll a , and in |40j . comparing with the condition 
pa || jja ^ ne difference between them being 0(S 2 )); or 
the work in jlOOj comparing Mathisson non-helical solu- 
tion with Dixon- Tulczyjew, are not a waste of time and 
make sense indeed. 

In the electromagnetic case the situation is different. 
As explained above, in general Ph idEM 7^ 0, so that 
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to dipole order one does not have P a ~ mU a (nor 
that Mathisson's non-helical solution is approximated by 
Dixon's Tulczyjew). As for P£ idl , it cannot also be ne- 
glected by arguments in the line of the ones above, and 
should in principle be kept, because when the particle has 
a net charge, the Lorentz force, first term, of Eq. (391, 



will yield a contribution to the acceleration independent 
of /t (or S), hence S Xjjd will be of first order in S. If, on 
the other hand, one considers the magnetic moment jl to 
arise from an uncharged body (for instance, the current 
loop depicted in Fig. [3^,), then P^ 
we focus on the dominant first term of (|46[)): 



hidi isjypically (again 



\S x v a\ 



-\{E u ) i:j \S 2 ~ mav\(E u ) l3 \R 2 ~ vma, 



(C7) 



To get an estimate of the quadrupole contribution to the 
hidden momentum, we contract Eq. ( 109 1 with Up, yield- 
ing 2m} a pt j.F^^ p Ui3 which involves terms of the type 

{h u )«J 5 ^2E Sl ~ am\{E u ) ZJ \R 2 ~ \S x v a\/v (C8) 



^jup~i j g Qiegneoi m ( 112 ); in this case we base our 
estimate in the moments of the spatial current, i.e., 
(h u )2J afSa ) where we used (h u )*J Sia ~ (iR ~ amR 2 . 
Hence, for this special case (the one considered in 



Sec. Ill), arguments analogous to the given above for the 



gravitational case apply. 



3. Impact of the Dixon- Tulczyjew condition on the 
gravito-electromagnetic analogies 

a. Analogy based on tidal tensors 

If one uses Dixon- Tulczyjew condition, then = 
e^ TX (S p ) T Px/M (with Pirani's condition, it would be 
g»u = € ht*S t U x ), and therefore Eq. (gib) of Table g 
would read instead: 



DP° 
~dr~ 



(C9) 



where H a p = -kR a ^ iv U v P» jM . Thus the force is still of 
course given by a contraction of a rank 2 tensor R a p with 
the spin 4- vector S a ; however, this new "tidal tensor" 
does not coincide with the magnetic part of the Riemann 
tensor (H u ) a p — -kRa^^u^u" as measured by some ob- 
server u a , because R a 7 results from a contraction with 
two different vectors (P a and U a ). The analogy in Table 
1 will no longer be exact, since by taking the trace and an- 
tisymmetric parts of Ti. a fi , one does not obtain equations 
of the type |lj2b)-([l]3b). The equations one will obtain 
will no longer be projections of Einstein field equations, 
and they will no longer be analogous to their gravito- 
electric counterparts Eqs. (3b)-(4b) of Table 1 of [57], nor 
to their electromagnetic counterparts (3a)-(6a) of [57] . 
since no longer will they be equations involving only tidal 
tensors and sources. Also, the tensorial structure of H a p 



(unlike H Qj g) is not similar to its gravito-electric counter- 
part, i.e., it is not spatial, nor it has to be symmetric in 
vacuum. 

The reason for all this being, as stated above, the fact 



that R a p (i-e., Papapetrou Eq. (C9) above) involves both 
U a and P a , so one cannot save the (exact) analogy by, 
e.g., projecting in terms of P a /M instead of U a . Let us 
compute explicitly the deviation from the analogue form. 
Firstly we rewrite Eq. ( C4 ) as: 



U a = U^ + UZ; U^ = ^P a ; UZ = (h p ) a p pP 



(CIO) 

where and U" are, respectively, the projections of U a 
parallel and orthogonal to P a . We may then write H a /3 
in the two alternative forms: 



m 1 

"M Waf3+ M 

■-^(U^^ + ^+R^P^; 



note that 



M 2 



hid 



(Cll) 
(C12) 

(C13) 



If p a P ^O. UI'm the second term of (|C4j); it is useful to 
write it in terms of tidal tensors: 



M 2 M 2 + (¥ p )^(S p ) x (S p ) a 



(C14) 



where (H p ) a7 = *R afil& P !i P & j M 2 and (F p ) a7 = 
*i?* Q/ 3 7 5 pPp 5 /M 2 denote, respectively, the gravito- 
magnetic tidal tensor and the "F" tensor [3H 1 101 j mea- 
sured by an observer of 4- velocity u a — P a /M. Thus 
we see that indeed the analogy in Table |T] is not exact 
if one uses the spin condition S a/3 P/3 = 0; the differing 
terms being 0(S 2 ) if one considers the purel y gravita- 
tional case (F afi = 0), as shown by Eqs (jCll|)-( |C14| ), or 
O(S) if one considers the gravitational force Fq in the 
general case that the two interactions are present, and 
assuming /x' 
due to Pi 



a/3 



crS a P, since in this case P£ d ~ O(S) 



hidEM ' 



pa 

(70) of 



cf . Eq. ( 25 ) (for an explicit expression for 



hid under Dixon- Tulczyjew condition, see Eqs. (35) and 



b. Analogy based on GEM inertial fields 



The expressions for the momentum ( 29 1 , force [3 El] , 
and spin evolution ( 24 1 , obtained within the gravito- 
electromagnetic analogy based in the 1+3 formalism, dis- 
cussed in [31], also rely on the condition S a ^Up — 0. 

Hidden momentum. - If the reference worldline is 
the one given by condition S a @ Pr = 0, the momentum 
no longer takes the forms (27)- (29 1. The inertial hid- 



den momentum P^ i(il 



U p lW a JdT in (25]) is no longer 
S a "ap, thus no longer exhibits an exact formal analogy 
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with its electromagnetic counterpart P^ [dEM = —^ a ^Ep. 
The analogy can however be recovered at an approximate 
level [IS]. Using P a = mU a + Pf lld and S^Pp =0wc 
can write: 



1 

m 



dm 
~dr 



(C15) 



i.e., P h Q idI = S af3 ap + 0(S n ), being n = 3 if F Qf9 = (and 
there are no other external forces), as follows from Eqs. 
| C13| -( |Cl4, o r n = 2 if F a ^ ^ 0, as follows from Eqs. 
(25) and (C13). Hence, at least to first order in S, the 



analogy holds. 

Force on the gyroscope. — The force Eq. (5) of [5], 
when it holds, is the same as ([T] lb) , only written in a 
different formalism; hence the discussion in Sec. |C 3 a| ap- 
plies. 

Spin precession. — As discussed above, only if Pirani's 
condition holds, a gyroscope (in the absence of electro- 
magnetic field) undergoes Fermi- Walker transport. If one 
adopts the condition S a "Pp = 0, we have (cf. Eq. (7.11) 
of [22]): 



DSa =SM" DU ° 



dr 



dr 



1 DP" 7 
M* ba dr 



(C16) 



(U a = P a /M) which is a different transport law, pre- 
serving the orthogonality to P a (instead of U a ), dubbed 
in [55] ("M transport"). Note that there is no conflict 
with Eq. (23) because the transport is along a different 
worldline. In the CM frame we have 



DS l 

~dT 



^ S °^ 



Thus, in the comoving (congruence adapted: SI 



uj, see 



B 



frame e^, the expression corresponding to (24) (with 
0) reads with this spin condition 



dS l ct ft 1 X a DP ° 



(C17) 



showing that the term breaking the analogy with the elec- 
tromagnetic counterpart comes from the hidden momen- 
tum. In the absence of electromagnetic field and ot her ex - 
terna l forces and torques, P£ ld ~ 0(S 2 ), cf. Eqs. (C13)- 
(C14), and thus the extra term is of order 0(S^. If 
the electromagnetic field is present (and fj, a ^ — <rS a P), 
Pu,^ ~ O(S) and the extra term is of order 0(S 2 ), cf. 



Eq. (25) 



DP a . 



dr 



dM 

~dV 



-W a {S p )^A a 



(U a = P a /M; M 2 = P a P a ), t he la st equality being 
obtained using Eqs. (C12) and (C14). M (not m) is 
thus the conserved quantity, and has the interpretation 
of energy of the body relative to the P l — frame. 

In the electromagnetic case, from Eq. ([8]), we have, for 
a magnetic dipole: 



dM 

~~dT 



DP a 
dr 



which is non zero if the derivative of F a P is not zero 
along P a . Therefore, in this respect, the conclusions are 
not very different comparing with Pirani condition: with 
the latter m is conserved for a gyroscope in a gravita- 
tional field, and in the electromagnetic counterpart its 
variation reflects the phenomenon of electromagnetic in- 
duction; with Dixon's condition, we have that M is con- 
served in the gravitational case, and that it varies in the 
electromagnetic case, again due do induction effects (only 
this time they are measured relative to the P l = frame) . 



Appendix D: Electromagnetic momentum and 
hidden momentum 

Let {T tot ) a0 = (T mattcr ) Q ' 3 + 6^ be the total energy 
momentum of the system, with (T mat tcr) Q ^ being the 
energy-momentum of the matter and 9 Q/3 the electro- 
magnetic stress-energy tensor: 



9 



a/3 _ 



47T 



(F tot ) Q7 (F tot )^ - \g afi (F tot )^(F tot r» 



(Tt t) al3 obeys the conservation equations: (T to t) a/3 fl 



(Dl) 



0; and 9 



af3 



-(Ftot) afi ip (cf. Eq. (7.36) of [7D|, 
Eq. (12.118) of 63 ). Therefore 

{T tot ) a % = => (T matter ) Q ^ = -Q a % = F°#3 P (D2) 

The electromagnetic stress-energy tensor splits in three 
parts (Eq. (36) of [3S]): 



Q aP = (9 cxt )^ + (O x ) afj + (e sclf ) afj 



where 



Conserved quantities with Dixon- Tulczyjew 
condition and electromagnetic induction 



(9 0Xt ) 



a/3 



4-7T 



Let us start by the purely gravitational case. With 
Dixon condition, it is in the frame of zero 3-momcntum 
(P l = 0) that the time projection of the force on a gyro- 
scope vanishes: 



(9 



l-[F^(F p )^ + (F p )^F^ 
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(e solf ) 



47T 



g^(F p ), u (F p r 



Since we are considering test particles, strictly speak- 
ing, one was supposed to neglect the electromagnetic 
field (F p ) ai produced by the test particle. However that 
would be inconsistent with taking into account its dipolc 
structure, as the cross term (0 X ) Q ^ above yields, as we 
shall see, a contribution to the momentum of the type 
[2 x E: this cannot be consistently neglected in an ap- 
proximation where couplings of the type B ■ p, are taken 
into account. This is in agreement with the treatment in 
Sec. VB of [H]. 

The part of O"' 9 which is to be neglected in this ap- 
proximation is the term (Q S eif) a ^ , which is quadratic in 
the test particle's field, and is related with the electro- 
magnetic self force [3^]. Hence, in this paper: 

e a ^ (e oxt )^ + (e x ) Q0 . 

If the only currents present are from the internal struc- 



ture of the test particle: j a = (j f 



then (P oxt ) 



Ct0 



0. 



and therefore (Q e xt) a % — 0. That is: the stress energy 
tensor of the external electromagnetic field is indepen- 
dently conserved (cf. [39] pp 7-8). Also, if the only mat- 
ter present is the test particle's (herein we are interested 
in a system consisting solely on a test particle plus an 
electromagnetic field): (T matt0 r) Q/J = T afi (T a P denotes 
the energy momentum of the test particle); then it follows 
from dD2| 



na/3 



ft"/ 3 



-G 



(6> 



NQ/3 



(D3) 



which is Dixon's equation (1.1) of [22] (with a relative 
minus due to the different signature) , Eq. (5) of (102] , or 
Eq. (10) of [25] (with /„ = 0). 

By obeying (D3), T a/3 encodes the hidden momenta 



of the particle. In order to see that, consider for sim- 
plicity flat spacetime and Lorentz coordinates. In this 
case we have T a0 = F a ^j . Noting that T a i = 



(T a P x "') i p - T Q ^2 7 , we have, from Eq. |4j) (work in the 
frame comoving with the test particle): 



T 0i d 3 x 



The second equality holds in stationary conditions (so 
that / j l x k d 3 x = J j^x k ^d 3 x = fi kl by conservation of 
charge) which is the problem at hand. This equals P mdEM 
[last term of Eqs. ( [25] ) or (27)], as expected, since for a 
stationary configuration U* = a 1 — and therefore 





Figure 7: A magnetic dipole subject to the Coulomb field of 
a spherical charge. Although the dipole is not moving, it has 
a non-vanishing hidden momentum Phid = PhidEM — ft x E 
which exactly cancels out the electromagnetic field momen- 
tum Pem = —ft x E, as one concludes from the conservation 
equations = 0. 



pl = ^hid = p L}EM, cf - E qs. (26), (28). In vec- 



tor notation, Eq. (D4) reads P = fx X E; this is pre- 



cisely minus the electromagnetic momentum generated 
by the presence of a magnetic dipole in an external elec- 
tromagnetic field, as is easy to show. In order to see that, 
first note that such momentum reduces to the cross term 
P% ee J s(u) (e x ) a/3 dS,9, as (@ se w) a0 is to be neglected in 
accordance with discussion above. Also one notes that, in 
the particle's frame (U l = 0), the electric field produced 
by the magnetic dipole is zero; therefore, the spatial mo- 
mentum is given by (see e.g. [53] p. 286, [5"9"]): 



Px = 



1 

in 



E x B d ip ieci 3 2; = E X p., (D5) 



= [ {T m x l ) fi d 3 x + [ {T 0k x l ) ik d 3 x - [ F Qk j k x l d 3 x momentum, since it follows from equation (T, 

J J J fViof fni~ n c+nli n*n nvn fi nil i~n i in m P l 



i.e., the volume integral of the Poynting vector p = 
-Exi3dipoie/47r. In order to convince ourselves that the ex- 
istence of the hidden momentum Phid = — Px is a require- 
ment of the conservation equations, consider the simple 
example depicted in Fig. [7] a magnetic dipole in the 
Coulomb field of a spherical charge. Although the con- 
figuration is stationary (from the viewpoint of the static 
observers), the dipole possesses a non- vanishing spatial 





that, for a stationary configuration, ± tot 



tot J . Q 

PL, = 0. Since 



{T w x l ) fi d 3 x+q)T" k x l n k dS- I F >k j k x l d 3 x. 



For a stationary configuration, the first term vanishes; 
the surface integral in second term also vanishes in a sur- 
face enclosing the test particle. Hence, neglecting the 
variation of F a/3 along the dipole, and using definitions 
^ and (10), we have 



matteiiXi EM ' uilcli 1 matter — ^EM' ^ n thlS Case, 



+P£ M , then PI 



pi j?0k 



J Jk x l d 3 x = -P°V^ = E^ttp . (D4) 



■pa — r>a 
^tpt — r m. 

P^ given by ( D5 ) is the total electromagnetic momentum: 
P^ = P EM - Thus we know that the matter must possess 

a hidden momentum P ma tter = — Pem — M x E; and this 
momentum must be carried by the current loop, since the 
static charge can carry no momentum. 

Thus, we succeeded in showing that in stationary 
conditions, the conservation equations imply — Pg M = 
P]^ dEM . For more general conditions we cannot follow the 
same procedure to obtain such equality (see also e.g. [29"] . 
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p. 20). But from Eq. (271 we can see that indeed Dixon's 



equations (which come precisely from the conservation 
Eq. (D3)) imply that the equality always holds. 



1. Physical model of electromagnetic hidden 
momentum. Inertial analogue. 

In this work we dubbed the contribution P{? idEM to 
Pj°j d as hidden momentum "of electromagnetic origin", 
because it only exists when an electromagnetic field is 
present; as explained above, for stationary configura- 
tions it equals minus the momentum of the field. It is 
however purely mechanical in nature [29] [62J I103j , and 
can be explained by different models (see [103] p. 520, 
[5"9l R)2"] ). all yielding the same results. For our purposes 
it suffices to consider the more simple and intuitive one, 
given in Fig. 9 of [621 : a square current loop, under 
the action of an electric field (which, to this accuracy, 
is taken constant along the body). The electric field 
pointing upwards, accelerates the charges on the right 
segment, and decelerates the ones on the left segment; 
thereby causing an increase in the concentration of pos- 
itive charges of in bottom segment (not the upper seg- 
ment, as one might naively expect!). That can be seen 
recalling that the flux of particles nv (i.e., the number 
of particles per unit time) is constant throughout the 
loop (otherwise, charge would be piling up somewhere), 
i.e., n bo tVbot = «top«to P ; therefore, if v top > v hot , then 
n bot > n top- n is equal for the left and right segments. 
Since the relativistic density of momentum is p = nvym, 
non-relativistically (7 w 1), the total momentum of the 
loop would be zero, as the momenta from the top and bot- 
tom, and the left and right segments, would cancel out. 
In the relativistic case, however, due to 7 , the top seg- 
ment has larger momentum than the bottom one; which 
means that the loop, in spite of being at rest, has a non- 
vanishing momentum [6"2"I ll03| P hid = m(^ top -^ hot )H /q 
(note that II /q = N top v top = N hot v hot , with N = num- 
ber of charges in a segment). That is, hidden momentum. 
Since m(7t op — 7bot) = qEw (i.e., the gain in kinetic en- 
ergy going from the bottom to the top is the work done 
by the electric field), and /i = Iwl, we obtain, for the 
setup in Fig. 9 of [62], PhidEM = fix E. 

Inertial analogue. — Replacing the charges in the loop 



of Fig. 9 of [62 by masses, and the electric field E by the 
gravitoelectric field G = —a, as measured in the frame co- 
moving with the CM (i.e., a is the CM acceleration), one 
obtains analogously a hidden momentum S x G, in this 
case originated by the inertial force G = —a. This obser- 
vation makes natural the similarity between the "electro- 
magnetic" hidden momentum fixE (i.e., /j, ea Eg) and the 
term axS (i.e., S a6 ag) coming from the Mathisson-Pirani 
condition. However it must be noted that this model can 
only be used to interpret the hidden momentum term 
given by Mathisson-Pirani condition; i.e., when the CM 
is computed in its rest frame. For any other spin con- 
dition/choice of observers u a relative to which the CM 
is evaluated, P" id ^ S a6 ag (in general it is not related 
to the acceleration of the body); also, consistently, the 
CM Xq M (u) (thus the apparatus in Fig. 9 of [52]) is not 
at rest in the frame u % = 0. This is the case even if one 
computes the CM relative to an observer u a momentarily 
comoving with the proper CM Xq M (U), but with differ- 
ent acceleration: the CM is measured to be at the same 
position Xq M (u) — Xq M (U), but their velocities are dif- 
ferent, as they depend on the acceleration, see Eq. (27) 
of [IS]. 

Hence, albeit the analogue model in Fig. may be used 
to interpret the term P{^ d j and the reason for its formal 
similarity with P^em, one must bear in mind that the 
fundamental, generic, explanation of the inertial hidden 
momentum is the one given in Sec. |IID| above, and Sec. 
VI of 125 

: x cm( u ) depends on u a ; if u a changes along 
the curve, the CM velocity changes in a non-trivial way, 
causing it to move in the P % = frame; conversely, in 
the CM frame U l = 0, P l is not zero, and that is hidden 
momentum. 
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